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Supersonic Wing-Body Lift 


GEORGE MORIKAWA* 
Hughes Aircraft Company 


ABSTRACT 


The representation of lift for a fairly general family of wing-body 
afterbody) in a linearized supersonic flow 
entire range of the The 
are obtained simply using the most ele- 


configurations (with no 
is considered over the 
bracketing limits of lift 
mentary considerations 
is recognized. By the proper choice of the independent param- 
eters and the proper representation of the lift, reasonable ap- 
proximations for the behavior of lift for intermediate values of the 


parameters. 


when the interpretation of these limits 


parameters are obtained. 

Further generalization to other wing plan forms with straight 
edges and other body shapes and the representation of other 
aerodynamic variables such as drag, moment, and center of 
pressure will be evident but not considered in detail 


NOTATION 


radius of body 

ratio of tangent of tip trailing-edge 
w, to tangent of Mach angle yu 

aspect-ratio of “exposed wing”’ 


angle 


root chord 

lift-curve slope for wing-body combination 
based on “exposed wing”’ area 

lift-curve slope for body portion of combina- 
tion based on “‘exposed wing”’ area 

lift-curve slope for wing portion of combina- 
tion based on ‘‘exposed wing”’ area 

lift-curve slope for ‘‘exposed wing’’ alone 

pressure coefficient 

pressure coefficient along centerline of fi 
delta wing 

complete elliptic integral of second kind 

Cia/Cra, total lift ratio 

C1 ms B) 

Cre’? Cra, Wing lift ratio 

ratio of 
to tangent of Mach angle u 

free-stream Mach Number 

free-stream dynamic pressure 

semispan of wing-body combination me: 


Cre, body lift ratio 


tangent of leading-edge angle 


ured from body centerline 
Presented at the Aerodynamics Session, Annual Summer 
Meeting, I.A.S., Los Angeles, July 12-14, 1950. 
* Missile Aerodynamics Section, Research and Development 


Laboratories. 


217 


projected area onto plane of wing (x-z 
plane) upon which body lift acts 

“exposed wing”’ area 

spanwise coordinate 

vertical coordinate normal to plane of wing 

chordwise in free-stream direc- 


tion coincident with body centerline (cf. 


coordinate 


Fig. 1) 

angle of attack 

geometrical angle of attack of both body and 
wing of combination 

V M?—1=ctnyg 

a/s, diameter-span ratio 

Mach angle measured from free-stream di- 
rection 

leading-edge 
stream direction 


angle measured from free- 
tip trailing-edge angle measured from free- 
stream direction (cf. Fig. 1) 


INTRODUCTION 


N THE AERODYNAMIC DESIGN of vehicles for super- 
sonic flight, it is essential that the bracketing limits 
lift, 


I 


of the important aerodynamic variables 
-as functions of 


i.€., 
drag, moment, and center of pressure 
the design parameters be known. Frequently, an 
“in the large’ analysis of this nature may be initiated 
using the most elementary considerations with the 
proper choice of the independent parameters and the 
proper representation of the aerodynamic variables, 
and from this study some insight is obtained for the 
more difficult analysis over the entire range of the param- 
eters. The main consideration in the present paper 
is the representation of lift for a fairly general family of 
stationary linearized 


] 
a 


wing-body configurations in < 
supersonic flow. Further generalization to other wing 
plan forms and body shapes and the representation of 
other aerodynamic variables will be evident but not 
considered in detail. 

The considerations of this paper are a natural and 
direct consequence of previous related work, notably 
by Beskin,! Lagerstrom and Van Dyke,’ and Spreiter.* 
Grateful acknowledgment is made of helpful discus- 
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(c) CLIPPED-DELTA WING-BODY, Oa,/m#!I 


FIGURE |. 
FAMILY OF WING-BODY CONFIGURATIONS WITH NO AFTERBODY 


sions with A. E. Puckett, T. F. Coleman, and J. L. 
Fredrick, of Hughes Aircraft Company; and H. J. 
Stewart, P. A. Lagerstrom, and C. R. DePrima, of the 
California Institute of Technology. 
REPRESENTATION OF INTEGRATED VALUES: LIFT, 
DRAG, MOMENT, AND CENTER OF PRESSURE 


The principal assumption made in the following dis- 
cussion, aside from the linearization of the equations of 
motion, is that the flow field immediately upstream of 
the wing is cylindrical (independent of the flow direc- 
tion)—that is, the influence of the nose of the body is 
assumed to be negligible. In particular, an infinitely 
long circular cylinder represents the body.''* 4 Only 
the lifting problem where both body and wing (flat 
plate) are at an angle of attack a» with respect to the 
free-stream direction is considered. It is remarked 
that this problem may be replaced by an equivalent 
problem*: * where the body angle of attack a = 0 and 
the wing is regarded as twisted with 


a@ = ao[1 + (a?/x?)] 


The coordinates are shown in Fig. 1. The main dis- 
cussion is confined to an important family of wing 
plan forms with straight edges and trailing edge normal 
to the free-stream direction; the body ends at the 
trailing edge. 

Referring to Fig. 1 and using Doris Cohen’s nota- 
tion® for wings, the basic parameters for this family of 
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“‘clipped-delta” wing-body combinations are \, m, and 
a,, where \ = a/s is the diameter-span ratio,t m is the 
ratio of the tangents of the leading edge angle to the 
Mach angle, and a, is the ratio of the tangents of the 
tip trailing-edge angle to the Mach angle. Then 


m = tan wo/tan uw = 6B tan wo 
a, = tan w)/tan uw = £6 tan a, 
a,/m = tan w,/tan wo 


where 8 = VM? — 1. 
parameters are taken to be X, a,, and a,/m, where 
0<A<1landO <a,/m <1. Then, a,/m = | de. 
fines the particular family of delta wing-body combina- 
tions (Fig. la), and a,/m = 0 indicates the family of 
rectangular wing-body combinations (Fig. lb). Clearly, 
if the trailing edge is swept (either back or forward), 
an additional parameter is introduced 


For convenience, the three 


namely, the 
ratio of the tangents of the trailing-edge angle to the 
Mach angle. 

In considering a representation for lift of wing-body 
configurations, it is reasonable to regard the ‘exposed 
wing”’ as the principal lift-carrying agent and any addi- 
tional lift (plus or minus) contributed by the body 
adjacent to the wing to be purely supplementary. To 
regard the “exposed wing”’ alone as the basic wing is 
the natural consequence of attending to the limiting 
cases of X = O anddA = 1. 
on the wing for the wing-body configuration and the 
“exposed wing’ alone are similar—i.e., for \ = 0 the 
adjacent body becomes negligible; for \ = 1, the ad- 
jacent body acts as a vertical plane reflector identical 
to the plane of symmetry of the ‘‘exposed wing”’ alone. 
Thus, the general point of view is taken that, given the 
integrated values for wings, the integrated values for 
wing-body configurations are directly obtainable— 
exactly for limiting values of the parameters and to good 
approximation for intermediate cases. A remark is 
appropriate here concerning Spreiter’s representation’ 
for lift in the limiting case exhibited by slender body 
theory, where the basic wing is taken to be a fictitious 
wing that passes through the body. For a study over 
the entire range of the design parameters, such a repre- 
sentation suppresses and, hence, obscures the behavior 
of wing-body configurations for values of \ near 1— 
i.e., for low span-diameter ratio configurations. 

Clearly, the preceding discussion holds for the repre- 


The boundary conditions 


sentation of other integrated values such as drag, 
moment, and center of pressure. In the following sec- 
tion the lift representation for delta wing-body con- 
figurations with no afterbody is discussed in detail. 


LIFT REPRESENTATION FOR DELTA WING-BoDY 
CONFIGURATIONS 


As proposed in the previous section, the lift for the 
delta wing-body configurations is represented by the 


t It should be noted that \ = 0 defines the limiting case of all 
wing and \ = 1 defines the limiting case of all body. 
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ratio of the total liftt of the combination to the lift of 
the “exposed wing”’ alone, 


K = Ge Cra, = fa, m) (la) 


4 function of the two parameters{ \ and m. For the 
a + 


basic “exposed wing”’ alone, 


on m hte ¢ 
Cia, = yy (forO <m <1 
“ 8B E( Vi- m?) (1b) 
= 4/8 (for m > 1) 
where E( V1— m?) is the complete elliptic integral of 


the second kind. It is instructive to separate the total 
lift C,, into that portion acting on the wing C,,°" and 
that portion acting on the body C,,‘"’, where 


. = ( 2a) 


Cr. si + Cre‘ e 


and 


K = KW + K® (2b) 


t To conform to engineering usage, this is expressed as the 
ratio of lift curve slopes based on the same (‘‘exposed wing’’) 
area. 

t For consistency with the more general clipped-delta wing- 
body family, the parameter a; instead of m should be used for the 
delta wing-body configurations, but m is the more commonly 


used parameter. 


(1) Limiting Total Lift 


WING-BODY 
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LIFT 


Note that this is an actual experimentally measurable 
breakdown of lift and generally includes interaction 
and tip effects. Thus, the wing lift C,,‘" generally 
includes lift due to (A) the angle of attack ao of the 
wing, (B) the upwash ao(a?/x*) of the body upstream 
of the wing, (C) the reflection off the body and the 
plane of symmetry (the y-z plane) of the disturbances 
due to (A) and (B), and (D) the reflection of wing tip 
disturbances 

The lift for limiting values of the parameter m may 
be obtained simply when the interpretation of these 
limits is recognized: (A) m = O implies either that 
the free-stream Mach Number / = 1 (wo # 2/2) or 
the leading-edge angle w = 0 (WM # ©); this is 
clearly the limiting case of “slender body theory.” 


(B) m = © implies either that 1 = © (wo) # O) or 
w) = 7/2 (M #1); thisis the limiting case** of “‘strip 
theory.’ In both limiting cases, the analysis is simpli- 
fied by the reduction of the original three-dimensional 


problem to two dimensions. From Spreiter’s results® 
and some elementary integration and algebraic manip- 
ulations, the results may be tabulated, letting a/s = 
\ for convenience. 

** Note that the body carries no lift form = © (e.g., for fixed 
A, the chord approaches zero as wo approaches 7/2); and 


K= KWM= assw)| f (a/ay) dS 
OW 


m=0: K = (1+)? (3a) 
m= oe: K = [1/(1 — d)2][(1 — A?) + 2A? In A] (3b) 
(2) Limiting Wing Lift 
; 2 1 ee I d , 
m=0: KW = (1 + d?*)? sin7! j =}—-(1-A*%) -- (4a) 
(1—A)? Lar Vi+ ™ 
m= o: K® = [1/(1 — d)2][1 — A2) + 2A2INnA] HK (4b) 
{3) Limiting Body Lift 
K 1 + A‘) sin~! : 
m = 0: oO = ( sin =} — 2d? sin— ] ' A(1 — d? (5: 
r(1 — )? ‘1+. ite é 
m=o: K® =O (5b) 
The lift for m = © can also be obtained directly from and 
reference 2 for this iguration. Eqs. (é i , , 
ce or this configuration. Eqs. (3), (4), and dK(0, 0)/OA = dK(0, ©)/dd = 2 (6b) 


(5) are plotted in Fig. 2. 


Some preliminary observations of interest may now 
be made. (A) The total lift K(A, m) for the limiting 
cases of m = 0 and m= = and also the slope 0K(A,m) + 
Od are identical at \ = 0 


K(0, 0) = K(0O, ~) = 1 


i.e., 


(6a) 


and the difference between A(A, 0) and AK(A, ©) in- 
creases monotonically with \ to a maximum at A = l, 
at which 


K(1, 0) — A(1, ~) = 2 (7a) 


The slopes at A = 1 are 
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FIGURE 2. DELTA WING~BODY LIFT WITH NO AFTERBODY, o,/m=1.0 chore 
OK(1, 0)/OA = 4 (7b) and the difference [K“ (A, 0) — K‘(X, ©)] increases 
- monotonically with \ to a maximum at A = 1, at which 
a 
‘ a K‘®(1, 0) — K®(1, ©)] = 2 (12a) 
OK(1, ~)/OA = 2/3 (7c) | ) 
(B) The difference between the wing lift A‘ (A, 0) and The slopes at A = | are 
K)(\, ©) for the limiting values of m is small—i.e., oK(1, 0)/OA = 2[1 + (4/3x)] (12b) ‘ 
Since 
[KW (A, ©) — KOA, 0)] < 0.15K (A, 0) (Sa) OK'(1, ~)/ = 0 (12c) 
for all \ and is zero at) = Oandd = 1, Thus, it is clear that the large difference between the ” 
; ‘ : an 
[K™ (0, ©) — K (0, 0)] = total lifts, [A(A, 0) — K(A, ~)], for small span-diameter ” 
bias sind ratio configurations—i.e., \ near 1—is due primarily to 
_ Ww = . . ° r 
[Km(1, ©) — KOMI, 0)] = 0 (8b) the large differences between the body lifts [K‘® (A, 0) — 
The slopes for \ = 0 and \ = Lare K®(,, ©)]. Initially then, it appears imperative to 
. )/ , ‘ , obtain at least approximate expressions for body lift then 
Ww) = * — (° Di e ° Ps 
oK™(0, 0)/OA = 2[1 — (2/m)] (9a) for intermediate values of m,0 < m < @. 
OK (0, ©)/OA = 2 (9b) F K() 
Approximate Body Lift for 0 < m< = 
and “ ; . se a 
From previous studies of pressure distributions and 
OK (1, 0)/OA = 2[1 — (4/3r)] (10a) lift? * 7 on wing-body configurations for intermediate 
ee srsonic Mach Numbers, it is evident that exact 
OKW(1, ©)/2% = 2/3 (10b) supersonic Mac ibers, it is evid bi 
calculations are, at the very least, difficult and prob- 
(C) The body lifts K (A, m) for m = 0 and m = ~ ably impractical. Even the approximate calculation 
are both zero at \ = 0, of local pressures for m > 1 (supersonic leading edge) 
K®(0, 0) = K®(0, ©) =0 (11a) is tedious, the difficulties becoming more apparent on 
the body and more so for 0 < m < 1 (subsonic leading ory 
the slopes being edge). However, for m > 1, the calculations made im 
: reference 4 for the zero leading-edge sweepback case 
o0K(0, 0)/0A = 4/x (11b) ee ee =e pan” 
and later unpublished calculations for m = 1 and Ve 


0K (0, ~)/OA = 0 (lle) indicate that a reasonable approximation for the pres- 
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sure on the body (not influenced by the wing tip) is the 
asymptotic pressure at a large distance downstream 
from the leading edge; this approximation as an aver- 
age value on the body appears reasonable even near the 
leading edge. The asymptotic pressure is identically 
the centerline pressure on the corresponding “exposed 


wing’ alone, 


oem 41 m - 
~ = - 7 sec~! m (13) 
a ™BVm? — 1 
on the lower surface. 
First Approximation for m > 1. 
pressure on the body is constant (for fixed M and m) and 
Then, 


Assume that the 


given by Eq. (13). 


2 ‘ i 
CJ" «= — | f ( ") dS 
Sw Sp\ (14a) 


= 2(C,,/a0)(Ss/Sw) 


where the integration is performed over an area Sz, 
which is projected onto the x-z plane, and Sw is the 
“exposed wing’ area. The body area Sz is bounded 
by (A) the projection onto the plane of the wing of the 
helices (lead angle equal to the Mach angle) emanating 
from the leading-edge-body junction, (B) the root 
chord, c, and (C) the base. For 


c T a l 
) < < Boor AS a | 
a 2 Ss 1 + (17/2)m 
ae eo z dz 
Sp/a*?*= 2 / [ — cos ( ) (14b) 
J 0 Ba a 
| 1— 2X _ fl -rA&\, 
— 93 ( — sin 
mr mx 
Since 
c/a = B[(1 — A)/mA] 
and the “‘exposed wing’’ area 
oo. efl—-dr\ BLA: 
Syw/a? = = ( (14c) 
a a m dA 


then the body lift 


i m 


r 
= sec" me ( ) x 
7 Vm? — 1 od. 
mx . 1—A7A 
[ — ( ) * sin ( )| (14d) 
1—A mx 


form > 1 and dy > 1/[1 + (2/2)m]. 


Kis) = 


l 
and for 0 < c/a < (1/2)8 le se 
[1 + (2/2)m] 


. ‘ 4 C,. ‘ pa Zz dz 
CC," /e* = — 1 — cos ae > 
Sw \ ao 0 Ba nBa/ a 


< ’ 
AS 1y 
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Similarly, for c/a > (x/2)8 or 0 SA < 1/[{l + 
(1/2)m], 


S,/a? = 28{[(1 — A)/m\d] — 1} (14e) 


and the body lift 


; | m A 
KK = sec”! m x 
— | 
— 


T Vm? - J 
| mr : 
| (, ) (14f) 


The curves represented by Eqs. (14d) and (14f) are 
shown in Fig. 2c for m = 1 and m = 2 and define the 
lower bound of the cross-hatched areas in the figure 
for these particular values of m. The slopes obtained 
from Eqs. (14d) and (14f) for’ = Oand A = 1 are 


OK‘ (0, m7) 4 m 


= sec! m 
On = Vm* — 1 
(for 1 <m < o) (15a) 
and 
OK ‘® (1, m) 2 | 
=— sec! m 
OX ar mM Vm? — |] 
(for 1 <m < @) (15b) 
Note from Eq. (15a) that, for \ = 0 and m = 1, the 


slope is 4/7, which is identical to that for m = 0 given 
by Eq. (11b), but the slope increases with m. Thus, 
in a small region near \ = 0, A‘ (Xd, m) for 1 <m < @ 
lies slightly outside the region bounded by m = 0 and 
m = o in Fig. 2c; this is attributable to the approxi- 
mation used. In order to have some measure for esti- 
mating the possible deviation of this approximation, a 
second approximation that overestimates the body lift, 
at least for 1/[1 + (7/2)m] < AX < 1 forO < c/a <s 
(1r/2)8] is given. 

Second Approximation for m > 1. 
average pressure on the body varies linearly from (C,, + 


Assume that the 


a) at the leading-edge-body junction to ('/2) (Cy, = 
a) at a distance (7/2)8a downstream and is constant 
(1/2)(C,,/ao) over the remainder of the body [cf. Eq. 


(13) ]. 


(=) - (=) ( Me z ) 
a ~ ao 1 Ba 
T 


(for0 < ~ < ~~ B) 
cf 


) 
(2 Sw) / f iG. a) aS 
SR 


Then, 


(16) 
2 


" @®» oa 
Ca 


(17a) 
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and 
: 8 m? r 1-—A 1 /l1—A 
Ks) = - Fea =sec™* m( ) 1( )f _ ( )| = 
Wet = 1 1—r/ \\ m 2r mr , 
E (5) Jn FE) +2 =>) 1) 
_ sin —-? (17 
«\ my mX . my nf (17) 


and for c/a > (7/2)B {or0 < XA < 1/[1 + (x/2)m]}, 


sec! m (, 


The curves represented by Eqs. (17b) and (17c) are 
shown in Fig. 2c form = 1 and m = 2 and define the 
upper bound of the cross-hatched areas for these par- 
ticular values of m. The slopes obtained from Eqs. 
(17b) and (17c) for \ = O and A = 1 are 


oK (0, m) 4 m 
- ON 2 


q 2 
K® = Pe) mM 
™ Vm? — | 


- sec~! m 
™ Vm? — 1 


(forl1 < m < ~) (18a) 


which is identical to that obtained by the first approxi- 
mation [cf., Eq. (15a)], and 


OK (1, m) 4 1 


—_ + 1 
= sec~! m 
Or 


3 m Vm? — | 


(for l1 <m< @) (18b) 


which is twice the slope obtained by the first approxi- 
mation [cf. Eq. (15b)]. 

Reasonable approximations for the subsonic leading- 
edge case 0 < m < | are difficult to obtain. However, 
for m near zero—i.e., either M near 1 or w) near 0— 
one expects® the spanwise pressure distribution over the 
body to resemble that for m = 0 except near the helical 
boundary—i.e., near the leading-edge-body junction. 
It appears reasonable then to assume that for m near 


: Xr 2 l 2 
K®=E (vi _ m)( -) ( ) x 
1 No 1—A 
2 : No do l 
(1 + Ao‘) sin , = — sin? | —— 
7 V1 +d? us V1+ 


The curves represented by Eqs. (19a) and (20) are 
shown in Fig. 2c for m = 0.05 and 0.10 and define the 
lower bound of the cross-hatched areas in the figure 
for these particular values of m. 

Second Approximation for 0 < m < 1.—A second ap- 
proximation may reasonably be that the linear shift is 
less than that defined by Eq. (19a) but is a mean dis- 


tance given by so/a,= [(s/a) — m] or 


do = A/(1 — md) (19b) 


Then the body lift AK“ is given by Eqs. (19b) and (20) 


(17e 


A WG) +2-5- 
—A \\ 2m 8 2 


zero the chordwise distribution of pressure may be rep. 
resented by that for m = 0 but with a prescribed linear 
shift downstream to indicate the formation of the zero 
pressure region between the x-y plane (the crossflow. 
plane at z = OQ) and the helical boundary. 


First Approximation for 0 < m < 1.—Assume that 
the linear shift is of length (7/2)8. Then the body 
lift C,,‘” may be obtained immediately from the slender 
body theory result for m = 0 by replacing s/a for m = 


0 by so/a = s/a — (r/2)m forO < m < lor 


d 


ho = a 
. 1 — (4 2)my’ 


(19a) 


for O<m <landaA < | 
1 + (7/2)m 


Note that this approximation is not valid for the chord 
length 


a 


c T l 
< — § orA > a 
2 1 + (2/2)m 
in fact, the approximation states that the body lift 
is zero for 0 < c/a < (r/2)8. Then, with Ao given by 
Eq. (19a), 


2 I 
Ao(1 — do” forvA < (20) 
= * T | c 1 + (x == | 


for\ < 1/(1 + m). 
equations are shown in Fig. 2c for m = 0.05 and 0.10 
and define the upper bound of the cross-hatched areas 
for the corresponding values of m. It is noted that 
for both approximations a definite maximum is found 
and that this family of curves for 0 < m < 1 resembles 
the family of curvesform > 1; the curves for 0 < m < 
1 were propitiously faired to K“(1, m) = 0 over the 
small region where Eq. (20) is invalid at and near A = 
1 

With the approximate representation of the body 
lift K®, there remains the representation of wing lift 


The curves represented by these 
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SUPERSONIC 


kK for 0 <_m < © before the total lift A for all values 


of m can be discussed. 


(5) Approximate Wing Lift for0 <m< 


Meaningful approximations for the wing lift are dif- 
fcult to obtain since the difference [AK (A, ©) — 
K(x, 0)] is uniformly small for all A, 0 < A < 1. 
However, graphical integrations of the pressure distri- 
bution (unpublished) on the wing calculated by the 
methods of reference 4 yield the curve shown in Fig. 
% for m = 1, A > 1/(1 + wm). Calculations for 
m = 3 give a nearly identical curve. 

Some additional observations may be made of the 
variation of OK‘? /OX with m at the end points \ = 0 
and A = 1. AtA = 1—i.e., at the limiting case of 
the effective angle of attack of the wing with 


all body 
the equivalent twist due to body upwash is linear, 

a/ao = 2[2 — (x/a)], [(x/a) -1]<1 (21) 
and the body acts as a vertical plane reflector.t For 


the supersonic leading edge case m > 1, the wing lift is 
simply obtained by Lagerstrom’s theorem? on the 
preservation of lift, which has also been given inde- 


pendently by Ward :° 


4 4 | z 
= — Wave = - = / [ a dS 
B Sw J Sw 


CL 
Then, 
Kw = 211 — {(1 — A)/3A]} (22a) 
and 
OK (1, m)/OrA = 2/3 (22b) 


which is identically that for m = © [ef. Eq. (10b)]. 
Thus, for m > 1, the slope at \ = 1 is independent of 
m. 

The slope at A = 1 for the subsonic leading-edge 
case 0 < m < 1 may be obtained by a somewhat 
lengthier calculation. However, an approximate cal- 
culation of the wing lift slope 0K‘ (1, m)/OA is made 
using Martha Graham's modified strip method® and 
assuming the twist distribution to be linear [exact at 
\= 1, cf. Eq. (21)]. Then, by a simple integration for 


\ near 1, 


i’ (w)/ 
K™(\, m) = 


2 E(V1— m’) fl —-d ' 
( (23a) 
oT Vm \ A 

and, at\ = 1, 

0K (1, m)/v =(4/3m)_LE( V1 — m?2)/ Vm] (23d) 


The strip approximation is exact for m = 1 and is in- 


t Actually, this last condition is sufficient and implies that, 
for any m > 1, the lift K™ and lift slope 0K‘")/ovA at AX = 1 


IS precisely that for m = ©. 
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FIGURE 3. WING LIFT SLOPE WITH RESPECT 


TO X} FOR DELTA WING-BODY AT A=! 


valid at m = 0. Fig. 3 shows the variation of the wing 
lift slope given by Eq. (23b) at \ = 1; the dashed por- 
tion is the interpolation of the expected exact curve to 
the correct value for m = 0 [ef. Eq. (10a) ]. 

At \ = 0, the wing lift slope 0K‘”(0, 1)/0A may be 
obtained for the sonic leading-edge case m = 1. As- 
suming that the total lift slope 0A(0, m)/Od is inde- 
pendent of m—i.e., identical to that for m = O and 


m = © (cf. Eq. (6b) ]—and, since the body lift slope 


OK‘®)(0, 1)/0A = OK‘(0, 0)/OA 


the wing lift slope for \ = 0, m = 1 is necessary identi- 
cal to the wing slope for \ = 0, m = O given by Eq. 
(9a). This implies that for \ = 0,0 < m < 1 (sub- 
sonic leading edge) the wing lift slope is independent 
of m. Now with these approximate representations 
of wing and body lifts for intermediate values of m, 
0 <m < ,a basis has been established for represent- 
ing the total lift for allA,0 < A < landallm,0 < m 
< o, The calculations indicate that (A) all the total 
lift curves lie in the region bounded by the total lift 
curves for m = 0 and m = o; (B) the total lift 
K(0, m) = 1, independent of m; and (C) the total lift 
K(1, m) = 2at\ = 1forO0 <m< ~. The observa- 
tion that there does not exist continuous dependence 
of the lift on the parameters—i.e., dependence of body 
lift on m at X = 1—is not startling, since for m = 0 
the wave equation is reduced [since the coefficient 
(M? — 1) of the second derivative with respect to the 
flow direction is zero] to the two-dimensional Laplace's 


equation. Sometimes this is called, mathematically 
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speaking, the boundary-layer phenomenon, stretching either = © or w, = 7/2 (this is the case given by 


the concepts a little. 


(6) Remarks on the Analysis 


Evidently in the foregoing analysis, two underlying 
assumptions are made concerning the behavior of the 


lift: (A) the lift for all finite m is bounded by the lift 
form = Oandm = o—i.e., all the lift curves lie be- 
tween the lift curves for m = 0 and m = ©—and (B) 


for fixed the variation of lift is monotonic in m. 
Mathematical verification appears difficult, but there 
exist strong physical justifications for this behavior. 
The principal justification, as in many physical prob- 
lems, seems to lie in experiment. 


LIFT REPRESENTATION FOR RECTANGULAR WING-Bopy 
CONFIGURATIONS 


The lift representation for the family of rectangular 
wing-body configurations with no afterbody may be 
carried out in a similar fashion to that for the delta 
wing-body configurations. The two parameters for 
this case are } and a; The parameter a, is propor- 
tional to the effective aspect ratio? of the ‘‘exposed 
wing,” 


BA.R. 6 fl- ‘) 
2 c/a ( 
Again, the lifts for the limiting cases of the param- 
eter a, are easily obtained, since (A) a, = 0 implies 
either M = 1 (w # 7/2) or w = O (this is the case 
given by slender body theory), and (B) a, = ~ implies 


“strip theory’). Then the lifts for a, = O are pre. 
cisely those given for the delta wing-body configura- 
i.e., given by Eqs. (3a), (4a), and 


© obtained by a simple inte- 


tions for m = 0 
(5a). The lift for a, = 
gration is 


K = K™ = (1 +A) (25a) 


and 


K® = @ (25b) 


The total wing and body limiting lifts for a, = 0 and 
a, = © are shown in Figs. 4a, 4b, and 4c, respectively. 
The limiting body lifts are identical to those for the 
delta wing-body configuration, the most evident de- 
parture being in Fig. 4b for wing lift; the difference be- 
tween the limiting wing lifts is uniformly smaller than 
for the delta case, 


[KO (A, ©) — KO, O)] < 0.05K (A, 0) (26) 


Hence, just as in the delta case, it is of first importance 
to obtain approximate representations of the body lift 
for intermediate values of a,. Beforehand, it is noted 
that there are distinct intermediate values of a, charac- 
terized for example by the number of tip leading-edge 
reflections off the body; for decreasing a,, the first 
of these is a, = 1 (the tip Mach element intersects the 
trailing-edge-body junction), there being an_ infinite 
number in a decreasing sequence—i.e., a, = 1, | 2, 
1/3,...,. In contrast, there is only one distinct inter- 
mediate a,(=m) for the delta case—namely, a, = | 
(sonic leading edge). 
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(1) Approximate Body Lift for0<a< = 


ING-BODY LIFT 


For a, > 1, assuming as in the first approximation for body lift in the delta case that the pressure on the body 


is constant (equal to the asymptotic value far downst 


(14f) are 


‘ ar 
K(® = 
( l 
l oman 
la, 
and 
rn ] ( r 
Ki® -_ - 
( | ) 1 — 
l — 
fa, 
where 
Cre, = (4 B){1 — (1/4a,)] (27c) 


Fig. 4c shows that the shapes of the curves defined 
by Eqs. (27a) and (27b) for a, = 1 and 2 are similar to 
the corresponding curves for the delta case (cf. Fig. 


2c). The body lift slopes at\ = 0 and A = | are 
O0K‘)(0, az) 1 ; 
v= (fora, > 1) (28a) 
on [1 — (1/4a,)] 
and 
OK ‘*)(1, a,) | ’ 
=— (for a, > 1) (28b) 
OA 2a,fa, — (1/4)] 
Foregoing the second approximation for a, => 1 be- 


cause the result is clear (cf. Fig. 2c), an approximation 
for body lift for 1/2 < a, < 1 is given assuming con- 
stant pressure on the body upstream of the boundary 
defined by the region of influence on the body of the 
tip leading edge and zero pressure downstream of this 
Fig. 5 shows the shape of the projection 
the 


boundary. 
on the x-z plane of this boundary for0 < A < 1; 
boundary projection is given by 


for 


and 


ar ; l 
[i _ ( )| [fora < | 
r 1—AzXz 1 + (x/2)a, 


ream), the body lifts corresponding to Eqs. (14d) and 


| N [ _ fl—-a : l m 
( ) |! — ( ) sin ( )| ford Pg | (27a) 
) 1—A~7A 1—A ar 1 + (7/2)a,; 


(27b) 


x I 


<= 14¥(1-2)h 
a =| + ( )| 
1—d2z Zz V1 — 2? 
tor ( ‘) < < | (29a) 
A Ba A 


and the helical extension of Eq. (29a) is given by 


x z z /) — 2 
= cos ( _ 2) (for > ud. A ) (30a) 
a Ba Ba Xr 


where 
oo = = = (* Sukh.4 — cos! s) (30b) 
Ba rd 
Thus, in the limit 4 = 1, the boundary is a helix; in 


the limit 4 = 0, the boundary projection on the x-y 
plane is the straight-line extension of the Mach element 
from the tip leading edge—i.e., as the span increases, 
the body, viewed from the tip leading edge, appears 
flat. Then the body lift is given by Eqs. (27a) and 
(27b) minus a correction 


K@e® = : (>) (29b) 
{1 — (1/4a,)] \Sw 


where 
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( ~ | 1+n[1 m2 (- “ai 
r? i iy TE = 3\ A + 


2d 


7 oo} — a V1 —}? (30c) 


V1 — r2/d < (1 — X/aad) < [(V1 — d2/A) + sin A] 


The range of validity is limited by the condition where 


the boundary defined by Eqs. (29a) and (30a) inter- 


sects the body base at the plane of symmetry (y-z 
For 1/2 < a, < 1, the body lift is plotted for 
Further approximations 


plane). 
a, = 1/2 and 3/4 in Fig. 4c. 
along these lines for 0 < a, < 1/2 do not appear fea- 
sible. Hence, the representation of wing lift for 0 < 
a, < © is next considered. 


(2) Approximate Wing Lift for 0 < a, < « 


Meaningful approximations for the wing lift are more 
difficult to obtain for the rectangular case than for the 
the difference 


[KA (rX, ©) — Kw (X, 0)] 


delta case; 


is uniformly about one-third of that for the delta case. 
For a, > 1, the wing lift can be obtained by graphical 
integration of the pressure distribution (unpublished). 
The results of the calculations for a, = 1, \ > [1+ 
(1 + 2a,)], are shown in Fig. 4b. By an analysis 
similar to that for the delta case resulting in Eq. 
(22b) and the paragraph that follows it, the behavior of 
wing lift at \ = 0 and A = | may be obtained; the 
results are analogous to the delta case. Thus it is now 
possible to indicate the approximate behavior of the 
total lift for intermediate values of a,. Intermediate 
curves for a, = 1 and a, = 1/2 are shown in Fig. 4a. 
For a, = 1/2, the wing lift was assumed to be well ap- 
proximated by the wing lift fora, = 1. The outstand- 
ing feature of the total lift representation for the 
rectangular case is the marked similarity to that for the 
delta case (cf. Figs. 2a and 4a): (A) for a, = O, the 
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FIGURE 5. COMPARISON OF REGIONS OF INFLUENCE OF TIP 
LEADING EDGE ON BODY FOR RECTANGULAR WING-BODY 


lifts are identical; (B) for a, = ©, the difference jy 
total lift between the two cases is uniformly less than 
10 per cent for 0 < A < 1, being zero at the end points 
(A = Oand A = 1); and (C) for intermediate values of 
a,, the shape of the lift curves are similar for the two 
cases and the end points are identical. It is emphasized 
that the similarity lies in the comparison of the /if 
ratios. 

REPRESENTATION FOR CLIPPED-DELTA WING- 
Bopy CONFIGURATIONS 


LIFT 


The lift representation for the three-parameter fam- 
ily of clipped-delta wing-body configurations, which in- 
cludes the limiting two-parameter families of delta 
wing-body (a,/m = 1) and rectangular wing-body 
0) configurations, is presented to indicate the 
forms with 


(a,/m = 
possible to other plan 
straight edges. Again, the details of the analysis are 
not carried out so completely as for the delta case, but 
The three parameters for this 


=), a, and a,/m, where 


generalization 


the extension is clear. 
case are taken to be a/s 
0 < a,/m < 1 (cl. Fig. Ic). 


Then the lifts for a, = 0 are preciselyt those given 
for the delta case for m = 0 [ef. Eqs. (3a), (4a), and 
(5a)]. Fora, = ©, the lift obtained by simple inte- 


grations is 


K = KW =1+4 


: r+ %( d? )+%( A yi NIG 
-_ 2 — n « } 
[2 — (a,/m)]L m\1—\AX m\1—A 


and K‘#) =(Q, which reduce to the delta case for a,/m = 
l and to the rectangular case for a,/m = 0. The 
limiting lift curves are plotted in Fig. 6 for a,/m = 


V3/3. 
(1) Approximate Body Lift for 0 < a,< 


For a, > 1 (also implies m > 1)—i.e., the condition 
that the Mach element from the leading-edge-body 
junction intersects the trailing edge—the first approx! 
mation for body lift, assuming that the pressure 1s col 
stant on the body [ef. Eq. (13)], is obtained. Ther, 

7 A recent discussion with Lagerstrom indicates that a proof 


for more general plan forms including camber exists for this 


equivalence.!° 
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after the procedure following Eq. (14a), the wing area the body area [cf. Eq. (14b)] is 
icf. Eq. (14c)] is 
~ \io —— | 2 29. . 1—A r 1—A . - 
Sw/a? = (B/a,)[2 — (a,/m) ][(1 d)/A] (32a) S;/a? = 28 ee. ef (32b) 
air ar 
For 
A> 1/[1 + (7/2)a,] and 


4/4 , ag a, l nN ar : 1—A ; 
Cie‘ —_ 3 - sec 1 = l _ sin ae (32c) 
BLr Va? — (a,/m)? a,/m [2 — (a,/m)] \l —A 1—A ar 


For \ < 1/[1 + (2/2)a,], the body area [cf. Eq. (14e)] is 
S,/a? = 2B{[(1 — A)/aA] — 1} (32d) 


: 4|4 ay a; ] A air 
Cre? = ; sec! "eo (32e) 
8B Lx Va? — (a,/m)? a,/m/ | [2 — (a,/m)] \l — A l—i 


where the basic lift of the “‘exposed wing’’ alone" is 


3 4 l \ " ay / = 
Cra = 1 — [1 — (a,/m)}? > (for a; > 1) (32f) 
( 


Ba, m[2 — (a,/m)] } 1, — (a; m)\ 


and 


For m > 1 (supersonic leading edge) and 1/2 < [1 — (a,/m)] < a, < 1 (for a,/m < 1/2) or 1/2 < a, < 1 (for 
a,/m > 1/2)—i.e., the condition that neither the Mach element from the tip leading edge nor the side-edge reflec- 
tion of the Mach element from the leading-edge-body junction hits the body—the body lifts C,,“ are given by 
Eqs. (32c) and (32e), but the basic lift is 
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FIGURE 6. CLIPPED-DELTA WING- BODY LIFT WITH NO AFTERBODY, a, /m= 
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0 B a,/m[2 — (a,/m)] \ a Tt Vai + (a,/m) 
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cos , _ cos”! > (32g) 
1 + (a,/m) Tt Va. — (a,/m) 1 — (a,/m) f 


form > landa,<1. The body lifts from Eqs. (32c), 
(32e), (32f), and (32g) are plotted in Fig. 6c for a, = 
2,1, and \/3/3. A second approximation for body lift 
may be made similar to the delta case [cf. Eqs. (17b) 
and (17c)], and an approximation for body lift for 
m > land 1/2 <a,< [1 — (a,/m)] (for 0 < a,/m < 
1/2) may be made similar to the rectangular case [cf. 
Eqs. (29b), (29c), and (30c)]; but, since nothing essen- 
tially new is contributed, these approximations are not 
presented. 

Also, for 0 < m < 1 (subsonic leading edge), ap- 
proximate expressions for body lift similar to those 
given for the delta case [cf. Eqs. (19.2), (19b), and (20) | 
can be obtained but are not presented here. 


(2) Approximate Wing Lift for 0 < a, < « 


For a, = 1, the wing lift can be obtained by graphical 
integration of the pressure distribution; but, since the 
calculations are more tedious than for either the delta 
or rectangular cases, a faired curve for a, = | is shown 
in Fig. 6b. The behavior of the wing lift and slope at 
\ = 0 and A = 1 can be investigated as in the delta 
case, but the results are clear and the details are not 
presented here. 

The representation of the total lift is now possible 
for the clipped-delta case over the entire range of the 
parameters. Curves for two intermediate 
1 and 1/3/3 are shown in Fig. 6a, where the 
1 is taken for both values of a,. 


values of 
a, = 
faired wing lift for a, = 

The extension of the analysis to other cylindrical 
body shapes (e.g., elliptical bodies) and to other wing 
plan forms with straight edges (e.g., swept trailing 
edge and/or inclined side edge) is now clear. Also, 
the representation of drag, moment, and center of 
pressure may be initiated along lines similar to those 
presented in this paper for lift. 


EXPERIMENTAL VERIFICATION 


Some experimental results that verify the present 
analysis in part are available, but the data do not cover 
a sufficiently large range of the parameters in order to 
determine completely the behavior of the lift ratio. 
Suggestions for future experimental work include the 
(A) Reliable measurements of the basic 
A fairly extensive series of 


following : 
wing lift curve slope C;,,- 


tests on delta wings have been carried out.'? (B) j 
systematic series of tests on delta wing-body configura. 
tions without afterbody with particular attention to 
the validity of the parameter m = B tan wo. (C) Inde. 
pendent measurement of the body (or wing) lift, and 
(D) measurement of afterbody lift. Also, since the 
theoretical calculations are difficult, pressure measure. 
ments on the body adjacent to the wing are of particular 
interest. 
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The Interaction of Shock Waves with 
Boundary Layers on a Flat Surface’ 


F. W. BARRY,? A. H. SHAPIRO, ano E. P. NEUMANN** 
Massachusetts Institute of Technology 


SUMMARY 


The development of supersonic compressors, supersonic 
diffusers, and high-speed aircraft points to the increasing im- 
portance of the interaction between shock waves and boundary 
layers. 

The experimental work reported here is intended to (1) pro- 
yide a better understanding of the nature of the shock-boundary- 
layer interaction, (2) serve us a guide and stimulus to theoretical 
work, and (3) develop an empirical method for predicting the 
effects of the interaction 

Experiments were performed on the reflection of an oblique 
shock from a boundary layer on a flat surface at a Mach Number 
of 2.05. The effects of shock strength and boundary-layer régime 
laminar, transition, or turbulent) on the interaction phenomena 
were explored. 


The results are in the form of schlieren photographs, constant- 
density contours found from interferometer photographs, and 
static pressure distributions at the plate surface. 


NOMENCLATURE 


Fig. 1 shows the meaning of some of the symbols and the 
regions to which the numerical subscripts refer. 


lL; = see Fig. 1 

M = Mach Number 

p = pressure 

R; = Reynolds Number = Vix1pi/m1 
R, = Reynolds Number = Vx:p1/m1 
J 


. = free-stream velocity 


x = horizontal distance from leading edge of plate, in. 
(see Fig. 1) 

6: = boundary-layer displacement thickness = 3.83x:/V R; 

€ = angle between stream line and plate surface, negative if 
flow is toward plate, deg. 

« = angular deflection of apparent edge of boundary layer 
(see Fig. 1), as measured from photographs, deg. 

@ = shock generator angle (positive as shown in Fig. 1), 
deg., used as a measure of shock strength 

uw = free-stream viscosity 

p = density 


The flow is from left to right in all figures and photographs. 
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INTRODUCTION 


General Considerations 


| aesce sap gn WORK WITH subsonic flows has 
justified the assumption that one may usually 
analyze a boundary layer independently of the main 
flow. Experimental work with transonic and super- 
sonic flows, on the other hand, has shown that for such 


flows one may not analyze a boundary layer inde- 
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Theoretical relations for ‘regular reflection’’ from a 
plane wall at a Mach Number of 2.05. 


Fic. 3. 


pendently of the main flow, except for those few cases 
where shock waves do not meet a boundary layer. 
That is, shock waves have been found to interact with 
boundary layers to cause significant departures from 
the theoretical flow, both in the boundary layer and in 
the main flow outside of it. 

The simple theoretical solutions for the reflection of 
a shock in a parallel, uniform, supersonic flow are too 
well known to need more than brief mention. 

A weak oblique shock reflects from a plane wall as an 
oblique shock that turns the flow parallel to the. wall 
(Fig. 2). The overall pressure rise is about twice that 
of the incident shock. Fig. 3 shows some theoretical 
relations for this type of reflection at a Mach Number 
of 2.05. 

An oblique shock reflects from a surface of constant 
pressure as an expansion wedge to preserve uniform 
pressure at the surface (Fig. 4). The overall deflection 
of the flow is about twice that across the shock. 

No account is taken of the boundary layer in the 
simple cases mentioned above. However, the subsonic 
layer next to the wall has a marked effect on the flow 
because pressure disturbances, which cannot propagate 
upstream in a supersonic flow, do propagate upstream 
i flow. Furthermore, the enormous 


in a_ subsonic 


1951 


pressure gradient that occurs in the incident shock jp 
the supersonic part of the flow cannot occur in a syb. 
sonic flow. Therefore, near the wall, the region jn 
which the pressure rises will extend both upstream and 
downstream of the point of incidence. 
the subsonic layer will thicken and deflect the stream 
lines in the supersonic region away from the wall. This 
deflection will give rise to “‘reflected’’ compression 
waves that originate upstream of the point of incidence, 
as shown in Fig. 5. These waves cross, and thereby 
modify, the incident shock. 


As a result, 


The incident shock will become more nearly normal 
as it traverses the supersonic portion of the boundary 
layer and will produce a series of reflected waves, !:? 
If the boundary layer thickens considerably or sepa- 
rates, there will be a thick subsonic layer next to the 
wall. Hence, the incident shock must reflect sub- 
stantially as an expansion wave (Figs. 4 and 5), since 
the subsonic layer cannot support a large longitudinal 
pressure gradient. The flow following the expansion 
will, in general, be inclined toward the wall. As the 
boundary-layer flow turns parallel to the wall, a series 
of compression wavelets are propagated into the main 
stream (Fig. 5). 

The system of reflected expansion and compression 
waves coalesces into a single oblique shock at some 
distance from the point of incidence. If the flow after 
this reflected shock is parallel to the wall, the reflected 
shock will be identical to the theoretical reflected shock 
of a “regular reflection.’’ Therefore, the overall 
pressure rise should be nearly the same as the theoretical 
rise (Fig. 3), even near the wall, because the differences 
in entropy increase for different parts of the shock sys- 
tem are not large. 

As a result of all these effects, it is evident that the 
shock modifies considerably the flow in the boundary 
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layer and that the changes in boundary-layer flow not 
only alter the incident shock but also give rise to a 
reflected wave system far more complex than the 
“regular reflection’’ from a plane wall illustrated in 
Fig. 2. The essence of the process, therefore, is a 
mutual interaction between the shock and the boundary 
layer. 


Previous Work 


The quantitative results of the present report are 
an extension of the qualitative investigations by the 
same authors! and contain the principal results of the 
investigation that is reported in detail in a thesis by 
Barry.’ 


Fage and Sargent‘ studied the interaction of normal 
and oblique shocks with a turbulent boundary layer 
on a flat wall at Mach Numbers up to 1.47. They 
present impact and static pressure surveys in the 
stream, wall-surface-pressure traverses and wall-friction 
traverses below the shock waves, and schlieren and 
Shadow photographs. They found that the shock 
bifurcates near the boundary layer when its strength 
exceeds a definite value, that the surface drag coefficient 
falls rapidly between the two bifurcated limbs, and that 
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Sketch of the apparatus. 


the shock-wave angles are influenced by interaction 
with the boundary layer. 

Liepmann, Ashkenas, and Cole,® * in an investiga- 
tion at a Mach Number of 1.25, found that a weak 
shock reflects from a laminar boundary layer on a flat 
surface as a combination of shocks and expansion 
waves and from a turbulent layer as an oblique shock. 

A few analytical studies have been published. Tsien 
and Finston’ applied linearized potential flow theory 
to a model that consists of a uniform subsonic stream 
bounded on one side by a wall and on the other side by 
a uniform semi-infinite supersonic stream. They found 
that the distance that pressure disturbances (caused 
by a weak wave in the supersonic stream) are trans- 
mitted upstream is proportional to the thickness of the 
subsonic layer and is a function of the Mach Numbers 
of the subsonic and supersonic streams. Their results 
can hardly be compared with the present experimental 
results because of the simplicity of their model. 

Lees* predicted the effects of Reynolds Number, 
shock strength, and Mach Number on the surface 
pressure distribution and on the location of the point 
of separation of a laminar boundary layer. He cal- 


* Further experimental data will presumably be given in 


reference 6. 
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Fic.7. Schlieren photographs, arranged in columns of constant Reynolds Number and in rows of constant shock strength; knife-edge 
parallel to plate; exposure: 1 microsec.; M = 2.05. 
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Fic. 8a (top). Illustrates effect of shock strength on the thickness of the boundary layer for R; = 600,000, traced from schlieren 


photographs. Fic. 8b (bottom). Illustrates effect of Reynolds Number on the thickness of the boundary layer for @ = 6°, traced 
from schlieren photographs. 
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SHOCK-WAVE AND 
culated the growth of the boundary layer up to the 
separation point by means of a modified Pohlhausen 
method and assumed a somewhat arbitrary pressure 
distribution to calculate the remainder of the flow. 
The experimental results of the present investigation 
are only in partial agreement with the theoretical 


prediction of Lees. 


EXPERIMENTAL PROCEDURE 


The investigation was conducted in the 8- by 8-in. 
tunnel of the M.I.T. Gas Turbine Laboratory at a 
Mach Number of 2.05. Fig. 1 is a cross section 
through the tunnel showing the shock generator and 
leading part of the flat plate, and Fig. 6 is a schematic 
view of the apparatus in the tunnel. 

Experiments were performed at four Revnolds 
Numbers—namely, about 110,000, 250,000, 600,000, 
and 1,200,000. For the highest Reynolds Number 
the boundary layer seems to be in the transition region. 
In one run, discussed separately, a wire was glued to 
the plate near the leading edge to trip the boundary 
layer. Shock deflection angles up to 7° were employed. 
Schlieren photographs were taken for every combina- 
tion of variables, and interferometer photographs were 
taken at the two intermediate Reynolds Numbers. 
Pressures along the plate surface near the mid-plane of 
the tunnel were measured for the three lowest Reynolds 
Numbers and for the tripped turbulent boundary 
layer. 

The boundary layer was formed on a flat plate that 
spanned the tunnel at zero angle of attack near the 
center of the tunnel. The foremost 4!/. in. of the 
plate were visible through the optical flats set into the 
tunnel sidewalls. 

The shock-generator airfoil could be adjusted through 
a range of positive and negative angles of attack while 
the tunnel was running. Most of the data were taken 
with the generator arms shown in Figs. 1 and 6 
where x; was about 1.9 in. By means of a second 
set of arms it was possible to make x, about 3.8 in., 
thus permitting runs at Reynolds Numbers up to 
1,200,000 and also a pair of runs with the same Rey- 
nolds Number (600,000) but different values of x. 
With each set of arms the Reynolds Number was 
altered by changing the tunnel pressure. 

Because it was impossible to keep the leading edge 
of the shock generator extremely sharp and because 
of boundary-layer formation on its surface, the gener- 
ator produced a compression wave followed by a 
rarefaction wave even when it was set at zero or 
negative angles. Therefore, for angles less than 
15°, the nominal angular setting of the shock generator 
is an unreliable measure of shock strength. 

The schlieren photographs were taken with a 5-in. 
schlieren apparatus’ at an exposure time of less than 
1 microsec. The knife-edge was parallel to the plate 
for all the photographs, which means that the photo- 
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graphs show density gradients normal to the plate 
surface. Because of the large length of light path 
through the tunnel (8 in.), refraction at high tunnel 
pressures produced such optical illusions as an ap- 
parent bending of the plate surface and an apparent 
kinking of a wire outside the tunnel (see Fig. 7). In 
some of the photographs the incident shock appears 
to be composed of as many as six parallel shocks. 
This, of course, is not possible. The departures from 
two-dimensional flow due to an observable oscillation 
of the shock and the boundary layers on the sidewalls 
account for at least part of this effect. Streaks of oil 
on the tunnel glass walls may also be seen in some of 
the photographs. 

The interferograms are time exposures taken with a 
Mach-Zehnder interferometer."° 


RESULTS WITH LAMINAR BOUNDARY LAYERS 


Qualitative Results 

The qualitative results obtained with the large tunnel 
described in this report are in accord with the results 
obtained previously with a 1.3- by 1.3-in. tunnel’ and 
are summarized below. 

The schlieren photographs of Fig. 7 illustrate the 
remarks of this section. 

The reflected wave system is seen to comprise a 
series of compression waves caused by the thickening 
of the boundary layer upstream of the point of inci- 
dence, an expansion region originating at the point of 
incidence, and a downstream compression region. 
Thus, the observed flow pattern is qualitatively similar 
to that predicted analytically and shown in Fig. 5. 

Consider any column of photographs in Fig. 7. For 
a small shock strength—i.e., for photographs near the 
top—the components of the reflected wave system are 
spaced closely together. The boundary layer seems 
to thicken slightly near the point of incidence and to 
become turbulent downstream of the region of inter- 
action. 

As the shock strength is increased—i.e., for photo- 
graphs near the bottom of the column—the first re- 
flected shock and the point at which the boundary 
layer begins to thicken move upstream. The down- 
stream boundary layer becomes much thicker and 
highly turbulent. The zone of interaction therefore 
spreads both upstream and downstream as the shock 
strength increases. The reflected expansion wave 
and subsequent shock still occur. 

The effects of a change in Reynolds Number, on the 
other hand, may be seen by considering any particular 
row. A decrease in Reynolds Number spreads out 
the region of interaction along the plate but does not 
alter the reflected wave pattern appreciably. 

Fig. 8 shows the effects of Reynolds Number and 
shock strength on the boundary-layer thickness to scale. 
Quantitative data on the behavior of the boundary 
layer are discussed below. 
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Effect on Boundary Layer ae 
; , a ip ee —--— io 
Inspection of the schlieren photographs (Fig. 7) Ri 1,200,000 xj ~=3.8 
shows that the boundary layer begins to thicken at a v —-— Rj ~ 600,000 
rather well-defined point. The location of this point, a _ 
point. iS po v Ri = 600,000 
as measured on the schlieren photographs, is given ie 
in Figs. 9 and 10. 4 —— Rj = 250,000 Pi =1.9 
Fig. 9, for example, shows that L,/x, increases rapidly o ---- Rj = 110,000 
with shock strength at first and less rapidly as the 5 
shock strength becomes large. The ratio L,/x, in- | = a a 
creases with decreasing Reynolds Number, thus 4 SY | 
showing that the effect of the shock is felt farther / | 
upstream as the boundary-layer thickness increases. i: lap. Les 
e a ia : ee oe 
Increasing x;, without altering R;, seems to increase 3 l 
both L,/x; and L,/é, slightly. The curves for the €t a —y> 
qv 
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0.8 | ; | | | 1 
awe Ba ono on | Bal 
site = | - 
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0.2 + boundary layer upstream of the point of incidence as a function 
4 a Sale = ae of shock strength for various values of Reynolds Number 
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Fic. 9. Plot of the ratio /;/x; as a function of shock strength for Xj =3.8 
various values of Reynolds Number. v Rj ~ 600,000 
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Ps ’ Fic. 12. Ratio of the thickness of the boundary layer after 
v / and before the region of interaction as a function of shock strength 
/ 
ly Pd S for various values of Reynolds Number; measured from photo- 
30 yi , graphs. 
/ 
/ 
/%/ 
20 Fi highest Reynolds Number (1,200,000) are dissimilar 
Y/ 79 to the other curves, especially in Fig. 10, probably 
A Pd because the boundary layer has entered the transition 
10 7Y region before the region of interaction. 
7 =e ° ° ° 
Le In Fig. 11 it is shown that for weak shocks the 
° at angular deflection of the apparent edge of the boundary 
- , P . . ? , layer seems to be approximately equal to the deflection 
- DEGREES = Sian : 
| ‘ ; " ident shock - but for 
Fic. 10. Plot of the ratio L;/é; as a function of shock strength of the flow through the incident shock (8); 
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SHOCK-WAVE AND 
layer deflection angle remains substantially constant. 
The maximum value of ¢, increases with decreasing 
Reynolds Number and, for a Reynolds Number of 
600,000, increases aS x, increases. This interesting 
phenomenon of a maximum value of boundary-layer 
deflection is probably associated with the maximum 
pressure gradient that a boundary layer in supersonic 
flow can tolerate. 

Fig. 12 indicates the increase in boundary-layer 
thickness through the region of interaction as a function 
of shock strength. This figure should be regarded as 
only qualitative, because the measurements from the 
photographs are subject to errors due to refraction and 
to uncertainty regarding the position of the edge of the 
turbulent boundary layer after the interaction. 


Pressure on Plate Surface 


Some selected pressure distributions measured at the 
plate surface are shown in Fig. 13. The short hori- 
zontal lines on the right indicate the theoretical pressure 
rise for a ‘“‘regular reflection’ (Fig. 3). In every case 
the measured pressure rise is close to the theoretical 
one. 

The short vertical dashes on the curves indicate the 
location of the point where the boundary layer begins 
to thicken (x,) and of the point of incidence (x;), both 
as measured on the schlieren photographs. The pres- 
sure rise begins close to the point of thickening (x,), 
thus showing good correlation between the two methods 
of determining the extent of the influence of the shock 
on the boundary layer. 

For the lowest shock strength (@ = 1°), the pressure 
rises and then falls slightly to a uniform value. The 
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Fic. 13. Static pressure distribution on the plate surface for 
selected values of shock strength and Reynolds Number; the 
measured values of x; and x, are shown by vertical lines; the 
theoretical downstream pressures for “regular reflection” are 
shown on the right; p; is free-stream static pressure. 
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Fic. 14. Pressure on the plate surface at the point of inci 


dence as a function of shock strength for three Reynolds Num 
bers; ~; is free-stream static pressure. 


slight drop is caused by the expansion that follows the 
incident shock at this low shock strength. 

For the larger shock strengths, the pressure rises 
for a short distance and then tends to level off at a value 
that depends on the Reynolds Number. The re- 
mainder, and indeed the largest part, of the pressure 
rise occurs after the point of incidence and is pre- 
sumably associated with the downstream reflected 
compression waves. 

Fig. 14 shows that the measured pressure at the 
point of incidence is fairly constant for the larger shock 


strengths and increases with decreasing Reynolds 
Number. This result is in accord with the variation of 


boundary-layer deflection angle shown in Fig. 11. 

By plotting p;/p, against R, for the regions of Fig. 
14 where the curves are substantially flat, it is found 
that the measured pressure rise to the point of incidence 
for large shock strengths varies approximately as 
1/(R,)'*, in accord with Lees’ theory. 


Density Fields 


Figs. 15a, 16, 17b, and 18 show the contours of 
constant density for four combinations of R,; and @. 
The constant-density lines are labeled with the corre- 
sponding density ratio, referred to the free-stream 
density. A band shift to the left in the interferograms 
indicates a decrease in density. 

The overall density ratio (ps/p,) and the density ratio 
across the incident shock (p2/p,) agree to within 3 per 
cent of the theoretical values for “regular reflection” 
(Fig. 3), except for Fig. 17b, where the interferogram 
does not extend far enough downstream to include the 
region where the flow is parallel to the wall. 

Fig. 17a was taken at a low tunnel density to ensure 
that refraction would not be troublesome. This 
interferogram, like the schlieren photographs, shows 
that a region of low density gradient lies next to the 
surface. It is believed that the boundary layer has 
separated from the plate, although there is insufficient 
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Fic. 15a (top). 
Fic. 15b (bottom). 


experimental evidence to show whether backflow exists 
at the surface. 


A Simple Model of the Wave Pattern 


The simplified model of the wave pattern discussed 
in this section is of interest for two reasons. First, 
the flow field of the model is at least a first step toward 
a means for predicting the effects of the interaction 
and is also suggestive for further theoretical work. 





Lines of constant density ratio (referred to free-stream density) obtained from Fig. 15b; 6 = 3°; 
Interferogram; @ = 3°; Ry = 610,000; AZ = 2.05; 1-sec. exposure. 
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R; = 610,000 


Second, it was desired to obtain information about 
the flow which could not be obtained directly from at 
experiment. In particular, the flow directions in the 
reflected wave system cannot be measured directly by) 
any available technique but may be found from the 
model. The accuracy of the model cannot be deter: 
mined from experimental data on flow directions but 
must be determined from other data, such as densities. 


However, if the model is accurate with respect t0 
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SHOCK-WAVE AND 
densities, it is also accurate with respect to flow direc- 
tions because waves of only one family, except the 
incident shock, exist, and, therefore, the density is 
related to the flow direction. If the shocks are ex- 
tremely weak, the flow is essentially isentropic and 
this relationship is unique. 

The method described here is not intended to be 
correct in details, since the compression regions are 
replaced by single oblique shocks for simplicity; nor 
is the method intended to describe the flow in the 
boundary layer where the assumption of a relationship 
between flow properties is invalid. Furthermore, the 
method is restricted to cases within the range of the 
experiments with laminar boundary layers, since it 
requires the use of experimental data. 

The following procedure is used to calculate a wave 
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cedure. 


strength and Reynolds Number. The location of the 
point where the boundary layer begins to thicken, f¢, 
is obtained from Fig. 9. A single oblique shock, 

is assumed to originate at this point and to have a 
deflection angle equal to the measured deflection angle 
of the edge of the boundary layer, e, (Fig. 11). The 
properties of the flow in region 4 may be calculated by 
assuming that the overall deflection angle and pressure 
ratio of shocks a and d equal those of shocks ¢ and 6. 
The flow in region 5 is found by assuming that the 
incident shock d reflects as an expansion wave to make 
bs = ps. The properties of the downstream shock 
e and of the flow in region 6 may be found by making 
the deflection angle of this shock equal to —«, and thus 
the downstream flow will be parallel to the wall. The 
thickness of the boundary layer upstream of the inter- 
action (6,) may be found from R, and x,, and the thick- 
ness downstream (4s) may be estimated with the aid of 
Fig. 12. In the region of interaction the edge of the 
boundary layer is assumed to be formed of straight lines, 
and thus the geometry of the figure is completely de- 
termined through a combination of model assumptions 
and experimental data. 

Except for small differences due to irreversibilities 
in the shocks, the flow conditions behind the region of 
interaction, as calculated by this method, should be the 
same as those for a “regular reflection’ (Fig. 3). 

Table 1 lists some pertinent flow parameters of four 
typical wave-pattern together with corre- 
sponding experimental data. Fig. 19 is a sketch to 
scale of the model whose parameters are listed in the 
right-hand column (@ = 6°, R,; = 120,000) and for 
which the actual wave pattern is shown in Fig. 7, 
bottom left. 


models, 
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Fic. 21. Static pressure distribution on the plate surface with 
a tripped turbulent boundary layer for two shock strengths. 
Two curves for a laminar boundary layer taken from Fig. 13 are 
shown for comparison; the values of the measured x; are shown 
by vertical lines. The theoretical downstream pressures for 
“regular reflection’’ are shown on the right; ); is free-stream static 
pressure 


TABLE 1 
6 3° 6° 6° 6° 
—6§00,000-—— 630,000 —260,000— 120,000 
R; Model Exp. Model Exp. Model Exp. Model 
Pa/pi 1.10 1.06 1.10 1. 1.14 1.16 1.20 
63 2°22’ 2°22’ 3°28’ 4° 38’ 
ps/p1 1.22 1.20 1.39 1.38 1.43 1.41 1.46 
& —1° 5’ —3° 38’ —2° 40’ —2° 
M, 1.84 1.82 1.75 1.67 1.61 
ps/p1 L410 61.10 1.10 1.20 1.14 1.16 1.19 
¢ —3° 38’ —9° 36’ —8° 42’ —7° 26’ 
M; 1.93 1.90 1.96 1.88 1.80 
1.27 1.22 1.59 1.54 1.57 1.53 1.57 


P6/P1 


By and large, the agreement between the parameters 
calculated according to the model wave pattern and the 
experimental values suggests that the calculated flow 
directions of Table 1 are reliable. 


RESULTS WITH TRIPPED TURBULENT BOUNDARY LAYER 


The results presented in the preceding section were 
obtained with boundary layers that had formed in a 
normal manner and which extended into the transition 
régime (for the highest Reynolds Number). In order 
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to extend the investigation to include turbulent 
boundary layers, the boundary layer was tripped by q 
wire glued to the plate surface near the leading edge, 
The data for this turbulent boundary layer wer 
obtained at a tunnel density equal to that used pre. 
viously for the case R; = 600,000 and x; = 1.9 in. 
Fig. 20 is a typical schlieren photograph. The 
reflected wave consists of a single shock rather than 
the complicated pattern shown in Fig. 7 for laminar 
boundary layers. The boundary layer 
thicken appreciably in passing through the interaction 


does not 


region. 

The pressure on the plate surface is shown in Fig, 21, 
together with two curves for a laminar boundary layer 
for comparison. A single, abrupt rise in pressure js 
observed at the point of incidence rather than the two 
successive increments in observed — with 
laminar boundary layers. If the 
boundary layers are extrapolated to higher Reynolds 


pressure 
data for laminar 
Numbers, they approach the curve observed for a 
turbulent boundary layer. 

Thus, there is a significant difference between the 
interaction phenomena associated with laminar bound- 
ary layers and those associated with turbulent boundary 
layers. It is not clear whether this difference is due 
mainly to the fact that the thickness of the subsonic 
layér is probably less in a turbulent boundary layer 
than in a laminar boundary layer or whether it is due 
chiefly to the different shear stress distributions in 
laminar and turbulent boundary layers. 


CONCLUSIONS 


The experiments described provide insight into the 
nature of the interaction process. This understanding 
should be helpful to the designer of supersonic flow 
equipment. 

Since none of the available theories adequately 
explains the complex shock-boundary-layer interaction, 
further theoretical work is necessary. The results 
presented here may serve as a guide to such work. 
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On the One-Dimensional Theory of Flame 
Structure 


RAYMOND FRIEDMAN* ann EDWARD BURKE 


7) * 


Westinghouse Electric Corporation 


ABSTRACT 


The dependence of the temperature distribution in a premixed 
gas flame on burning velocity, flame thickness, heat-release 
pattern, and the physical properties of the mixture is calculated 
from the fundamental equations. By an analogous procedure, 
the local concentration of a chemical species created in the flame 
isobtained; both distributions are shown to be of the same form. 
The general equations are solved for a variety of heat-release 
patterns (or species formation patterns), and numerical values 
are presented over a wide range of parameters, so that the in- 
fluence of the pattern on the distributions may be seen. Solu- 
tions are obtained both under the assumption of constant trans- 
port coeflicients—i.e., thermal conductivity and diffusion coeffi- 
cient—and under the assumption that these coefficients vary 
with temperature. The results, presented in tabular and graphi- 
cal form, may be used directly to relate experimental information 
of reaction-zone thickness, burning velocity, species concentra- 
tions, and temperatures in the flame zone. The flame measure- 
ments of Wolfhard are discussed in relation to the equations de- 


veloped in this paper. 


NOMENCLATURE 


= position parameter (dimensionless) 


¥ 

6 = temperature parameter for preheat zone (dimen- 
sionless) 

nN thermal conductivity (cal. cm.~! sec.~! °K.~!) 

p = density (Gm. cm.~*) 

o = concentration parameter (dimensionless) 

T = temperature parameter (dimensionless) 

¢ = flame parameter (dimensionless) 

Q = number of collisions per molecule in reaction zone 
(dimensionless) 

C = weight fraction of a species (dimensionless) 

C,, C. = constants of integration 

D = coefficient of diffusion (em.? sec.~!) 

b = reaction zone width (cm.) 

£ = mean free path (cm.) 

M = local Mach Number (dimensionless) 

Q = heat release rate (cal. cm.” * sec.~!) 

: = absolute temperature (°K.) 

W = rate of species formation (Gm. cm. ~* sec. ~!) 

( = mean specific heat at constant pressure (cal. Gm.~! 
oK ~1) 

m = temperature exponent of transport coefficient (di- 
mensionless) 

S = gas velocity (cm. sec. !) 

0 = mean molecular velocity (cm. sec. ') 

x = distance (cm.) 

Subscripts 

a = average 

b = completely burned gas 

c = constant transport coefficients 
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D = concentration gradients and diffusion flow 

0 = plane in which reaction commences (except for Case 
IV) 

i = temperature gradients and heat flow 

u = unburned gas 


variable transport coefficients 


INTRODUCTION 


U’ TO THE PRESENT, the equations for the struc- 
ture of one-dimensional, steady-state, premixed 
gas flames have been derived in four ways: 

(A) Rigorous simultaneous differential equations 
for heat flow and for diffusion flow of each component 
are formulated. ' 

(B) Heat flow is assumed to be the limiting process, 
and diffusion is ignored.?~* 

(C) Diffusion, particularly of free radicals, is as- 
sumed to be the limiting process, and heat flow is 
ignored.°—* 

(D) Both heat flow and diffusion are considered, but 
assumptions are made which restrict the generality of 
the treatment.*~'* 

No one of these approaches has yet been conclusively 
demonstrated to be applicable to flames of practical 
interest, although approach (A) offers the most promise 
of leading to an ultimate solution. For flames such 
as the hydrogen-bromine and hydrogen-oxygen sys- 
tems, where the chemical kinetics is moderately well 
understood, application of this rigorous approach may 
demonstrate that either heat flow or diffusion can be 
ignored under certain conditions. However, for hydro- 
carbon flames even a qualitative knowledge of the 
chemical mechanisms is lacking, whereas the theory re- 
quires quantitative values for the reaction rate con- 
stants, as well as for their temperature coefficients. 
Since an immediate application of the rigorous ap- 
proach to hydrocarbon flames cannot be expected, it 
is of interest to develop approaches (B) and (C) as far 
as possible. 

In this paper, approaches (B) and (C) are used to 
derive the relations between the flame variables in such 
a way as to be applicable to models based either on heat 
flow or on diffusion while investigating the influence of 
various patterns of heat release or of active species for 
mation. Furthermore, solutions are obtained both for 
the usual assumption of temperature-independent 
transport properties—e.g., thermal conductivities and 
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diffusion coefficients—and for a more realistic assump- 
tion wherein these coefficients vary with temperature. 
The results, presented in tabular and graphical form, 
may be used directly to relate experimental informa- 
tion on reaction-zone thickness, burning velocity, and 
temperatures and species concentrations in the flame 
zone. 
The following dimensionless parameters will be used: 
r =(T —-T,)/(% — T,) 
g@ =(C — C,)/(G — C,) 
~ = eo L 
Gr = S,LCp,/Xo 
Yp = SyLpu/ Dopo 
h = QL/[supuc (Ty, — Ty)] 
i = WL/[supu (Co — C,)]) 


Here, the subscript wu refers to the unburned and 6 to 
the burned gas; 7’ is temperature; C is weight fraction 
of a particular species; x is distance (gas flow being in 
positive x-direction); JL is reaction-zone thickness; 
Sy is burning velocity; A is thermal conductivity; c 
is average specific heat at constant pressure; p is den- 
sity; and D is diffusion coefficient of a particular spe- 
The reaction is considered to take 
place between x = 0 and x = L. If Q denotes the 
local heat-release rate per unit volume and JV the local 
rate of species formation per unit volume, then the 
variation of Q with x between x = 0 and x = L will 
be referred to as ‘“‘the heat-release pattern,’’ appearing 
in dimensionless form as h, while the variation of W 
with x in the same region will be referred to as ‘‘the 
species-formation pattern’ and appears in dimension- 
The subscript 0 in the definitions refers 


cies in the mixture. 


less form as j. 
to the plane x = 0. 

The steady-state equations for conservation of 
energy and of mass may be formulated, and explicit 
solutions may be obtained by assuming an arbitrary 
dependence of Q or W upon x between 0 and L. When 
transposed to the dimensionless groups, such solutions 
will be of the form 


II 


0 
0 (2) 


F(¢gr, 7, ¥) 
G(¢p, 0, 7) 


¢, which contains the burning velocity, may be com- 
puted only when some value of 7 (or a) other than zero 
and unity is known. As an example, the r-value cor- 
responding to y = 0, written 7, is the temperature at 
which reaction is represented as beginning, 0) being an 
analogous critical concentration of free radicals or other 
species required for reaction. Substituting y = 0 into 
the pair of Eqs. (1) and (2) gives 
F(gr, to) = 0 (3) 
G(¢gp; oo) 0 (4) 
The Mallard-le Chatelier equation’ is an example 
of Eq. (3), while the Gaydon-Wolfhard equation® may 
be shown to be an instance of Eq. (4), when it is recog- 
nized that their assumptions are equivalent to o) = 


1/e. 


I 
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Equations of the forms represented in Eqs. (1) to 
(4) are obtained in this paper for four heat-release anq 
species-formation patterns. The effect of temperature 
dependence of \ and (Dp) is investigated for two of these 
cases. 


GENERAL EQUATIONS AND BOUNDARY CONDITIONS 


The steady-state, one-dimensional heat balance may 


d dT dT 
rN ) — Spc : + Q = 0 (5) 
dx 


dx dx 


be written 


Here, specific heat is assumed constant. The product 


sp is constant by the equation of continuity 
SP = SuPu = SoPo (6) 


The equation for conservation of a particular species 
may be written! 


d (v ~<) dC +w=0 a 
dx f dx ” dx 7 " 


This equation would be rigorous only if partial pres- 
sure were used instead of weight fraction C in the first 
term, but no appreciable error is introduced if either of 
the following conditions is fulfilled: (1) if the con- 
centration of the species is small, or (2) if the molecular 
weight of the species is not far different from the aver- 
age molecular weight of the mixture. 

The following pair of conditions on Q and W may be 
obtained by integrating Eqs. (5) and (7) between minus 
infinity and plus infinity and by noting that the first 
derivatives vanish at these limits: 


o 71 
/ QO dx = / Odx = sy,c(T,—T,) (8) 
=% J0 
> 71 
/ W dx = / W dx = syp,(C, — Cy) (9) 
=e J0 


These equations may be transformed to a dimension- 
less basis, yielding 


d (; *) dr i 0 (10) 
—_ Por ‘al 1 = 
dy \Xo dy vie dy _ 


d (52 “e) do 4 oi=0 (11) 
dy \Dopo dy - dy “= ; 


1 
fi dy (12) 
0 
*1 
I jdy=1 (13) 
0 


Since Eqs. (10) and (12) are strictly analogous to Eqs. 
(11) and (13), only the equations of heat flow will be 
carried through from this point. Eq. (10) may be 
integrated directly to give 


Ad 
Lerten fhay=G 
Xo dy 


II 
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ONE-DIMENSIONAL 


The boundary condition required to evaluate the con- 
stant of integration C is obtained from the heat balance 


across the plane x = 0, 


[A(dT/dx)],=0 = Subuc(To — Ty) (15) 
or, in dimensionless form, 
(dr dy), 0 = ¥TTo (16) 


The integration constant is then found from Eqs. (14) 
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and (16) to be given by 
C, = or( Sh dy),=0 (17) 


It will be assumed that the transport coefficients in the 
region where y is between 0 and | may be represented 
by the relations 


Upon substitution of Eq. (18) into Eq. (14), the latter may be integrated a second time, yielding 


(1 in ie . (Ab “~ > (Cy —— f h dy)dy _ ) 
T= — : 2¢ 
—  * UT IL 4+ y (a/de) — yx] trertr/aamay Tf 


The new constant of integration C, is obtained from the boundary condition that r = 1 at y = 1. 


d/o = 1 + y[(Ao/Ao) — 1] (18) 
Dp Dopo => l + v[(Dop» Dopo) —_ 1] (19) 
(20) 


A strictly anal- 


ogous equation may be written exhibiting o as a function of y, ¢p, j, Dopo, and Dypy. 
For the simpler case where the transport coefficient is taken as independent of temperature, Eq. (14) leads to 


T = exp (¢ry) | / (c — er fi iy) exp(—¢ry)dy + c.| (21) 


with an analogous equation for o. 


EQUATIONS FOR PREHEAT ZONE 


In the preheat zone (y < 0), # and j are defined equal 
to zero, so that the terms in Eqs. (10) and (11) con- 
taining these variables drop out, and solutions are 
therefore independent of the pattern of heat release or 
species formation. 

When ) is treated as constant, another term drops out 
of Eq. (10), and the resulting equation is readily inte- 
grated, with the boundary conditions rt = 7) at y = Oand 


dr/dy = r = Oat y = — ~, to yield 
dr/dy = grt (22) 
T = 7 exp (¢ry) (23) 
Similarly, 
o = oo exp (¢pY) (24) 


When \ is treated as variable, the first integration of 


Eq. (10) gives 
(A/Xo) (dr/dy) = grr (25) 
Similarly, 


(Dp/ Dopo) (da/dy) = gpa (26) 


In order to integrate Eq. (25), it is necessary to assume 
a relation between thermal conductivity and tempera- 
ture. Experimental data for air indicate that the fol- 


lowing relation would be satisfactory : 
A\/XAo = Dp Dopo = (T, To)’ f (27) 


In Eq. (27), T and Ty denote absolute temperatures. 
When Eq. (25) is integrated using the three-quarter 


power law, it is found that 


. = 3 [1-(F) | + 2m fe 
Gry = 3 A U9 ) an 
; | "+ 1\ /@“ — 1\)) 
, 1 4 . ‘ 
tan @’*+ 9 log (7 — ‘) = 4 I (28) 


where 


19," = 


r(T7, —_ i ee 
9 = 4 


a (=7/T,) 
r. 


5 = 


to(7, — Ty) — 
T +3 (=70/T.) 
u 
An analogous solution for Eq. (26) is obtained upon re- 
placement of r and ¢7 in Eq. (28) by o and gp. The 
latter solution is valid, however, only insofar as the ap- 


proximation 
T + + (T, seg T.,) | _ o - ise (1, ead T.) | 


= ae (29 
to + [Tu/(To-— Tu)] 0 + [Tu/(T — Ti) ] 


is satisfactory. This approximation is necessary be- 
cause Dp is, in a physical sense, a function of tempera- 
ture, whereas concentration appears as the dependent 
variable in Eq. (26). The magnitude of the error in- 
volved in Eq. (29) may be seen from the following ex- 
pression obtained from Eqs. (23) and (24), which is 
valid only for constant \ and Dp: 


t/t = (a/o0) exp [(¢p — ¢r)¥] (30) 


The approximation is seen to be good as long as either 
(¢p — ¢r) or y is small. 
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Fic. 1. Heat-release and species formation patterns. 


SOLUTIONS FOR VARIOUS HEAT-RELEASE AND SPECIES- 
FORMATION PATTERNS IN REACTION ZONE 


Any analysis of the heat-release pattern demands 
consideration of the exponential temperature depend- 
ence of chemical reaction rates; this suggests an ex- 
ponentially increasing heat-release pattern in the re- 
gion0 <x < L. However, interdiffusion modifies this 
distribution, making it less steep. Finally, the heat 
release would tend to decrease in the downstream por- 
tion of the flame, as the reactant concentrations be- 
come depleted—i.e., the function must possess a maxi- 
mum. 

The free radical formation pattern is more difficult to 
predict, since it depends on the specific chain mecha- 
nism that is effective. It is even possible that the free 
radical concentration goes through a maximum because 
the rate of recombination might exceed the rate of for- 
mation in the downstream part of the flame. 

The procedure in this paper will be to choose arbi- 
trarily certain forms for the heat-release and species- 
formation patterns and then to investigate the corre- 
sponding relations between the variables. The follow- 
ing cases were chosen and are sketched in Fig. 1: 

(1) Q (or W) constant for 0 < x < L and zero for 
x <Oandx > L. 

(II) Q (or W) directly proportional to x for 0 < 
x < Land zero forx < QOandx > L. 

(III) Q (or W) directly proportional to [1 — 
cos (24x/L)| for 0 < x < L and zero for x < 0 and 


o> LL. 
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(IV) Q (or W) directly proportional to exp(x/L) fo, 
—o <x < Oand zero forx > 0. 

The symbols Ic, IIc, IIIc, and IVe are used for those 
cases in which ) or Dp is assumed to be constant, and 
the symbols Iv and IIv are used for those cases in which 
they are assumed to vary with temperature. The 
proportionality constants introduced in each of these 
four cases are evaluated from Eqs. (8) and (9). [It 
must be noted that one of the boundary conditions for 
Case IV is different from the other cases, so that cer. 
tain minor modifications must be made in the pre- 
viously derived equations in order to treat this expo. 
nential case. 

Solutions for the various cases in the reaction zone are 
listed below. In every case, it is necessary only to 
replace rt and ¢7 by o and ¢p in order to transform an 
equation for temperature distribution into the corte- 
sponding one for species distribution. 


Case Ic.—The general solution of Eq. (21) is found 
to be 
r= {1 + gry — exp[—er(1 — y)]}/er (31 
for0 < y < 1, and the solution for y = Ois 
™m% = [1 — exp (—¢7)]/¢r (32) 


It is interesting to compare Eq. (32) with the Mallard- 
le Chatelier equation,’ which is, when written in terms 


of 7) and ¢7, 


l 1 — [1/(¢r + 1)] ~ 
T = = (33) 
¢r t+ | ¢r 
Case IIc.—The general solution for0 < y < lis 
9 22,2 
2 ery 
r2— dy + ery + — = 
¢r° \ 2 


© ce 
(gr + 1) exp [—¢7(1 — Wy (34) 


and the solution for y = O is 
(35 


tT) = (2/g7?) [1 — (gr + 1) exp (—¢7r)] 


Case IIIc.—The general solution for 0 < y < lis 


l Y7 ] 
T= +7+| : <n exp[—¢r(1 — )] - 
YT gr’ + 41° Yr 
(gr?/2mr) sin (2ry) + gr cos (27y) ,. 


- E (50) 
gr? + 4x? 


and the solution for y = Ois 


| = | exp ( ) 1 (37) 
"m= 2 “Tas XP \— Gr) — ' 
: gr + 4n° Yr 


Case I Vc.—Since in this case the reaction zone is in- 
finitely wide, L is defined as the distance in which the 
function Q (or W) decreases to 1/e of its maximum 
value. The solution for —» < y < Ojis 


: | 
= exp iy) = exp aie (38) 
or — l Dr a 
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ONE-DIMENSIONAL 


and 
_ gr exp (gr — 1) — 1 
(gr — 1) exp (¢7) 


(39) 


In the cases where \ varies, a new variable 
The solution of Eq. (20) for 0 < 


Case Iv. 
\,»/Mo is involved. 
y < lis 


No [1 + (ro/Ao)¥ — ¥ | TO4/%)) — 1) 
1+ ¢ry — 


Xo Ap/Ao 
l + Gr~— (Ap do) 


(40) 


and the solution for y = Ois 


1 — (Ap/Do)E O0/M0) — Fer y/%y) =H 
_— (41) 
I ss Yr~— (Ay do) 


d, 7 / ((Ab/dv) — 1] 
1+ ¥ _~? | : 
aii d/o ‘iors 2| | 
| e x i) 


and, when y = 0, 


men fei ™ 
b 


The comments concerning Case Iv also apply here. 


DISCUSSION OF RESULTS 


Numerical values have been computed over a wide 
range of the variables for each of the above equations 
and are given in Table 1. It may be noted that the 
equations generally become indeterminate at g, = 0. 
However, in each case, by a suitable expansion in series, 
it may be shown that the limit of r¢ as gr approaches 
0isl. The details will not be presented. 

Fig. 2 shows the solution for Case Ic; all of the other 
cases possess solutions of this same general form. The 
shape of the curve for y = 0 reveals the inverse nature 


of the relation between 7) and ¢7. As previously noted, 
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The ratio \,/Ao may be expressed in terms of 79 as fol- 


lows: 
a ie (re eee” 
do Dopo To af UT, (T, = T.) ] + tT! 


Plotting Eqs. (41) and (42) over the range of the 
variables makes it possible to eliminate \)/A» between 
the two plots by preparing a third plot relating ¢7, 
7m, and T,/(T, T,,). It may be noted that, for 
actual flames, 7,,/(7;, — 7.) is almost always between 
0.075 and 0.5. 

The corresponding relations between gp, 7, and Dypp + 
Dopo are seen to be rigorous only when 7) = oo because 
of the dependence of Dp on 7 instead of on oo in Eq. 
(42); fortunately, only small errots are introduced 
when oy differs moderately from 7». 


(42) 


Case IIv.—The solution for 0 < y < lis 
) No _ (*) ; 
, ”™ ™ Fag : 
Xp Jams ] as , ro a ee ae 
No i > No ig 
No No : 
2\i ty — vi}er ity, 9 ?7 
0 0 
l + -_ (43 
Np No ) 
r =— i — o 25 —2— oy 
ro» \? 
© Bs 
—_ — 
a ? r» a 2 a & 
7 ro 2 7 
l 2 gr Yr 
a bee = y 
(*-1) on tm oe oe ” 
Xo ro 7 Xo ~ ¢T 


the curves in Fig. 2 also express the relation between 


o and gp. 


The solution is again shown in Fig. 3, where curves of 
constant ¢ are plotted instead of constant y; this plot 
shows the gradient through the flame. Figs. 4 through 
8 are analogous to Fig. 3, except that they represent 
Cases IIc, IIIc, IVe, Iv, and IIv, respectively. 


A comparison of Cases Ic, Ic, and IIIc is given in 
Fig. 9 showing gr vs. 7) at y = 0. Fora given value of 
¢r, it may be seen that Case Ic yields the highest, and 
Case IIc the lowest, value of 7). The solution of the 
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Plot of r vs. y for various values of gp. 
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TABLE 1 


Values of r for the Various Cases 





} Tu 
Case Tb — Tu |y ~~} 0 0.4 0.8 1.3 1.6 
} 9 1.000 0.371 0.140 0.053 0.020 
| | = 1.000 0.555 0.312 0.176 0.100 
le 25 1.000 0.678 0.465 0.320 0.223 
0 1.000 0.828 0.695 0.584 0.497 
0.25 1.000 0.904 0.818 0.744 0.688 
} 0.5 | 1.000 0.956 0.915 0.879 0.8 
—2 1.000 0.342 0.120 0.042 0.015 
a 1.000 0.511 0 267 0.140 0.075 
IIc an —0.5 1.000 0.625 0.398 0.256 0.166 
0 | 1.000 0.763 0.595 0.466 0.371 
0.25 | 1.000 0.844 0.726 0.632 0.551 
| | 0.5 1.000 0.924 0.857 0.796 0.744 
| <i 1.000 0.369 0.137 0.051 0.019 
=i 1.000 0.551 0.304 0.169 0.095 
IIIc : —0.5 | 1.000 0.674 0.454 0.307 0.211 
0 1.000 0.823 0.677 0.560 0.470 
0.25 1.000 0.904 0.820 0.751 0.679 
0.5 1.000 0.970 0.943 0.920 0.898 
—3 1.000 0.455 0.265 0.166 0.121 
1Ve —2 1.000 0.658 0.469 0.360 0.291 
| 1.000 0.881 0.781 0.701 0.643 
—0.5 1.000 0.967 0.936 0.904 0.876 
0.14 —2 1.000 0.280 0.021 
0.14 —1 1.000 0.518 0.239 0.093 0.031 
0.14 —0.5 1.000 0.666 0.438 0.291 0.188 
Iv 0.1 0 1.000 0.824 0.700 0.611 0.540 
0.3 0 1.000 0.821 0.694 0.605 0.531 
0.5 Q 1.000 0.818 0.688 0.599 0.522 
0.14 0.25 1.000 0.908 0.832 0.776 0.730 
0.14 0.5 1.000 0.960 0.927 0.900 0.876 
0.14 2 1.000 0.267 0.028 ae 
0.14 1 1.000 0.492 0.218 0.081 0.03 
0.14 0.5 | 1.000 0.633 0.398 0.248 0.15 
IIv 0.1 0 1.000 0.790 0.935 0.521 0.435 
0.3 0 1.000 0.784 06.625 0.509 0.42 
0.5 0 1.000 0.778 0.615 0.497 0.41: 
0.14 0.25 1.000 0.870 0.779 0.708 0.644 
0.14 0.5 | 1.000 0.910 0.850 0.809 0.779 


Mallard-le Chatelier equation is shown to be reason- 
ably close to that of Case IIc for ¢7 less than 2. 

Figs. 10 and 11 show comparisons of the results ob- 
tained according to assumptions of constant and vari- 
able transport properties for t) as a function of ¢7; 
Case I is treated in Fig. 10 and Case II in Fig. 11. It 
will be observed that a family of curves appears for 
each of the variable cases corresponding to the pre- 
viously discussed 7j{7, — T,). The 
curves are seen to be insensitive to the value of 7, + 
T,-—T,). Fora given ¢, an increase in this parameter 
corresponds to a slight decrease in 7». 

Of the various flame parameters occurring in the 
treatment, fairly good experimental data are avail- 


parameter 
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able for s,, A, py, c, and 7. However, few values, and 
these of uncertain reliability, have been reported for L 
and 7) and none at all for concentrations of active 
species except in the region y > 1. 

Values for L and Ty are given by Wolfhard, et al.,® ' 
fora lean acetylene-oxygen flame at 4.3mm. Hg. The 08 
consistency of these measurements may be determined 
by comparison with the previously developed relations 
between yg, and 7». From Wolfhard’s results, 7) = 0.6 
800°C., 7, = 2,950°C., and 7,, = 30°C., so that 7 is 


found to be 0.26. To compute yz, values must be ob- 


? : es ; 4 
tained for Xo, py, ¢, L, and s,. The first three quanti- ° 
ies were computed by averaging published data for 
acetylene and oxygen. JL is taken to be 12 mm., the aa 


value given for the total thickness of the luminous 
zone. The value of s, was estimated as recommended 
by Wolfhard" by considering the flame to be flat and r) 
equal in diameter to the burner tube; since the volu- 


metric flow is given, s, may be computed then to be Fic. 
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260 cm. per sec. Utilizing these values, yg; is found 
to be 3.8. 

From Figs. 10 and 11, corresponding to Cases Iv 
and IIv, values of 7) corresponding to the above gr are 
found to be 0.31 and 0.21, respectively, with 7,,/(7, — 
T,,) = 0.10. The values of 7) = 0.31 for Case Iv and 
7 = 0.21 tor Case Ilv may be compared with the ex- 
perimental value of 7 = 0.26. Since Cases Iv and 
IIv probably represent two extreme forms for the heat 
release, the values of 7) and L given by Wolfhard, 
et al., seem to be mutually consistent. 

Values of 7> and reaction zone thickness have been 
determined from spark ignition energies by Lewis and 
von Elbe.” This topic has been discussed in a pre- 
vious note by the writers. '® 

When more precise measurements are made for a 
variety of flames, the relations developed in this paper 
will aid in the understanding of the mechanism of com- 


bustion. 
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Appendix 
MOLECULAR INTERPRETATION OF THE RESULTS 
From the kinetic theory of gases, 
\/Cpg =~ D ~v£ (45) 
where'p, is an average density, £ is mean free path, 


and v is mean molecular speed. From Eg. (45) and 


the relation s,pq = S,p,, it follows that 


Sq/t = ¢ (L£/L) (46 
But s,/0  M, so that 
M = o(L£/L) (47 


This interesting equation states that the local Mach 
Number of a flame is always far below unity, since 
the mean free path is small compared with the flame 
i.e., only a small fraction of the molecular 
This concept 


thickness 
collisions in a flame lead to reaction. 
may be developed further: The number of collisions 
undergone by a molecule traversing the reaction zone, 


to be denoted by Q, is given by the product of collision 


rate (v/£) and residence time (L/s,). Thus, 
2 = 7/52 (48) 
A combination of Eqs. (45), (47), and (48) yields 
U~ Vo/2 (49 


The term 1/Q2 is a measure of the probability that any 
particular collision of a reactant molecule with any other 
molecule will be the critical collision producing reaction. 
The term ¢ is seen from Figs. 10 and 11 to be near 
unity as long as 7) is near 0.5. Thus, the Mach Num- 
ber of the flame is of the same magnitude as the square 
root of the reaction probability. 

It must be noted that the residence time upon which 
Eq (49) is based is no longer accurately given by L/% 
when the molecule in question is consumed by reaction 
The following reasons are given for this: (1) The mole- 
cule will react before traveling the entire distance L, 
so that L/2 should be used for Case I and Case III 
flames and L/+/2 for Case II flames. (2) The reactant 
molecule will possess a diffusion velocity in the down- 
stream direction superimposed upon the gas velocity. 
The magnitude of this latter effect may be estimated. 
Such an estimate, based on Case I, shows that, for gp = 
1, the velocity of reactant molecules is about three 
times the average velocity. Thus, both of these effects 
act to make the originally calculated residence time 
and, therefore, 2 too large. More refined calculations 
would give further insight into the physical processes 


taking place. 
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The Laminar Boundary Layer on a 
Rotating Blade’ 


LAURENCE EUGENE FOGARTY? 


Cornell University 


SUMMARY 


Boundary-layer equations of motion for a blade rotating 
steadily in an incompressible fluid are developed. The equations 
show that, at stations several chord lengths from the axis of rota- 
tion of a thin blade, the chordwise velocity is independent of the 
spanwise velocity due to rotation; hence, the separation line is 
unaffected by rotation. The detailed calculation of the spanwise 
velocity distribution in the boundary layer is carried out for two 
Spanwise velocities in the boundary layer are found 
No large effects 


examples. 
to be small compared with chordwise velocities. 


of rotation are found. 


INTRODUCTION 


2 )UGH THE BOUNDARY-LAYER THEORY formulated 
by L. Prandtl in 1904 is of great importance in 
practical aerodynamics, the difficulties presented by the 
nonlinear boundary-layer equations have limited the 
classes of solutions to flows of special symmetry. Most 
of the cases that have been treated to date are two- 
dimensional in character; solutions of three-dimen- 
sional boundary-layer flows are rare. 


One of the practically important types of three- 
dimensional boundary-layer flows is the flow on rotating 
blades, such as turbine, helicopter, and propeller blades. 
In this case, the centrifugal and Coriolis forces due to 
rotation, combined with the pressure gradients and 
viscous forces, cause the flow to be three-dimensional 
even in the absence of effects of finite span. von Kar- 
man’s solution of the problem of an infinite rotating 
disc' is probably the only solution of a problem of this 


type. 


As part of a program of investigation of three-dimen- 
sional boundary-layer flows, the laminar flow over a 
rotating cylindrical blade is treated here. Although the 
present investigation deals only with that portion of a 
blade of infinite span which is several chord lengths 
from the axis of rotation, it is believed that the results 
are of considerable interest. The results should be 
directly applicable, to good approximation, over a large 
portion of the span of high-aspect-ratio blades. Fur- 
thermore, it appears that treatment of the problem of a 
finite blade, or of a blade in the presence of a wall, for 
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example, will make use of the results of the present 
investigation. 


(1) BouNpARY-LAYER EQUATIONS 


The time variation can be eliminated from the 
Navier-Stokes equations by fixing the axes with re- 
spect to the rotating blade. Let the Z axis be taken 
along the axis of rotation, the Y axis along the span 
of the blade, and the X axis so as to form a right-handed 
Cartesian system (see Fig. 1). The Navier-Stokes 
equations for steady motion referred to the rotating 
reference frame are” 
u'uy’ + v'uy’ + w'uz’ — 2v'R — XR? = 

—(1/p)px + v(uxx’ + uyy’ + Uzz’) (1) 


u'vy’ + v'vy’ + w'vz’ + 2u’R — YR? = 
—(1/p)py + v(vxx’ + vyy’ + 22’) (2) 


u’'wy’ + v’'wy’ + wwe’ = —(1/p)pz + 
v(Wx x’ + Wyy’ + Wzz') (3) 
uy’ + vy’ + wz’ = 0 (4) 


Here, u’, v’, and w’ are the velocity components in the 
directions of the rotating axes, X, Y, and Z, respec- 
tively; R is the angular velocity of the blade; p is the 
density; p is the pressure; and » is the kinematic vis- 
cosity. 

Following Goldstein,* one introduces boundary-layer 
coordinates (Fig. 1) in which x is a curvilinear coordinate 
parallel to the blade surface, z is normal to the blade 
surface, and the x-z plane is normal to the blade lead- 
ing edge, which is taken as the y axis. For those parts 
of the blade where the radius of curvature is large and 
nearly constant and where the angle between the z 
axis and the axis of rotation is small, the equations of 


motion are 
uu, + vu, + wu, — 2vR — xR? = 
—(1/p) pe + v( ther + thyy + ez) (5) 


uv, + vw, + we, + 2uR — yR? = 
—(1/p)py + v(Yrr + Vy + V2) (6) 


UW, + vw, + ww, = —(1/p) Pz + 
V(Wrr + Wyy + Wee) (7) 
Ue, + vy + Ww, = 0 (S) 


If the usual boundary-layer approximations for 
fluids of small viscosity are made in these equations, 


the boundary-layer equations of motion referred to the 
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Fic. 1. Sketch showing} boundary-layer coordinate system on 


rotating blade. 


rotating coordinate system result: 


uu, + vu, + wu, — 2vR — xR? = 


—(l/ pips Tr Vis (9) 
uv, + vy + we, + 2uR — yR? = 
—(1/p)py + vz, (10) 
pz: = 0 (11) 
Ur + vy, + w, = O (12) 


The boundary conditions are that there be no fluid 


velocity relative to the surface at the surface and that 


the flow approach the potential flow far from the sur- 
face. 
u=v=w= 0, z= 0 
u—> Uy , 
(13) 
vU—>?) sa © 


. a 
Ww-—, WwW, |} 


ICAL SCIENCES—APRIL, 1951 


where 1, 2, W; are the velocities in the potential flow. 


It has been shown by Sears‘ that the external flow 


over a rotating cylindrical blade is given by 


uy VReix(X, Zz) 
Yy = R[gi(X, Z) aie 2X | (14 
w, = VRgz(X, Z) 


/ 


where ¢i(X, Z) denotes the potential for plane, steady 
flow past the cylinder placed in a parallel stream of unit 
speed in the X direction. The assumption is made 
that the relative orders of magnitude of the spanwise 
and chordwise velocities in the boundary layer are de 
termined by the external flow. For example, since the 
blade was assumed to be thin when Eqs. (5), (6), 
(7), and (8) were obtained, the potential velocities cay 
be written: 


um = yR+ um", 1% = —xR + 0,” 


where the perturbation velocities 1,” and v7,” are small 


in comparison with yR. Therefore, v/m%, = O(x/y), 
or, taking the maximum value of x to be the chord, 


C, 
V1 My = O(c ¥) 
The assumption mentioned above is that 
v/u = O(v/m1) = O(c/y) 


Similarly, by examining the potential velocities (14), 
it is seen that 

uy, = O(u/y), uy = O(u/c), etc. 
The relative orders of magnitude of the terms of the 
boundary-layer Eqs. (9), (10), and (12) have been ex- 
The 
term wu, is arbitrarily taken to be of unit order of mag- 


amined and are shown in the following scheme. 


nitude. 


uu, + Vuy + wu, — 2vR — xR? = —(1/p)P, + vu., (3 
7 ») 
(> twyyy= €1) (e/s)* /s)* (1) (1) 
uvr + wy, + wr, + 2uR — yR? = —(1/p)py + w.. (16 
) 
ke/y) ‘te/y)* (c/¥) (c/y) (c/y) (c/y) (c/y) 
Uz + Wy Ww, = @ (17 
(it) te/y* ) 


For stations far out on the blade such that (cy)? < 1, 
the approximation is made that terms of relative order 
of magnitude (c/y)* can be neglected. 


uu, + wu, = —(1/p)p, + mu-:: (18) 
uv, + we, + 2uR — yR* = —(1/p)p, + w.. (19) 
p. = 0 (20) 


ur + wW, = O (21) 


The pressure gradients p,, p,, are calculated from the 
Bernoulli equation for this case :4 


(p/p) + (1/2) (mum? + 2,2) = 
(R2/2) (x2 + y*) + const. 


(22) 


(w, vanishes in this coordinate system. ) 
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LAMINAR BOUNDARY 
Pr/p = XR? — uj, — wz (23) 
py/p = yR? — mm, (24) 


The term 221; is of relative order of magnitude (c/y)?, 
as is the term xR*, which would cancel the similar term 
in Eq. (15) in any case. 


The equations of motion now read 


UU, + WU, = UyMyz + VUz2 (25) 
uv, + we, + 2uR = Mth + v2 (26) 
Ur +w, = 0 (27) 
The boundary conditions are 
u=v=w =), z=0 | 
u —> yR¢iz(x) t (28) 
v— Ri[gi(x) — 2x) > seo | a 
w—> 0 } } 


Notice that Eqs. (25) and (27), together with the 
boundary conditions (28) on u and w, are precisely the 
same as if the flow were two-dimensional with free- 
stream velocity yR in the positive X direction. There- 
fore, Eqs. (25) and (27) and their boundary conditions 
can be solved by any of the methods available for 
solving two-dimensional laminar boundary-layer prob- 
lems, and, with w and w known, the linear equation for 
vcan be solved separately. 

This behavior of the flow is similar to that found by 
Prandtl,® Jones,® and Sears’ for the case of the yawed in- 
finite wing. In particular, as for the case of the yawed 
wing, it is clear that rotation of the blade does not affect 
the position of the line of separation of the chordwise 
flow component from the surface. 

Notice also that, if the blade is being translated along 
the axis of rotation, the equations of motion are not 
changed, and the translation appears only in the 
boundary conditions. It has been shown by Fogarty 
and Sears* that the potential flow over a rotating cylin- 
drical blade that is being translated along the axis of 
rotation with speed W is given by 

uw = VYReix + Weox 
v1, = R(g, — 2X) (29) 
VYReiz + Werz 


where 9,;(X, Z) and g(X, Z) are the potentials for plane 
steady flow past the cylinder in a parallel stream of unit 
speed in the positive X and negative Z directions, 
respectively. If W/yR < 1, W/yR = O(c/y)? say, 
the effect of the translation may be neglected in this 
approximation. Thus, the solutions of Eqs. (25)—(27) 
will apply to rotating cylindrical blades and also to 
blades of constant cross section lying in the helix gen- 
erated by steady rotation and translation of the blade, 
provided that W/yR = O(c/y)’. 

In this paper, the application of the theory to the 
calculation of spanwise flow in the boundary layer will 
first be illustrated by means of two cases for which the 
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two-dimensional boundary-layer flow is known. These 
are, respectively, the flat plate and a certain symmet- 
rical cylinder. For these cases, the solutions of Eqs. 
(25) and (27) are available, and it remains only to con- 
struct solutions to Eq. (26). Following this, the treat- 
ment of more general cases will be discussed briefly. 


(II) THe RoratinG PLATE 


Consider a flat plate of chord c and infinite span rotat- 
ing steadily with angular velocity R about the Z axis. 
For this case, the potential flow is 


uw, = yR, y= —xR, w, = 0 (30) 
The equations of motion are 
UU, + Wu, = VUz, (31) 
uv, + we, + 2uR — yR* = w,, (32) 
ury+w, = 0 (33) 
and the boundary conditions are 
u—yR 
yv—> — xR> s—> @ 
(34) 
w>0 | 
u=v=w=0, 2:=0 } 


Eqs. (31) and (33) and the boundary conditions (34) 
on « and w were solved by Blasius,? who introduced a 
function F(r) such that 


u = yR[F'(r)/2] 


; — F(r) tF’(r) |] + (35) 
w= (Re) “yR —— + —=— | | - 
: ae 2a ) 
where 
rT =(1 2)[(Re)' *(z/c)/(x/e)' 
Re = yRc/v (36) 
¢ =x/c 


If the above expressions are substituted for uv and w 
in Eq. (32) and the assumption is made that 
v = xRG(r) (37) 
the following ordinary differential equation for G re- 
sults: 

G" + FG’ — 2F°G = 4(F’ — 1) (38) 

The boundary conditions are 
G= 0, 7 = 0 


G—> —-l, T—> © 


(39) 


Only tabular values of F and F’ are known, so it is 
necessary to solve Eq. (38) numerically over the range 
0 <7+< o. The tabular values of F’ show F’ & 2, 
4<7< ». We see that the solution of Eq. (38) 
which satisfies the boundary condition G—~> —1, 7 —> 
is 


Ge~=-l, 49°50 (40) 
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5 (III) FLow over A BLADE WITH THICKNESS T 
out | 
The boundary-layer flow over a certain class of cy- 
i ae 3 : ’ L 
lindrical blades of chord c whose two-dimensional po. 
4 tential flow velocity is of the form 
< + 
u, = ByR ——}, CS 2 fc 
c co 
3 ) ; ’ 
f is next considered. Here, Bisa parameter related to the whe 
thickness of the blade. This external velocity dis- 
$ tribution produces a “‘separation’’ type of profile of the 
2 - 7 chordwise velocity component in the boundary layer 
\ ce vk that is, a profile having negative values of u, at the “ai 
X Pa blade surface downstream of the separation line. '° 
a The boundary-layer equations for this case are 
' —< sii ' 
N Le 
te ea B?y? R2 /x x3 ; x? 
ol uu, + Wu, = — 1-3 =] + vue, (42) The 
Cc Cc c* c* ° 
—_ tio! 
r@) . -3.2 | (& 
: ere x } fs 
9 2 4 & 8 1.0 uv, + we, + 2uR = B?yR? ( aon a + w., (43) | 
VELOCITY FUNCTION —~ . ' 
Fic. 2. Velocity profiles in the boundary layer of a rotating : 
flat plate. v, and ve denote the radial and tangential components ur + w, = 0 (44) | 
of the velocity relative to the plate. x denotes the distance from Su 
the leading edge, y the radius, and R the angular velocity. ¢ is The boundary conditions are computed by the - 
the dimensionless distance from the plate, defined in Eq. (49). ‘ ae . g il 
Note: For Vr/XR and Vo/YR in this figure, read v,/xR and method given in reference 4. ‘ 
v,/yR, respectively. ing 
x sa ‘ 
: j uw, = VR giz = BYR — (45) Gy 
We therefore approximate to the actual boundary con- 1= YR Ou = Be. & =) 
dition G— —1, 7 — © by taking G = —1,7 > 4, and 
~ 90 . 2 -4 
solve Eq. (38) numerically over the range 0 < 7 < 4. = io B e Bx 2x) (46) 
ne : “ : : “y= gi — <x) = PSR el Si ae 
The equation for G was integrated numerically over 2 tc8 wl 
the range 0 < +r < 4. Since the calculated velocities The ' hit} 
° “ P ° > aryc 1S are 
are referred to the rotating Cartesian coordinate sys- ihe hountary commitions are 
tem, the results are apt to be confusing to an unwary cuwemeonut’é — 
reader. For greater ease of interpretation, plots of . ¥! 
radial and tangential velocities, rather than spanwise u—> ByR (“ = =) 
and chordwise velocities, are presented. If the angle Bx? Bx" a (47) | rl 
between the y axis and the radius vector to the point v—>cR & -_ a -_- ) ; 
x, yon the blade is small, - 4c : 
w—> 0 
v, Av+ (x/y)u i 5 F , 
are rhe two-dimensional problem represented by Eqs. 
ve 4° pa 
(42) and (44) and the boundary conditions (47) on u i. 
The quantities v,/xR and v/yR are plotted in Fig. 2. and w was originally solved by Blasius,’ who introduced I 
Also of interest is the direction of the flow near the the stream function W to satisfy Eq. (44) and obtained q 
surface. If @is the angle of the limiting streamline with 4 power series solution of Eq. (42). Following Gold- } fi 
the x direction, stein,* let th 
(! 
F v : xRG a . se sal ) 
tan@é = lim = lim a uy, = Bix + Bax? + Box Fees | (48) q 
soO0% +0 9RF’/2 YW = Fix + Px? + Be +... 
This limit has been evaluated numerically. where, in this case, 
tan 6 = (2x/y) (1.323/1.328) (41) ByR ByR 
‘Sa = 3 = = @eP R=... 89 
The writer suggests that it is more than a numerical c 
coincidence that and the F's are functions of z. f 
hm. (G/P') = 1:00... One substitutes u = V,, w = —W, in Eq. (42) and p 
‘gia The resulting set 


however, he has not yet been able to demonstrate that 
this fact has any significance. 


equates coefficients of powers of x. 
of ordinary differential equations for the F’s has been 
solved, and the results are tabulated in reference 3. 
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The solution of Eq. (43) for v is conveniently carried 
out in nondimensional form. 


Let 
x y eee 
g _ ’ U] - -, ¢ -_ Re 
c c ’ i (49) 
t v Ww 
f= ’ ,= ’ h = Re’ 
J ByR . BcR ByR ) 
where Re = ByRc/v. Eq. (43) becomes 
2f ; 
Sg: + hg; = 9 8 = (§ — $9)? + Ley (50) 
and 
f _ filé = Af,’ &8 + 6hs’ ein te fe es (51) 


h = —f; + 12f3¢? — 30h5é4 + a a 


eg : ’ a 
The functions fi, fi’, fs, fs’ . . 
tions of ¢ in reference 3. 


. are tabulated as func- 
A power series expansion for 
g(, ¢) is assumed: 


co 


2(& $) = DY &"G, (6) 


n=O 


(52) 


Substitution of expressions (51) and (52) for f, h, and 
gin Eq. (50) yields an equation in powers of —. Equat- 


ing coefficients of powers of &, 
G,” + fiG,’ — nfi'Gy, — 12fsG, - + 
A(n — 2)fs'Gr_2 + 30h5G,_4’ — 


G(n — 4)hs’Gn4 +... = Qn (53) 
where 
Oo = 0 Os = 2 
—; = 2f,' 'B Os = 12hs'/8 
Q» =—| Os = —]1 
Q; = —S8f;'/B 
The boundary conditions are 
G, = 0 all n, ¢=0 
Gy > 0 G3 > 0 ee (54 
G,—> —2/8 Gym —1/4 aoe wore 


The tabular values of f,’, fs’, hts’, . . 
quantities are 


. show that these 


nearly constant and the quantities 


If 
these approximate expressions are inserted in Eqs. 


(53), 


fi, fs, As, . . . are linear functions of ¢ for ¢ 2 4. 


it is possible to obtain analytical solutions for the 


quantities G,, which are valid ¢ 2 4. The results are 


Gy. = 0 G, = —1/4 

G, = —? B Gs = 0 =~ 
’ . (2) 

Ge _ l 2 Ge = Q 

G3; = 0 


Eqs. (53) can now be integrated numerically over the 
finite range 0 < ¢ < 4. 
carried out for Gy through Gs. 


The integration has been 


he profile of the spanwise flow was calculated for 
8 = 3, which gives a maximum value of u,/yR of 1.155. 
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Fic. 3. Velocity profiles in the boundary layer of a rotating 
symmetrical cylinder for which u, = 3yR (€ — &). Here, & de- 
notes x/c, cis a reference chord, and other notations are the same 
as in Fig. 2. Nore: For Ve/CR and Vr/YR in this figure, 
read v,/cR and v,/yR, respectively. 


The corresponding value for an NACA 0010 airfoil is 
about 1.158. The potential speed, according to the 
formula used here, 


u, = ByR(E — €) 
exhibits a stagnation point at & = 0, rises to the maxi- 


mum value, and decelerates. 
the flow around 


Thus, it is typical of 
the forepart of an airfoil at small or 
zero incidence. Curves of radial and tangential ve- 
locity profiles are presented in Fig. 3. 

It is seen that there is no appreciable radial flow near 
the separation line for this case. For larger values of 
8, there would be a larger pressure gradient along the 
line of chordwise flow separation, and greater spanwise 
flow in the boundary layer might occur. One must 
remember, however, that the present theory is limited 
to slender shapes. 


(IV) PROCEDURE FOR ARBITRARY CYLINDERS 


The method of solution used here—namely, expan- 


sion in powers of x—has been particularly convenient 
because the two-dimensional solution was available in 
that form. This will not be the case in general, and a 
more convenient procedure will be needed for calcula- 
tion of both uw and v for arbitrary u(x). 

One of the practical methods of approximate solu- 
tions of the boundary-layer equations is that 


von Karman and Pohlhausen.!° 


of 
A simpleex tension 


of this method can be used in the present case. If an 
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analysis similar to that of reference 11 is made, a first- 
order ordinary differential equation is obtained for the 
ratio of the spanwise and chordwise boundary-layer 
thicknesses. While this method is simple in theory, 
it becomes laborious in practice, because of the large 
number of terms in the differential equation which must 
be solved. An approximate check of the validity of the 
method for the case of the blade with thickness treated 
above has been made. It is concluded that the method 
should yield a good approximation in this typical case. 
(V) CONCLUSIONS 

The striking feature of the results of the two cases 
investigated is that the effects of rotation as calculated 
are small. This is in disagreement with the engineer's 
common expression that three-dimensional effects on 
rotating blades are large. Possible explanations of this 
disagreement are: 

(1). The engineer's impression may be derived from 
observation of what happens after separation of the 
flow from the surface. The theory presented here is not 
valid beyond the separation line and so cannot predict 
such effects. 

(2) The profiles treated have weak pressure gra- 
dients. Comparison of the results of references 7 and 
11 shows that, for these two cases of yawed wings, the 
profile with strong pressure gradient had much stronger 
spanwise flow than the profile with weak pressure gra- 
dient. It is reasonable to assume that some such effect 
would be obtained with rotating blades. 
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(3) It is possible that the effect of finite span—that 
is, of a tip or wall—is to accentuate the effects of blade 
rotation locally. 

(4) Only the laminar boundary layer is considered 
here. It is possible that the effects of rotation on the 
turbulent boundary layer are more profound. 
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First U. S. National Congress of Applied Mechanics 


Plans have been formulated to hold the First U.S. National Congress of Applied Mechanics on 
The Illinois Institute of Technology will be host to the Congress. 

Members of the Institute, one of the sponsoring societies, are cordially invited to submit papers, 
| which should constitute original research in applied mechanics. 

words or the equivalent in equations, tables, and diagrams.) Abstracts must be submitted to the Chair- 

| man of the Editorial Committee before April 14, 1951; complete papers must be submitted prior to June 
| 
| 


Additional information about the Congress may be obtained by writing to N. M. Newmark, 
Secretary of the Congress, University of Illinois, Urbana, Ill. 


(Their length must not exceed 5,000 


All members are invited to attend. 
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A Flexible Nozzle for a Small Supersonic 
Wind Tunnel’ 


SATISH DHAWAN? ann ANATOL ROSHKOt 
California Institute of Technology 


SUMMARY 


The design of a small supersonic wind-tunnel test section (4 by 
The flexible 


10 in.) incorporating a flexible nozzle is outlined. 
Two screw 


nozzle consists of a high-strength stepped steel plate. 
jacks provide an easy means of continuously changing the nozzle’s 
shape according to the aerodynamic requirements. The bound- 
ary-layer compensation can also be varied during operation. 
Pressure surveys, together with schlieren and interferometric 
analysis of the test section, show the flow to be uniform over the 


operating range (7 = 1.1 61.5). 


INTRODUCTION 


an DEVELOPMENT OF experimental research in 
transonic and supersonic flow has increasingly 
demanded the use of a wind tunnel with continuous 
control over the flow parameters. Since the air speed 
in a supersonic wind tunnel is determined completely 
by the geometry of the nozzle, it is necessary to pro- 
vide a different nozzle shape for each different required 
speed. The use of a flexible nozzle permits a continuous 
change in the flow Mach Number and constitutes a 
considerable advantage and improvement over the 
alternative of using a series of fixed geometry nozzles. 
The idea of using a flexible nozzle is not new, and there 
are in existence and operation in this country several 
supersonic wind tunnels incorporating the device. 
However, as far as can be determined from published 
descriptions of these wind tunnels, the flexible nozzle 
with its shape-changing mechanism is rendered com- 
plex by the range and accuracy requirements. For 
instance, in order to change the nozzle shape the wind 
tunnel has to be shut down and the required changes 
made in small degrees of adjustment. During the past 
year or so there has been in use a flexible nozzle of ex- 
ceptionally simple design in the 4- by 10-in. Transonic 
Tunnel at the California Institute of Technology.f 
The supersonic range of the tunnel is limited from 
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t The test section described here is a modification on the orig- 
inal test section of the 2- by 20-in. high subsonic tunnel built in 
1944 under A.A.F. contract for high-speed flow investigations 
A study of boundary-layer and shock-wave phenomena is at pres- 
ent in progress under N.A.C.A. sponsorship; the design pre- 
sented here was executed in the course of this work. Results of 
recent research in the tunnel are being published elsewhere by the 


N.A.C.A, 


approximately Mf = 1.1 to 1.5, and this fact is fully 
exploited in the interests of simplicity to achieve a 
design that permits continuous shape changes in the 
nozzle during operation while retaining accuracy and 
uniformity of flow. Over the major portion of the 
operating range, the only operation required to change 
the flow Mach Number is the turning of a single jack- 
screw. This paper briefly outlines the interesting fea- 
tures of this tunnel with special attention to the design, 
construction, and use of the flexible nozzle. 


DESCRIPTION OF TEST SECTION 


The working section of the GALCIT 4- by 10- by 
48-in. transonic tunnel is sketched in Fig. 1 showing the 
essential features of the design. The floor block of the 
test section carries the one-wall flexible nozzle plate, 
together with the swiveling jackscrew controls. The 
floor is hinged just downstream of the main jack, and 
its downstream end can be raised or lowered by means 
of a screw jack. This allows alterations in boundary- 
layer compensation during operation. The ceiling 
block of the tunnel supports the traversing mechanism 
and contains a narrow slot for operation of the travers- 
ing arm. A pressure box, with access panels mounted 
on the ceiling block, encloses the traversing mechanism 
and seals it to the test section. Pressure-sealing of the 
test section is secured by means of rubber tube-in- 
groove seals between the sidewalls and floor and the 
ceiling blocks. The sidewalls are made in panels and 
allow the window to be located anywhere in the 48-in. 
length of the test section. The main flexible nozzle 
consists of a continuous spring steel plate 45'/, in. long 
of varying thickness, of which 36 in. forms the flexible 
nozzle and the remainder is the floor of the test section 
(see Fig. 1). It is anchored in the contraction with 
the downstream end also direction-fixed, but it is free 
to move horizontally on rollers when deflected by the 
jacks. A second flexible nozzle plate begins where 
the primary nozzle ends. The flexible second throat 
acts as a speed control during subsonic operation of 
the tunnel and may be used as an adjustable super- 
sonic diffuser during supersonic runs. 


DESIGN CONSIDERATIONS 


(1) Aerodynamic Nozzle Shapes 
The supersonic range of Mach Numbers in the tunnel 


is from IJ = 1.1 to 1.5. Fig. 2 shows the theoretical 
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nozzle shapes covering this range. These shapes were 
obtained graphically by the well-known method of two- 
dimensional characteristics.'. The fundamental pro- 
cedure in this method consists of assuming an initial 
expansion curve starting from zero inclination at the 
throat up to the nozzle expansion angle 0 (Fig. 3), 
which must be less than, or at most equal to, half the 
Prandtl-Meyer angle associated with the final Mach 
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10- by 48-in. transonic-tunnel test section 


The initial expansion curve is then divided 
The expansion wave- 


Nutnber. 
into suitable straight segments. 
lets generated by this convex initial curve determine 
uniquely the subsequent form of the nozzle for parallel 
flow at the exit. For a required Mach Number in the 
test section, there is thus an infinity of possible initial 
curve and expansion angle 8 combinations. Theoret; 
the shortest length nozzle would be obtained by 
using the maximum value of 0 (equal to half the 
Prandtl-Meyer angle) and an initial expansion of in- 
finite curvature—i.e., a sharp corner at the throat. 
Actually, the layer 
limits the expansion. 
radius of curvature at the throat of approximately two 


ically, 


however, presence of boundary 


rate of For design purposes, a 





to four times the test-section height may be used. In | 
choosing the aerodynamic shapes shown in Fig. 2, the 
initial expansion wave and the nozzle slope 0 were so 
adjusted as to give constant length nozzles (throat to 
This per- 


exit) for the various design Mach Numbers. 
mitted the use of a constant length flexible plate for the 
supersonic nozzle and simplified the mechanical design 
considerably. 


2) Boundary-Layer Compensation 


The theoretically obtained nozzle shapes do not take 
into account the growth of the boundary layer on the 
walls of the tunnel. In order to allow for this, it 1s 
usual to displace the physical walls at each point by the | 
“displacement thickness”’ of the boundary layer at that 
point. One method of obtaining approximate estt- 
mates of the boundary-layer growth is to assume the 


pressure gradient on the walls as obtained by the in- 
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— 
viscid flow computation and, then, to use boundary- 
laver theory to calculate the displacement thickness. 
in the design under discussion, the compensation for ‘ =. (a) 
boundary layer was taken into account in a more ele- 3 Ps mea rong 
mentary fashion. Experience in the previous test \* Pf 
section of the wind tunnel had shown that a linear cor- \ Pa 
rection of 0.021 in. per in. in the walls was sufficient to \ L-7 
allow for the boundary-layer growth. In addition, the — a ® 7 (b) 
present design includes provision for alterations in the 4 7 ’ SIMPLE BEAM WITH 
allowance with the tunnel in operation. Subsequent \ |" TWO LOADS 
results justified the estimates used. f 3 W Ps 
\ \ J 
, y 
Tue Nozz_e PLATE : ¥ / 
, \ ra 
TE 1) Requirements / \ 
(a) The nozzle plate should be able to reproduce, : p 
‘af during operation, any aerodynamic nozzle shape in the \ Fi a 
required range. . : Fi 
. (b) The shape-changing mechanism should be t ue oo Re omection 
simple and nozzle settings accurately repeatable. \ Fa - FIXED ENDS AND 
, \ y, 2 TWO LOADS 
(c) The plate must be free from local distortion : VIM) 
and vibrationally stable during operation of the wind ; \ y 
tunnel. \ / 
\ PORTION OF BEAM 
2) Plate Shape and Control Configuration ;: -| APPROXIMATING |. 
vided ‘ . ; - \ ITHE AERODYNAMIC 
Study of the aerodynamic shapes (Fig. 2) to be re- \ NOZZLE 
ore produced by the nozzle plate shows that, in general, 4 : _————} (d) 
— the aerodynamic shapes have an initial region of rela- ; ; "STEPPED" BEAM WITH 
rallel tively high curvature followed by an almost straight \ 7 DIRECTION FIKED ENDS 
» te - ge ; a tas \ | | /| AERODYNAMIC AND TWO LOADS 
Bs. portion containing an inflexion point (reversal in sign \ '| (M=1.365) 
utial of the curvature). The terminal section of the shapes \| IV 
oret; is seen to have relatively the largest length and a uni- \ 
d by formly decreasing curvature. Fig. 4 shows the curva- / 
P Z tures of various configur ations of a plate represented as Wes. 4. Canvetuses of clastic and anredyanmie shapes. 
an elastic beam in relation to the curvatures of a repre- 
ans sentative aerodynamic shape. It is clear that a mini- 
my ocean is required to produce the required and having direction-fixed ends. The location of the 
— = in sign of curvature in the plate. Fig. 4 also jacks and the step in the plate was determined by trial 
2 shows how the jacks for changing the shape of the to give the best overall reproduction of the aerodynamic 
ps nozzle should be located. Since the greatest curvature shapes over the required range. 
, in the plate is produced in the vicinity of the main jack, 
adie this must be situated close to the throat of the nozzle. 
oie It is also evident that it is impossible exactly to repro- ape 
per | duce the required curvatures with a continuous plate 7 
r the | of constant thickness. Fig. 4d shows the manner in J a [== ‘ = g 
Sign which a variation in moment of inertia of the plate a. wes ‘riTrmg ; 7 Ouse anu ROLLERS 
~— $TUDS LTS 
cross section may be used to advantage in bringing the rvs 
elastic and aerodynamic shape curvatures closer to- 
gether. The second integrals of these curves—i.e., the NOZZLE PLATE AND FITTINGS 
take actual shapes—would then be still closer. In order to x 
: = secure smooth entrance and exit conditions at all Mach > 2x 
it iS | Numbers, it is highly desirable to have the ends of the ooss" 3 018s" — as 
the nozzle plate fixed in direction. A comparative study 4% 7 -— ty rz] r™ 
that of several configurations along the lines indicated above ay eee eee is" RE o— 
est FT led to the adoption of the configuration schematically lf 36" + 
fs represented in Fig. 5. The nozzle plate consists of a 
» in- 


REPRESENTATION AS FIXED-END BEAM 


Stepped beam loaded at two points by screw jacks Fic. 5. 
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(3) Determination of Control Settings 


The restriction on the deflection curve of the plate 
when it is required to reproduce an aerodynamic shape 
is that it must have a prescribed maximum deflection. 
This is given by the unique height ratio h/h* associated 
with a desired supercritical flow in the test section. 
Here, / is the height of the two-dimensional test sec- 
tion, and h* is the height at the throat. Knowing the 
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dimensions, end conditions, and type of loading on the 
plate, the elastic influence curves of the configuration 
can be calculated. These permit a determinatiop 
of all the possible plate shapes (Fig. 6) with a givey 
maximum deflection in terms of the load ratio » = 
W:/Wi, where W; = load in pounds on the main jack. 
screw and W2 = load on the second jackscrew. 4 
superposition of the required aerodynamic shape and the 
corresponding plate curves serves then to determine 
the optimum value of v, and the control settings are 
computed from the deflection influence functions of 
the plate. A representative case of matching of the 
aerodynamic and elastic curves for M = 1.5 is shown in 
Fig. 6, the vertical scale beimg' exaggerated to increase 
the accuracy of matching. Fig. 7 shows. the final 
curves from which control settings for any Mach Nun- 
ber within the design range may be obtained. It will 
be observed that over almost the entire range of 
operation only the main jack setting is required to be 
altered. 


FABRICATION OF THE NOZZLE PLATE 


The problem in fabrication was to obtain a relatively 
thin, high-strength plate free from local deformation, 
particularly local deformations near the attachment 
fittings. 
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dbs 
Schlieren of flow past a cone. 
{ = 1.38. 





Fic. 9. Cone half angle = 15°. 


(1) Choice of Material 


From simple beam theory for a beam with iixed ends 
and given deflection curve and thickness distribution, 


stress ¢ ~ Et 


deflection load W ~ Et* 
where 
E = modulus of elasticity of the beam material 
t = some reference thickness of beam 


For 
must be kept as high as possible. 
between iron and copper alloys shows that, for a given 


the plate to be stable during operation, WV and / 
A typical comparison 
factor of safety, 

bsteel tbronze + 1.5 
and 


WW steel HW bronze = ¢ 


These considerations and problems of welding and 
heat-treatment dictated the choice of steel. SAE 
6150, an alloy spring steel, was found to meet the re- 
quirements. 


(2) Machining of Plate 


“step” and stiffening ribs was ma- 
by 5-in. piece of annealed SAE 6150, 


The plate with 
chined from a */s- 
the studs and fittings also being machined from the 
same material. About 0.05 in. 
left on the upper surface of the plate for final surface 


of extra material was 
grinding. 
(3) Attachment of Fittings 


This was done by copper brazing in a hydrogen at- 
mosphere furnace at 2,000°F. (above the heat-treat- 
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ment temperature). During this operation the plate 
was clamped in a steel jig. 


(4) Heat-Treatment 


The plate was heat-treated to give a Rockwell hard- 
ness of C-38, corresponding to a yield of 154,000 Ibs. 
per sq.in. During the normalizing, quenching, and 
tempering processes, the plate was clamped in a jig. 
The results obtained extremely satisfactory. 
There are no local distortions ia the plate. 


were 


(5) Surface Grinding 


The working surface of the plate was ground to a 
polish finish, eliminating all surface irregularities. 


(6) Stress Test 


Since the nozzle plate was to carry extremely high 
stresses in its fully deflected state, a strain-gage test was 
made to ensure safe operation. The maximum stress 
produced in the plate when set for MJ = 1. 
to be 110,000 Ibs. This gives 


safety of 1.4. 


5 was found 


per sq.in. a factor of 


EVALUATION AND CALIBRATION 


Figs. Sa and 8b show the results of pressure surveys 
of the test section made at the various design Mach 
Numbers. The vertical distributions, Fig. 8b, reveal 
the effect of waves originating at the nozzle and indi- 
cate uniform flow with negligible gradients over the 


entire operating range. Fig. 9 is a schlieren photo- 
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Wedge half-angle 


Interferogram of flow past a wedge. 


= 45°. M = 1.37. 
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graph, and Fig. 10 is an interferogram of the flow in the ) 
0.48 test section. The straight head wave in Fig. 9 and the 
SECOND THROAT CALIBRATION straight interference fringes in Fig. 10 confirm the uni- 
0.44 FOR form wave free flow indicated by the pressure surveys, 
SUBSONIC OPERATION The somewhat greater variations in the horizontal 
0.40 - | ] surveys (Fig. 8a) are due to the small disturbances 
introduced at the joints in the sectioned sidewalls. 
ed km The use of continuous panels for the sidewalls reduces 
F the variations in the horizontal direction to the same 
wW . . . . 
032} order as in the vertical surveys (see Fig. 8a). Figs, 7 
§ and 11, respectively, show the calibration curves for the Dr 
0.28 5 controls over the supersonic and subsonic range of oper- supel 
t ation of the tunnel. The change from supersonic and | 
0.24 2 to subsonic flow or vice versa can be accomplished sphet 
S| without shutting down the tunnel. _ 
o20}-4 | but 
| show 
CONCLUSION tion 
| 
0.16 t F . ” r . chan, 
A flexible relatively high-strength plate free from a 
a | local distortions has been successfully designed and used Flow 
| | | for the production of uniform shock-free flow over the actio 
pre | | | | | range M = 1.1to1.5. The plate is vibrationally stable ve 
a : : coe ie Ad 
| | during operation, and the repeatability of flows is ex- inne 
| | mpc 
> a ce cellent. mt 
TEST SECTION MACH NUMBER | 
| REFERENCE 
0 
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Sphere Drag Data at Supersonic Speeds and 
Low Reynolds Numbers’ 


E. D. KANET 
University of California at Berkeley 


SUMMARY 


Drag forces on spheres have been measured in a low-density 
supersonic wind tunnel. Mach Numbers were from 2.1 to 2.8, 
and the Reynolds Number (based on free-stream properties and 
sphere diameter) ranged from 15 to 800. Available experi- 
mental data for the same Mach Numbers as the present tests, 
but with Reynolds Numbers in the order of magnitude 105, 
showed that the drag coefficients could be correlated as a func- 
tion of Mach Number alone and were not altered by a sixfold 
change in Reynolds Number. The present data showed an in- 
crease of drag coefficient by a factor of 2'/2 over the 
Flow visualization photographs indicated a possibility that inter 


test range. 


action occurred between the sphere boundary layer and shock 
wave at low Reynolds Numbers. 

Additional experimental data are required before the relative 
importance of viscous, molecular flow, and shock-wave effects 
can be evaluated completely. 


SYMBOLS 


R = Reynolds Number = Vdp/u (dimensionless) 

R; = gas constant, defined by Rg = p/pT 

V = velocity (consistent units) 

d = sphere diameter (in.) 

p = gas density (consistent units) 

be = gas viscosity (Ibs. per sec. ft.) 

( = average gas molecule speed (consistent units) 

p = gas pressure (microns Hg) 

y = ratio of specific heats = 1.40 for air (dimensionless ) 

a = local sound speed in the gas (consistent units ) 

M = Mach Number = V’/a (dimensionless) 

l = molecular mean free path [defined by Eq. (1)] (in.) 

A = an area, defined by Eq. (4), = md?/4 for a sphere 
(sq.in.) 

Ky = drag coefficient for a sphere [defined by Eq. (5)] (di 
mensionless ) 

D = drag force on sphere (Gm. ) 

r = gas temperature (°F. abs.) 

fr) = boundary-layer displacement thickness (in 

a, B = constants (dimensionless) 

Subscripts 

¢ = compressible, nonviscous fluid [Eq. (12) ] 

1 = impact (pressure ) 

§ = static (pressure ) 

l = conditions behind normal shock wave 

2 = quantity measured at surface of cone (Probe No. 15) 


INTRODUCTION 


Mes REMENTS OF DRAG FORCES on spheres by 
many workers have served to clarify certain 
Presented at the 1950 Heat Transfer and Fluid Mechanics 
Institute, Los Angeles, June 28-30, 1950. 
* This investigation was aided by funds supplied under a 
contract from the Office of Naval Research. 


1 Associate Professor of Engineering Design. 


aspects of the general problem of resistance of bodies to 
uniform motion in a fluid. Most of the available in- 
formation for sphere drag forces at subsonic speeds has 
been summarized in reference 1. The data are corre- 
lated as a function of Reynolds Number over the range 
from approximately 0.04 to 10°. The lowest Rey- 
nolds Numbers were obtained for spheres immersed in 
liquids. Also available are the data of reference 2 for 
low Reynolds Numbers in air, under conditions that 
resulted in molecular flow effects. For supersonic 
speeds, the only published data appear to be those of 
Thomas* obtained by firing spheres 
The Mach Number range was ap- 
The Reynolds Number 
was from 


Charters and 
through still air. 
proximately from 0.3 to 4.0. 
variation, for Mach Numbers greater than 1, 
2X 10° to 1.3 X 10°. 
related by the Mach Number within the accuracy of 
experimental measurements, and there was-no Rey- 


The supersonic data were cor- 


nolds Number effect. 

The theory for the drag force on a sphere in uniform 
motion is incomplete. For incompressible fluids and 
low speeds, the analysis of Stokes,‘ obtained by neg- 
lecting inertia terms in the flow equations, is appli- 
cable. Comparison with experiment indicates that the 
analysis is valid for R < 0.5 and gives values too low 
at higher R. Oseen’ obtained a second approximation 
by considering the effect of some of the inertia terms, 
and Goldstein® computed a result from Oseen’s equa- 


tions which everywhere lies above the experimental 


data. This latter solution is within 5 per cent of experi- 
ment at R = 1.0 and deviates by more than 20 per cent 
at R = 10. Stokes’s equation can also be applied for 


spheres moving in a gas at sufficiently low speeds to 
give R < 0.5, 
of the gas is small compared to the sphere diameter. 
In Millikan’s*® experiments, for example, the data did 
not depart from Stokes’s equation until the ratio of gas 
molecular mean free path to characteristic body di- 
exceeded 


provided the molecular mean free path 


mension—in this case, the sphere diameter 
approximately 0.001. 

It is convenient to discuss the various results in 
dimensionless combinations of variables 


Karman’ pointed out 


terms of the 
used in fluid mechanics. 
that the Reynolds Number could be written 


R = (V/c) (d/l) 


von 


(constant) 


The first ratio is proportional to the Mach Number 


M = V 


a, since 
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Fic. 2. Drag balance arrangement. 


¢ = V(8/n) (p/p) = V8/ry a 


A result available from kinetic theory* is 


uw = 0.499 pel (1) 
It follows that (for air with y = 1.40) 
l/d = 1.49 (M/R) (2) 


Eq. (2) relates the parameter //d to the M and R com- 
monly used in fluid mechanics. 

Stokes’s law was extended by Basset,’ who used a 
modified boundary condition to account for ‘“‘slip”’ 
and thus obtained a correction term involving \//R. 
The resulting expression is plotted in Fig. 1. For 
comparison, an empirical equation that fits Millikan’s 
experimental data is included [the constants have been 
adjusted to conform with the value of / given in Eq. 
(1)]. It appears from Fig. 1 that Stokes’s equation 
modified for ‘‘slip’’ begins to deviate from the data at 
M/R = 0.02, while Millikan’s formula indicates that a 
slip correction of the form 
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1/(1 + Const. M/R) 


fits the data for values of 1//R up to about 0.12. 

For spheres moving in a gas under large mean free 
path conditions—i.e., / > d—theoretical expressions 
for the drag have been deduced by using the concepts 
of kinetic theory. Tsien’ has reviewed the work avail- 
able up to 1945, and some later papers!!—!* have ey. 
tended the calculations. 

Tsien” discussed the general problem of the effects 
to be expected at low gas densities and suggested cri- 
teria for estimating the boundaries of continuum, large 
mean free path, and transition flows. 
vestigation was concerned with measurement of sphere 


The present in- 


drag forces at supersonic speeds, over a Reynolds 
Number range that was expected to produce flows in 
the continuum and transition regions. 


EXPERIMENTAL 
PROCEDURE 


DESCRIPTION OF EQUIPMENT AND 


No. 3 Wind Tunnel 


The required low density supersonic gas stream was 
provided by the No. 3 Wind Tunnel at Berkeley." A 
single nozzle, designated the No. 2 nozzle, was employed 
in these experiments and gave a Mach Number range 
of approximately 2.1 to 2.8 and a static pressure range 
from 30 to 128 microns Hg. The free-stream tempera- 
ture (220° to 300°F. abs.) corresponded to assumed 
adiabatic expansion from an approximate stagnation 
condition at room temperature (70°-80°F.) to the indi- 
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Fic. 4. Probe No. 15. 
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cated Mach Number at the test section. The dimen- 
sions of the No. 2 nozzle are tabulated in reference 15; 
a photograph of the nozzle installed in the wind tunnel 
is shown in Fig. 2. 

The properties of the gas flow at the test section 
(where the drag forces were measured) were deter- 
mined by traversing with two probes, Nos. 14 and 15, 
in the test region and measuring the pressure asso- 
ciated with each probe (Figs. 3 and 4). The probes 
were connected by tubing to a Pirani gage and an oil 
manometer. The Pirani gage was a commercial type 
with platinum-nickel-alloy filaments and was inserted 
in an aluminum block that had been wound with re- 
sistance wire and covered with insulation. A copper 
constantan thermocouple was inserted alongside the 
Pirani tube in the aluminum block and connected to a 
potentiometer. By controlling the power to the re- 
sistance winding, the Pirani-tube wall temperature 
could be kept constant, as indicated by the thermo- 
couple, within +2°F. The power supply for the 
Pirani gage and the details of an arrangement for re- 
cording the Pirani-gage output on a Leeds and Northrup 
recorder are described in reference 16. The oil ma- 
nometer is a standard unit employed for low-pressure 
measurements at Berkeley. One leg of a U-tube 
manometer filled with butyl phthalate was evacuated 
by a diffusion pump backed by a mechanical pump. 
The pressure to be measured was connected to the other 
leg of the U-tube, and the difference in elevation between 
the menisci in the two legs was measured with a 
micrometer screw arrangement. Both the manometer 
and Pirani gage were calibrated frequently against a 


McLeod gage. 


Drag Balance 


The arrangement used for measuring drag forces is 
shown in Figs. 2 and 5. The sphere was mounted on a 
support rod of the size indicated on Fig. 6, the support 
rod standing vertically in the wind tunnel. The sup- 
port rod in turn was attached to the beam of an an- 
alytical balance. The brass support column of the 
original balance was changed to a plastic in order to 
insulate the upper part of the balance, and a motor 
drive was added so that the mechanism for raising and 
lowering the beam on its knife edges could be remotely 
controlled. With the balance located at the test sec- 
tion, with only the support rod and sphere extending 
into the flow region, the combined drag force on the 
sphere and support rod resulted in a moment about the 
main balance knife edge (Fig. 7). This moment was 
resisted by a quartz spring with one end hooked to one 
side of the balance beam and the other end hooked to a 
rod whose vertical position could be changed by turning 
a micrometer screw. The motion of this rod was also 
motorized for remote control, and the angular move- 
ment of the micrometer screw was geared to a counter 


that could be read at a distance. The gearing was 
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Fic. 7. Schematic of drag balance. 


arranged so that the counter indicated the vertical 
position of the bar directly in thousandths of an inch. 
In operation, the quartz spring could thus be extended 
a known amount until the moment due to drag on the 
sphere and support was balanced by the moment in- 
troduced by the spring force, and the balance beam 
thus returned to a null position. The null position was 
indicated by the original balance pointer and scale, 
viewed at a distance through a small telescope. In 
operation, the balance arm was frequently raised and 
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TABLE 1 
Summary of Experimental Results 
Flow Axial 
Sphere Rator Position * - — Pan Force 
Diameter (Scale (In.) (Microns (Milligrams) 

Run (In.) Reading) (+) Hg) M Total Tare Net Kp R 
2 1.0 § 80 0.10 93 2.40 2,812 2,488 0.480 464 
41 Lo. F 79 0.10 98 2.31 324 aS 0.494 448 
42 Los 139.5 0.10 128 2.56 4,280 : 3,759 0.474 753 
41 Lo FT 140 0.10 126 2.53 ; 521 0.482 768 
44-5 0.258 80 0.10 106 2.30 337 180 0.545 120 
47 0.26 T 80 0.10 86 2.55 157 ; 0.547 126 
47 0.25 T 80 0.10 100 2.33 157 0. 564 116 
44-5 0.258 8O 0.40 106 2.30 340 182 0.550 120 
47 0.25 T 80 0.40 86 2.55 158 0.551 126 
45 0.258 30 0.10 42 2.57 179 107 0.645 63 
47 0.25 T 30 0.10 56 2.24 72 0.648 59 
47 0.25 T 30 0.10 63 2.10 72 0.650 57 
45 0.258 30 0.40 42 2.57 181 108 0.6638 62 
47 0.25 T 30 0.40 56 2.24 73 0.656 58 
46 0.25S 18 0.10 31 2.45 121 76 0.695 4] 
47 0.25 T 18 0.10 31 2.46 45 : 0.689 $] 
47 0.25 T 18 0.10 44 2.08 15 0.680 39 
46 0.258 18 0.40 31 2.45 120 75 0.688 4] 
47 0.26 T 18 0.40 31 2.46 : 45 0.682 l 
47 0.25 T 18 0.40 45 2.05 45 0.677 38 
49 0.108 18 0.10 46 2.12 19.0 19.4 0.995 17 
51 0.10 T 18 0.10 36 2.23 29.6 1. ii) 15 
49 0.108 18 0.40 33 2.50 48.8 18.8 0.967 18 
51 0.10 T 18 0.40 38 2.40 30.0 1.13 4 
50 0.108 80 0.40 77 2.65 153 36 0.698 50 
51 0.10 T 80 0.10 94 2.34 117 0.738 45 
50 0.10S 8O 0.40 94 2.40 154 36 0.702 47 
51 0.10 T 80 0.40 93 2.37 118 0.726 45 
50 0.1085 30 0.10 51 2.37 74.3 23.6 0.873 5 
51 0.10 T 30 0.10 56 2.20 50.7 0.927 22 
50 0.108 30 0.40 50 2.42 74.7 23.7 0.857 26 
51 0.10 T 30 0.40 56 2.20 ° 51.0 0.926 22 
33 0.258 140 0.10 107 2.81 496 255 0.514 204 
54 0.25 T 140 0.10 113 2.71 241 0.522 194 
53 0.258 140 0.40 112 2.7 497 255 0.513 199 
54 0.25 T 139.5 0.40 117 2.67 ; 242 0.517 194 
53 0.258 120 0.10 109 2.61 466 246 0.561 170 
54 0.25 T 120 0.10 113 2.54 220 0.572 163 
53 0.258 120 0.40 109 2.61 470 256 0.583 169 
54 0.25 T 120 0.40 111 2.57 214 0.591 164 
53 0.258 100 0.10 99 2.54 403 218 0.575 143 
54 0.26 T 100 0.10 97 2.56 185 0.581 143 
53 0.258 100 0.40 101 2.51 414 229 0.608 141 
54 0.25 T 100 0.40 102 2.50 185 0.607 142 
55 0.508 70 0.10 Ti 2.53 917 590 0.509 221 
56 0.50 T 68 0.10 85 2.39 327 > 0.516 210 
55 0.50S 70 0.40 78 2.52 921 595 0.511 221 
56 0.50 T 68 0.40 say 326 

55 0.5085 80 0.10 88 2.49 1,011 640 0.499 243 
56 0.50 T 80 0.10 92 2.43 371 0.499 242 
55 0.508 80 0.40 94 2.42 1,025 652 0.500 238 
56 0.50 T 80 0.40 : ‘ 373 

55 0.50S 120 0.10 112 2.60 1,380 858 0.482 342 
56 0.50 T 120 0.10 104 2.69 522 ) 0.483 354 
55 0.508 120 0.40 111 2.61 1,383 858 0.482 344 
56 0.50 T 120 0.40 al , 525 , 

55 0.508 21.5 (L.F.) 0.10 133 2.2 1,764 1,082 0.465 461 
56 0.50 T 21.5 (L.F.) 0.10 131 2.74 682 0.465 465 
55 0.508 21.5 (L.F.) 0.40 131 2.74 1,761 1,074 0.462 466 
56 0.50 T 21.5 (L.F.) 0.40 687 : 
55 0.508 140 0.10 113* 2.91" 1,554 962 0.494 387 
56 0.50 T 140 0.10 ‘ a a 592 : 7 

55 0.50S 140 0.40 Li7* 2.67* 1,559 963 0.488 387 
56 0.50 T 596 


140 0.40 





* Estimated from Run 54. 


lowered from the knife edge during a reading so that 
possible movement of the arm on the agate surface 
could be corrected. 

In order to obtain a tare correction for the drag force 
on the support rod, every run was repeated with the 
arrangement shown in Fig. 5. The sphere was sup- 
ported by a rod extending downstream along the nozzle 


axis and connected to a stationary support arm. A 
dummy support rod, of the same exposed length and 
diameter as the original sphere support rod, was at- 
tached to the drag balance and placed in the same posi- 
tion as the original support rod. The clearance between 
the top of the dummy support rod and the bottom of 


the sphere surface varied from approximately 0.010 to 
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0.030 in. in the worst case, as compared to an exposed 
jength of the support rod in the nozzle ranging from 
9.96 to 2.71 in. The tunnel was then operated with 
the same flow conditions used for the original sphere 
plus support arm drag measurements, and the dummy 
support arm drag was measured with the balance under 
these conditions. Simple subtraction then gave the 
net drag force on the sphere alone, which was subse- 
quently used in reducing the data. Reference to 
Table 1 shows that the tare correction—1.e., the drag 
on the support rod alone—varied from 13 to 330 per 
cent of the sphere drag force (Runs 42, 41 and 50, 


51, respectively). 


Balance Calibration Procedure 


The balance was calibrated by adding weights to the 
left pan of the balance arm (see Fig. 7) and then ex- 
tending the quartz spring to bring the unit back to a 
null position. The counter reading indicating the 
spring extension and the corresponding calibration 
weight were then recorded. Three different spring 
sizes were used to maintain the desired accuracy over 
the whole range of drag forces. The calibration pro- 
cedure was checked generally before and after each 
run and, in every case, was also checked for ascending 
and descending values of the calibration weights in 
order to detect any hysteresis or creep effects. 


REDUCTION OF EXPERIMENTAL DATA 


The experimental results for each run are tabulated in 
Table 1. Pressure readings with the two probes were 
obtained at a number of points in the region to be occu- 
pied by each sphere in order to detect any variations in 
stream conditions over the volume occupied by the 
sphere and to allow determination of suitable average 
values of the flow properties for purposes of data re- 
duction. For Probe No. 15, it was assumed that the 
analysis summarized in reference 17, which was derived 
on the basis of a compressible nonviscous gas, was ap- 
plicable. This reference gives, in tabulated form, the 
relationship between the cone half-angle (5° in the case 
of Probe No. 15), the undisturbed stream static pres- 
sure p,, the undisturbed Mach Number , and the 
measured pressure at the surface of the cone, po. By 
a similar compressible nonviscous flow analysis, the 
pressure obtained with Probe No. 14, /i, is related to 
the stream static pressure p, and the Mach Number 
M by the Rayleigh formula.'* With these assump- 
tions, the two unknowns M and /p, could be obtained 
from the measured values of p2 and pi. 


Additional quantities recorded included the room 
temperature, the stagnation temperature and pressure 
in the settling chamber of the wind tunnel, and the 
arbitrary scale indication of a flow-rate meter. In the 
case of the drag balance readings, the centerline of the 
sphere and support rod always coincided with the 
centerline of the nozzle exit plane and was located 
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various distances downstream (along the nozzle axis) 
from this plane. The values of M and p, recorded in 
Table 1 correspond to the indicated axial positions. 

In order to put drag data on a nondimensional basis, 
it is customary to write the drag force as 


drag = Cp-(1/2) (pV? A) (4) 


where Cp is a dimensionless coefficient and the other 
quantities are in consistent units. In the case of a 
spherical shape, this type of equation can be further 
reduced to 

D = Kp-p V*d? (5) 
where Kp is then defined by 

Kp = (2/8) Cp (6) 


The quantity Ay» will be referred to as the drag co- 
efficient, and the experimental data were reduced to 
give the values of this drag coefficient in the present 
case. For convenience in the matter of units, the 
above equations can be transformed, by using the per- 
fect gas law and the definitions of the quantities in- 
volved, into 


D = 1.23 X 10-2 Kp: p,M°d? (7) 


where the units are given in the nomenclature list. 
Calculation of the quantity A, from the experimental 
data is summarized in Table 1. 

As a final step, the Reynolds Number corresponding 
to the conditions of the experiments was computed. 
The magnitude of the viscosity coefficient 4 was deter- 
mined from Sutherland's equation” 


‘ (F: + c) (7)’ 2 (8) 
KM = Ko 2. 2 T, ‘ 


where po = 10.0 K 10~* Ib. per sec. ft. at 450°F. abs.” 
and C = 202°F. abs." 

It was assumed that the flow in the wind-tunnel 
nozzle from the settling chamber to the test section, 
where the drag forces were measured, was adiabatic. 
If it is further assumed that the perfect gas law is ap- 
plicable, then the free-stream temperature can be de- 
termined in terms of the measured settling chamber or 
stagnation temperature 7» by the relation”! 


T/To = 1+ [(y — 1)/2] I? (9) 
The numerator of the Reynolds Number Vdp can be 


expressed in terms of M and p, by the following manipu- 
lation, using the perfect gas law and the definition of M: 


p = ps/ReT 
M? = (V/a)? = V?/yReT 
R = Vy/ReT (Mpd/p) (10) 


Finally, the Reynolds Number (for air) can be ex- 
pressed in a convenient system of units as 


R = 2.14 X 10-4 (Mp,d/uvV/T) (11) 
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The free-stream Reynolds Numbers calculated by this 
equation for each experimental run are tabulated in 
Table 1. 


SOURCES OF EXPERIMENTAL ERROR 


Drag Balance 


Forces in still air were measurable within +2 per cent 
maximum deviation, with the mean deviation less than 
+1 percent. These figures include creep or hysteresis 
effects, since calibrations were taken both before and 
after a run. During a run, the drag balance counter 
was read at least twice for each condition, the unit being 
balanced once, then moved well out of balance, then 
balanced again. The readings were consistent within 
+1 per cent. For Runs 50 and 51, where the tare 
correction was approximately 330 per cent of the sphere 
drag force, the mean deviation in the balance calibration 
data was less than +0.5 milligram. This would result 
in a maximum error of approximately +3 per cent 
in the measured sphere drag force for these runs. The 
error in measuring the drag force was less than this 
figure in all other runs. 


Pressure Readings 


The maximum shift in calibration constant for the 
oil manometer corresponded to +5 per cent variation 
about a mean value, and this figure may be taken as 
representing the accuracy of the impact-tube read- 
ings. The pressures associated with the cone probe 
were, in some cases, as low as 40 microns Hg, and in this 
range a variation in “‘zero’’ settings was the most serious 


source of error. Since the individual errors are difficult 
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to separate, Run 47 was carried out with pressures 
indicated by both the manometer and the Pirani gage. 
This run covered flow rator readings at 18, 30, and 80, 
where the manometer zero shift errors should be most 


serious. Ap and R were computed using both sets of 
pressure readings. In the worst case (Run 47, FR = 
80, d = 0.25 in., axial position = +0.10), the ma- 


nometer and Pirani gave, respectively, R = 126 and 116, 
Kp = 0.547 and 0.564 
cent in R and 3 per cent in Kp. 


or variations of about 8 per 


Tare Procedure 


Obtaining a tare correction for the drag of the sphere 
support required resetting tunnel conditions to corre- 
spond to those existing when the sphere plus support 
rod readings were obtained. This was accomplished 
by setting the flow rator to the same scale reading and 
by independently reading the indicated stagnation 
pressure. The least count of the flow rator scale was 
1 unit (in 18, 30, SO, 140 
perature variation over all the runs was under 10°F., 


see Fig. 8). The room tem- 
or less than 2 per cent change in temperature on the 
absolute scale. 

The small variations in adjusting the flowmeter, 
room temperature changes, and possible small distor- 
tions of the plastic nozzle would cause changes in the 
test-section flow conditions. One method of indicat- 
ing the cumulative effect of these errors, plus the pres- 
sure gage calibration errors, was to determine the Rey- 
nolds Number (per 1-in. characteristic dimension) vari- 
ation obtained at each flow rator setting. This could 
be done readily from the tabulated values of Table |, 
and the results are plotted on Fig. 8. This graph shows 
that the maximum variation in R at a given flow rator 
setting was from +4 to +8 per cent about a mean 
value. The resulting uncertainty in Ap does not ex- 
ceed +5 per cent. 


Interpretation of Probe Readings 


In addition to experimental errors, there was a syste- 
matic error involved in interpreting probe pressure 
readings to obtain the Mach Number and static pres- 
sure of the gas stream. As described previously, the 
probe readings were analyzed using theory applicable 
to a compressible, nonviscous gas, with standard bound- 
ary conditions for a continuum fluid. This procedure 
may lead to errors in the present case because of vis- 
cous and molecular flow effects. 


TABLE 2 


M M, R R, 
2.8 0.50 1,000 475 
2.8 0.50 100 48 
2.8 0.50 40 19 
2.1 0.56 1,000 605 
2.1 0.56 100 60 
2.1 


0.56 40 24 


(pi/ps) Eq. (3) 


(Reynolds Numbers R and R; Based on 0.3-in. Diameter Dimension of Probe No. 14) 


Per Cent 


(pi/ps) Eq. (13) Difference 


9.86 9.88 0.2 
9.86 10.02 1.6 
9 86 10.25 4.0 
6.17 6.18 0.2 
6.17 6.27 1.5 

é 


6.17 6.40 3 
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The viscous effects have been analyzed for an impact 


tube in subsonic flow and take the form,”* for a spheri- 


cal shape, 
V [3 — (83/110) M? 
b= pte 3 )M*] (12) 


d 1+ (0.457/~/R) 


p. is the impact pressure that would exist for a com- 
pressible nonviscous fluid. Eq. (12) shows that the 
measured impact pressure will be higher than the 
“ideal” value. For supersonic flow, reference 22 sug- 
gests use of Eq. (12), with p, being the stagnation 
pressure that would occur in a compressible nonviscous 
fluid after passage through a normal shock wave. The 
complete equation then becomes (for a spherical shape) 


[iy + 1)/2 Mey”O7-” 
[2y/(y+1)) M2=-[(y- 1D) /(v tt DJ? 
( 27 We Y- ye [3 — (83 110)Mi*}l 0,2 
vyt+1 y+i1/ Ri. 1+(0.457/~VR,) $ 


The actual impact tube (Probe No. 14) was source- 
shaped rather than spherical. However, Eq. (13) 
would be expected to give approximately correct re- 
sults for the actual probe. More significant than this 
possible discrepancy is the fact that Eq. (13) has never 
been checked experimentally, and the assumptions on 
which it is based may be in error. 

Only the first term of Eq. (13) has been used in re- 
ducing the experimental data. To estimate the mag- 
nitude and effect of the possible error in this procedure, 
the correction of Eq. (13) has been computed for the 
range of values of \/ and R covered in the tests and is 
summarized in Table 2. 

The following conclusions can be made regarding the 
probable effect of viscous corrections on interpretation 
of the impact tube pressure : 

(1) The correction is important only for the lowest 
values of R encountered in this investigation. For 
R = 40, the smallest value encountered in any run 
(Run 51), the correction on p;/ps; was less than 4 per 
cent. The corresponding error in M would be ap- 
proximately 2 per cent; that in Ap, about 2 per cent; 
and that in R, about 4 percent. At R = 100, the indi- 
cated error is less than 1 per cent and decreases rapidly 
for higher R. 

(2) The correction always consists of an additional 
term. This results in a value of M which is too high if 
the correction is neglected by utilizing Eq. (3) in re- 
ducing the experimental data, as was done here. 
Carrying through a sample calculation shows that the 
use of a value of M which is too high gives a value of Kp 
Which is too low. This means that the rapid increase 
in Kp for decreasing R (see discussion of results) would 
be accentuated rather than reduced if the correction 
of Eq. (13) had been utilized in reducing the data. 

This discussion has been concerned entirely with 
viscous effects on an impact tube, and the problem of 
interpretation of the readings of Probe No. 15 re- 
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Fic. 9. Test arrangement to determine effect of cross-stream 
support. 
mains. Probe No. 15 is a 5° half-angle cone, for 


which the nonviscous theory gives ratios of probe pres- 
sure, ps, to static pressure, p,, ranging from 1.10 to 1.16, 
for M from 2.1 to 2.8. While viscous effects undoubt- 
edly occur, probably in the form of a boundary layer 
that changes the pressure distribution from the ‘‘ideal”’ 
values, they would not be expected to have a large 
effect; a 50 per cent change in the pressure increment on 
the cone surface would result in less than 8 per cent 
change in the static pressure deduced from the experi- 
mental probe pressure. No experimental or additional 
theoretical conclusions are available at the present time 
for the problem of viscous compressible flow over a cone. 

Finally, the effect of slip or other molecular flow phe- 
nomenon on either probe,has not been evaluated. The 
largest M;/R, value (based on a 0.3-in. characteristic 
dimension) was approximately 0.03. From the indi- 
cations of Millikan’s data (see Introduction), it would 
be expected that the molecular flow effects would be 
small, even for this condition. For most of the experi- 
mental range, the small values of 14,/R, would indicate 


a negligible correction. 


Method of Sphere Support 


Support of the sphere by a rod of finite size will 
disturb the flow and cause an error in the drag coeffi- 
cient which is not completely eliminated by the method 
of tare correction outlined earlier. Support from 
the rear, by a rod aligned in the flow direction, is the 
recommended procedure. The cross-stream support 
in these tests can cause serious discrepancies in the re- 
sults when the Reynolds Numbers are in a critical re- 
gion that is of the order of magnitude of R = 10°. 
The effect has been traced to a transition from laminar 
to turbulent boundary layer, which changes the posi- 
tion of the boundary-layer separation point and alters 
the pressure distribution on the sphere. In the present 
case, however, the Reynolds Numbers are low enough 
so that the boundary layer remains laminar,** and the 
marked effect apparent at higher Reynolds Numbers 
would not be expected. 
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Fic. 10. 


In order to determine experimentally the effect of a 
‘ cross-stream support for the present operating condi- 
tions, the 0.5-in. sphere and support rod were re- 
mounted on the drag balance as indicated on Fig. 9. 
A rod, equal in diameter to the original 0.060-in. sup- 
port rod, was arranged so that it could be moved either 
out of the nozzle flow or directly against the top side 
of the sphere (with approximately 0.010-in. clearance). 
Drag measurements were obtained when the rod simu- 


TABLE 3 
(Data from Runs 71 and 72) 





Flow 
Rator 
Scale 
Reading fi fo f 
68 1.020 0.424 1.04 
1.011 0.412 1.02 
100 1.017 0.42* 1.03 
120 1.034 0.417 1.06 
140 1.032 0.416 1.06 
1.021 1 


* Estimated. 


04 


Experimental drag coefficient Kp vs. Reynolds Number R. 


lating the cross-stream support was out of the stream 
and then were repeated with the rod moved into the 
position indicated on Fig. 9. The experimental values 
of f; (the ratio of total drag force without cross-stream 
rod to total drag force with cross-stream rod in posi- 
tion) are shown in the second column of Table 3. By 
comparing the total drag force obtained for the same 
flow conditions during the original runs (Table 1) 
with the force obtained for the arrangement of Fig. 9, 
the tare force for the downstream sphere support of 
Fig. 9 was determined. The ratio of tare force to total 
force with cross-stream rod in position is entered in the 
third column of Table 3 (as fo). The net effect of the 
cross-stream rod on the sphere drag force then was com- 
puted and is entered in the fourth column of Table 3. 
The quantity f is the ratio of sphere drag force without 
a cross-stream rod to that with a cross-stream rod if 
position. 

The effect of the cross-stream support was to give 
values of the sphere drag coefficients /ower than those 
determined with a downstream support or sting. As- 
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suming that the sting support gives more nearly correct 
yalues of the drag coefficient, the drag coefficients in 
Table 1 should be increased by approximately 4 per cent. 
The experimental correction ranged from 2 to 6 per 
cent, based on the results with a 0.5-in. diameter sphere 
tabulated in Table 3. 

The 0.5-in. diameter sphere and 0.060-in. diameter 
support rod represented the largest ratio of support rod 
diameter to sphere diameter used in this investigation 
(Fig. 6). Evidence that the experience with the 0.5- 
in. sphere is representative of the behavior to be ex- 
pected from the other test configurations is gained from 
Fig. 10, which shows that drag coefficients obtained for 
spheres of different sizes at the same Reynolds Numbers 
gave the same measured drag coefficients. This was 
verified at three different Reynolds Numbers, where the 
(0.1- and 0.25-in., the 0.25- and 0.50-in., and the 0.50- and 


|-in. sphere diameter data overlapped. 


DISCUSSION OF EXPERIMENTAL RESULTS 


The values of the drag coefficient Kp are plotted as a 
function of the Reynolds Number R on Fig. 10, and 
the experimental points are seen to define a single 
curve with a scatter that is the same order of magni- 
tude as the experimental errors. Reference to Table 
1 shows that the Mach Number range was small in 
these tests, (J varying from 2.1 to 2.8. This change 
is not sufficient to permit conclusions regarding the 
Mach Number effect; the change in Kp at normal pres- 
sures over the same Mach Number range was less than 
3percent.* The striking fact is the significant variation 
in drag coefficient with Reynolds Numbers less than 
approximately 10%. 

The cumulative effect of the various sources of error 
discussed in the section on ‘‘Sources of Experimental 
Error’ can be estimated from Fig. 10. The standard 
deviation in Ay for all the data is approximately +5 
per cent. The largest variation from a faired curve 
through the experimental points occurred at R ~ 140 
and R ~ 165, for Runs 53 and 54 (0.25-in. diameter 
sphere) at flow rator settings equal to 100 and 200, 
respectively. These values are approximately 9 per 
cent higher than the faired curve. The reason for the 
large scatter at these pointsis not clear. Data obtained 
during the same runs, at R = 200, check well with the 
points obtained at the same R with the 0.50-in. diam- 
eter sphere, so that some gross error in calibration or in 
the drag balance adjustment is unlikely. 

In analyzing the drag forces obtained with aerody- 
namic shapes, experience has shown that the total 
drag force can be separated into two components. The 
first component is called wave drag and is the net drag 
force that results from integration of the pressure dis- 
tribution over a body, the pressures being measured 
normal to the body surface. The second component 
is a skin friction drag, which is the integrated effect of 
the tangential shearing stress components over the sur- 


face. The latter component often can be estimated 
from boundary-layer theory. In general, this estimate 
requires knowledge of the type of boundary layer 
whether laminar or turbulent—and the location of 
points where the boundary layer may separate from the 
body. 


A common procedure is to obtain experimentally 
a pressure distribution for a particular aerodynamic 
shape under given flow conditions, integrate the result- 
ing distribution to obtain the wave drag, and then add 
skin-friction drag computed from boundary-layer 
theory. This procedure, apparently, is satisfactory 
for subsonic speeds with streamlined bodies at high 
Reynolds Numbers, where the skin-friction com- 
ponent is a small part of the total drag and where sepa- 
ration of the boundary layer is present for only a small 
percentage of the body surface. In the case of a sphere, 
the problem is complicated by the fact that separation 
inevitably occurs (unless the Reynolds Numbér is so 
low that the flow is in the Stokes region), and the wave 
drag is greatly affected by the behavior of the boundary 
layer. However, this procedure—separating drag into 
two components—may be utilized to give a rough esti- 
mate of the nature of the aerodynamic phenomena 
underlying the results shown in Fig. 10. 

For the low Reynolds Numbers of the present in- 
vestigation, the boundary layer is assumed to remain 
laminar up to the point of separation. This assump- 
tion implies that the shock wave associated with the 
sphere at supersonic speeds is sufficiently far upstream 
so as not to affect the boundary layer. The flow be- 
hind the shock wave then may be characterized by a 
Reynolds Number R,, which is related to the free- 
stream R by the known process occurring in a normal 
shock wave. R,; was calculated, for known values of 
M and R, by determining the ratios Vi/V, p:/p, and 
T,/T from reference 25 and y,/u from Eq. (8). The 
R, derived in this way actually would be expected to 
give the correct value only on the stagnation line im- 
mediately behind the shock wave, since the shock 
has different curvatures at other points of the flow; 
but this single value is a plausible one for characterizing 
the subsonic region behind the shock wave, particularly 
for the limited Mach Number range covered in these 
experiments. 

With this hypothetical picture of the flow process, 
an explanation of the trend of the experimental values 
can be attempted. The drag coefficient is known at 
high Reynolds Numbers, where the boundary-layer 
thickness is negligibly small compared to the sphere 
radius. It will be assumed that the drag coefficient at 
low R can be determined by using the known coefficient 
at high R for a sphere that has an effective diameter 
(d + 26). 6 is some value, averaged over the sphere 
surface, of the boundary-layer displacement thickness. 
The product of the drag coefficient by d? [see Eq. (5) ] is 
then replaced by 
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computed behind a normal shock wave. 


(drag coefficient) X 
[1 + (26/d)]* d? 


(drag coefficient) (d + 26)? 


so that the effective drag coefficient Kp is 
Kp = Kp’ [1 + (26/d)]? 
From reference 3, 


R= 10° 


where K,’ is evaluated for 6/d = 0. 


Kp’ ~ 0.38, 2.1 < M < 28, 


Separating Ky’ into wave-drag and skin-friction com- 


ponents and assuming that the skin-friction component 
varies like 


a/VR, a constant 


then Kp can be written 


Kp = [0.38 + (a/VR,)] [1 + 2(6/d) |? 
Finally, from the results of reference 24, 26/d is taken 
in the form 


constant 


2(8/d) = B/V/R, B 
and 
Kp = [0.38 + (a/VRi)] [L + (8/VR)]? (14) 


The two assumed constants, a and 8, can be evalu- 
ated from the experimental data. Taking the values, 
determined by successive approximations, a = 0.52 
and 6 = 1.0 gave a curve that fitted the experimental 
data as indicated on Fig. 11. 

The equation naturally checks the known experi- 
mental value of Kp for large R; (~10°), because it was 
derived as an extrapolation from this value. More- 
over, the assigned numerical values of the constants 
a and 8 are the order of magnitude that would be ex- 
pected from boundary-layer theory. For example, the 
value of 8 is given in reference 24* as ranging from 0.45 
at the stagnation line to 3.2 at a point about 81° down- 
stream (measured between the flow direction and a line 
passing through the sphere center and the separation 
point). The value of 6 is nearly constant for the first 


* Actually only the boundary-layer thickness is given explicitly 
in this reference. The displacement thickness was calculated 
from its definition in terms of the boundary-layer thickness. 
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45° of the sphere and then increases rapidly, so tha 
8 = 1.0 seems a reasonable value to characterize ap 
average displacement thickness. The 
reference 24 also provides enough information to per. 
mit calculation of the approximate magnitude of ¢ 


analysis of 


This factor varies from zero at the stagnation point to 
approximately 1.8 at the 45° point, with a zero pressure 
gradient (and therefore decreasing local skin friction 
where 
While it is concluded from the 
boundary-layer analysis of reference 24 that the value 


at 75°, and, finally, a zero value again at 81 
separation occurs. 


a@ = 0.52 is correct in order of magnitude, a slightly 
higher value (say, approximately 2.0) would be ex. 
pected. This suggests possible molecular flow effects, 
with “‘slip’’ reducing the skin-friction coefficient. 

Further evidence of molecular flow effects is ob. 
tained by comparing these data with drag coefficients 
obtained in subsonic flow in liquids for the same Rey. 
nolds Number range. It was argued previously that 
the total drag can be separated into wave-drag and 
skin-friction components. At high R, where the skin- 
friction effects are small, the wave-drag components are 
different for supersonic versus subsonic Mach Numbers. 
If the energy loss in the shock wave at supersonic 
speeds remains essentially constant as R decreases (M 
being approximately constant), then one would expect 
the supersonic and subsonic drag coefficients to follow 
the same trend when compared on the basis of the same 
Reynolds Number (R; in the supersonic case), pro- 
vided there were no molecular flow effects. The data 
would be displaced because of the different wave drags 
in the two cases but should define essentially parallel 
curves. This statement neglects compressibility ef- 
fects, but these are small for cases where WM, < 0.5. 
The subsonic, incompressible drag data summarized in 
reference | have been plotted on Fig. 11. Examination 
of this figure shows that the present data do not in- 
crease so fast as the subsonic curve and, in fact, cross 
this curve. A possible explanation is that “‘slip” ef- 
fects are decreasing the skin-friction drag in air. 

If molecular flow effects are present, it follows that 
Eq. (14) should break down as R decreases below the 
limit of the present experimental data, since the equa- 
tion contains no function of M/R. Moreover, Eq. 
(14) gives a drag coefficient increasing without limit as 
R tends to zero, and it is known from kinetic theory 
that the drag coefficient becomes independent of R in 
the large mean free path region (M/R > 10). The 
values of Kp computed for this region" are included on 
Fig. 10 for 17 = 2 and 3 (these values are based on 
assumed equal temperatures of impinging and re- 
flected molecules). 

The assumption that the shock wave does not affect 
the boundary layer is particularly suspected. Fig. 12 
shows a typical photograph of the shock wave for 4 
0.7-in. sphere obtained by a glow method.”* 7” Similar 
photographs obtained with air and a l-in. diameter 
sphere indicate that the beginning of the shock wave is 
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approximately 0.1 in. in front of the sphere under con- 
ditions where the molecular mean free path is approxi- 
mately 0.006 in. in the undisturbed flow (Run 43, 
M = 2.6, R = 700). The thickness of a laminar 
houndary layer at the stagnation line (computed from 
the results of reference 24) would be 0.07 in. If the 
shock-wave thickness were of the order of magnitude 
of 10 molecular mean free paths, or 0.06 in., an inter- 
action between the shock zone and the boundary layer 
might occur. 

It will be apparent from the preceding discussion 
that the available data do not permit separation of 
viscous, molecular flow, and shock-wave effects. The 
simplified assumptions made in deducing the form of 
Eq. (14) imply that only changes in viscous effects are 
occurring, and this explanation seems to fit the data 
satisfactorily, at least for R > 40. According to the 
estimates of Tsien,'® naolecular flow effects should be 
present in this region (since M/+1/R = 0.4). A com- 
plete understanding of the relative importance of 
these various factors requires additional experimental 
work. 

Measurement of drag forces in air at low subsonic 
speeds (M < 0.3) and with R varied over a wide range 
from a maximum of approximately 400 would permit 
direct comparison with incompressible flow data. Since 
there would be no shock wave, deviations from the 
incompressible flow results would indicate the magni- 
tude of molecular flow effects. By comparing the re- 
sults of Fig. 11 with these new subsonic data, conclusions 
regarding interaction between the shock wave and 
boundary layer might be possible. Measured pressure 
distributions over the sphere surface would be valuable, 
but the experimental difficulties are formidable for 


small models. 


CONCLUSIONS 


(1) Drag coefficients have been measured for 
spheres in supersonic air streams at low Reynolds 
Numbers. These previously unavailable data cover 
a Reynolds Number range (based on the sphere diam- 
eter) from approximately S800 to 15, at Mach Num- 
bers from 2.1 to 2.8. The data show a marked increase 
(by a factor of approximately 2'/2 over the test range) of 
drag coefficient with decreasing Reynolds Number, 
whereas available information at supersonic speeds, for 
R = 10°, could be correlated as a function of Mach 
Number alone and was independent of Reynolds 
Number variations. 

(2) The precision of the experimental values of the 
drag coefficient Kp is estimated to be +5 per cent. 
This is the standard deviation of all experimental 
points from an average curve drawn through them. 
The use of a cross-stream support has been shown 
experimentally to give drag coefficients that are ap- 
proximately 4 per cent lower than the correct values. 
The effect of systematic errors introduced in the inter- 
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pretation of probe pressure readings cannot be com 
pletely evaluated at this time. 


(3) The data have been correlated, within the limits 
of the experimental errors, over the entire R range 
covered in these tests by an equation based on the con- 
cept of separation of the total drag into wave-drag and 
skin-friction components, with the form of the skin- 
friction component corresponding to an assumed lami- 
nar boundary layer up to the separation point. This 
formulation resulted in a function with two constants, 
which were determined to give the best,fit with the 
experimental data. Comparison with an available 
boundary-layer analysis for a sphere indicated that the 
values of the constants are the same order of magnitude 
as indicated by the theory. The equation fails for the 
limiting case of free molecular flow. The limited 
Mach Number range in the present investigation did 
not permit conclusions regarding the effect of variable 


M. 


(4) A photograph of the shock wave associated with 
a sphere indicates that the shock zone may interact with 
the boundary layer in the vicinity of the stagnation 
line of the sphere. If such an interaction occurred, 
the assumptions made in deducing the form of the pro- 
posed correlation equation would be incorrect. 


(5) Comparison of the data with available drag co- 
efficients for subsonic incompressible flow indicates 
that molecular flow effects “‘slip’’ may have been present 
for the lower part of the Reynolds Number range. 
Additional experimental data are required before the 
relative importance of viscous, molecular flow, and 
shock-wave effects can be evaluated completely. 
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| The Secondary Flow in a Cascade of Airfoils 


in a Nonuniform Stream 


H. B. SQUIREt AND K. G. WINTER 
Royal Aircraft Establishment, Farnborough, England 


SUMMARY 


A method is given of calculating the secondary flow in a cascade 
of airfoils in a stream with a variation of approach velocity along 
the span of the airfoils. The results are compared with measure- 
ments of the secondary flow at a cascade at a corner in a return- 
circuit wind tunnel, where the velocity gradient is obtained as a 
normal consequence of flow through a diffuser. The agreement 
between theory and experiment is satisfactory and indicates that 


the theory is sound. 3 


List OF SYMBOLS 


= Cartesian coordinates in plane normal to span of 


x, 9 
cascade 

a, 8 = orthogonal curvilinear coordinates in (x, y) piane 

3 = distance measured along span from one end wall 

h = the following: elements of length in (a, 8) coordinates 


are hia and hig 
r = radius at any point of passage between two vanes 
mean radius of passage 


- 
ll 


2b = width of passage at exit 

21 = height of passage at exit 

6 = angle of deflection at corner up to point a 

€ = total angle of deflection at corner 

v = velocity vector with components u, v, w in directions 
a, B, 2 (=y) 

U = velocity of approach to corner 

0), @, = values of v, w downstream of corner 

w = vorticity vector with components &, 7, ¢ in directions 
a, B,2(=y) 

¥ = stream function such that »; = dy/0z, wi = —Ody/dp 

p = static pressure 

p = density 


(1) INTRODUCTION 


les EXISTENCE OF SECONDARY FLOW downstream 
of a cascade of airfoils, in a stream of nonuniform 
velocity, is well known. It consists usually of a pair of 
vortices in each passage through the cascade with the 
character shown in Fig. 1. The extent of the second- 
ary flow depends on the degree of nonuniformity of the 
approaching stream, and two types of secondary flow 
can be distinguished. In the first type, the stream is 
uniform except for a thin boundary layer at the walls. 
The secondary flow is then limited to layers near the 
end walls and in the wakes of the cascades. This type 
of motion has been investigated by Carter and Cohen.! 
In the second type of flow, the stream is nonuniform 
over a large part of the span of the vanes, as will occur 
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if the cascade is placed at the end of a diffuser in a wind 
tunnel. The vortices shown in Fig. 1 thus extend over 
a large part of the cross section and eventually combine 
to form a pair of large vortices far downstream. 

This note deals only with the second type of flow, 
in which the region of nonuniformity is large com- 
pared with the gap between the vanes. In this case, 
although the nonuniformity of the stream may have 
been created by frictional effects, the development of 
the secondary flow system can be considered as being 
independent of friction and can be calculated from the 
equations of motion of an ideal fluid. 

There is an analogy between the secondary flow and 
the motion of a gyroscope. If a disc, which is rotating 
about its axis of symmetry, is turned about an axis in 
its plane, it develops an angular velocity about a third 
axis perpendicular to the first two. The analogy with 
the present problem is obtained by supposing the disc 
to lie originally in the plane containing the stream direc- 
tion and the spanwise direction of the cascades. Its 
original rotation is associated with the velocity gradient 
along the spanwise direction of the cascades. The disc 
and the stream are then turned, on passage through 
the cascade, about an axis along the cascade span. The 
rotation of the disc about its third axis corresponds to 
the secondary flow of the character indicated in Fig. 1. 
THEORY 


(2) 


(2.1) General 
For steady flow of an inviscid incompressible fluid, 
we have the equation® 


0 (1) 


curl (v X w) 
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FIG.1. CHARACTER OF SECONDARY FLOW AT 
OUTLET OF A CASCADE IN A DUCT. 
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FIG.2. ORTHOGONAL 
THROUGH CASCADE. 


where v is the velocity and w is the vorticity. With 
orthogonal curvilinear coordinates (a, 8, y), for which 
the elements of arc are h', 5a, h258, h3éy, the component 
velocities are (u,v, w) and the component vorticities are 
(é, n, ¢); the first component of Eq. (1) is 


(2) 


a ra) 
- [h3(un — vt)] = . [ho(wt — uf)] 


Also 


div v = 0, div w = 0 


which become, in (a, 8, y) coordinates, 
ra) eo. re) ; 
a (heh3u) + ry, (h3hyw) + oy (hyhow) = 0 (3) 


2 p D 


sea (hohzé) + a8 (h3hyn) + > (Ayhot) = 0 (4) 


a "y 


From Eq. (2), making use of Eqs. (3) and (4), we ob- 


tain 
(* o,2a9 ° (=) ‘ 
hyOa hOB hg Oy/\hy ns 
t 
(; 24324! 2 )(*) (5) 
hy Oa hy faye} hs Ov hy 


This equation and two similar ones are Helmholtz’s 
equations in curvilinear coordinates. If the motion is 
not steady, we have to add a term (0/0f) (€/h;) to the 
left-hand side of Eq. (5). 





NETWORK FOR’ FLOW 


We now suppose that the curvilinear coordinates are 
a(x, y), B(x, y), and z, where x and y are Cartesian co- 
ordinates in the plane normal to the span of the cascade 
and z is measured along the span. a is the velocity po- 
tential and 8 is the stream function of the two-dimen- 
sional motion that would be present if the approaching 
stream were uniform (Fig. 2). The flow in one of the 
passages between the airfoils and its extensions ahead 
of and behind the airfoils lies between two stream lines 
8 = constant, and we shall restrict consideration to a 
single passage bounded in this way. Then, for these 
particular coordinates, a, 8, and z (= y), we may put 


h, = he = h, he = | (6 





and the equation of continuity [Eq. (3) ] becomes 


= (hu) + : (hv) + - = 0 (7) 
Oa op. * Os 

The quantity 4 can be regarded as proportional to 
the reciprocal of the velocity of the two-dimensional 
flow at any point in the field. We choose the scale of 
a and £ so that h is equal to unity far ahead of the cas- 
cade; h will also be equal to unity far behind a cas- 
cade that gives a deflection without expansion or con- 
traction. 

The vorticity components are? 


] a) Oo low Ov 
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SECONDARY FLOW 


1 fo re) | ou 10w 
= hw) — — (hgw) | = — -=— 9) 
1 Ashi i ~— Oa (haw) oz hoa ( 


3) re) 
t= . | (hov) — (in) | = 
* — Inh, LOa 0B 


1fo re) 
— (hv) — — (hw 
h? E (ho) reye} (iu) | 


We now consider the flow past a cascade for which 
the approach velocity distribution is nonuniform. 
Then, if the degree of nonuniformity is not large, the 
two-dimensional flow referred to above can be regarded 
as a first approximation to the actual flow. We assume 
that this is so and shall calculate the actual flow by 
treating it as a linear perturbation from the two-dimen- 
sional flow. 

If the approaching stream has velocity U (where // is, 
in general, a function of 6 and z, such that the velocity 
variation is not large), then the vorticity components 
n and ¢ are small quantities of the first order and & is 
zero in the approaching stream. It will be assumed 
that, on passing through the cascade, the vorticity 
component £ becomes a small quantity of the first order. 


For the two-dimensional basic motion, we should 
have 


nm = 0, ,=wo = (10) 


and, since we are only considering small perturbations 

on this basic motion, we can regard v and w as small 

quantities of the first order; also (0/Oa) (hu) is a small 

quantity of the first order, since it is zero for the basic 

motion [Eq. (7)]. We shall neglect all quantities of 

the second order in Eq. (5), which becomes, on making 
EO 


uo (‘) (*) ‘ n O (‘) 
hda\h/  hda\h hop \h 
ra) () _ mo (*) 
da lu) u2d08\A 


The second Helmholtz equation gives, to the same 
order of approximation, 


use of Eq. (7), 


or 


(11) 


(0/d0a) (n/h) = 0 
Hence, from Eq. (9), 


n/h = 0U/odz (12) 


which is its value in the approaching stream. Substi- 
tuting from Eqs. (10) and (12) in Eq. (11), we get 


— 2h? (log h) 
ds 0B” 
The neglect of the term depending on 0U/O6 is in ac- 
cordance with the approximations introduced above. 
We now integrate Eq. (11) with respect to a from far 
upstream to far downstream, obtaining 
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FIG.3. SKETCH ILLUSTRATING CHANNEL USED 
TO REPRESENT ONE PASSAGE OF CASCADE. 
(hé) gu f- ne (1 h) d (13) 

i), = —2— h? og h) de : 
1 a dB Pe ada 


where the suffix 1 denotes the value far downstream. 
This formula enables us to calculate the component 
vorticity along the stream direction downstream of the 
cascade. 

Since the integration in Eq. (13) is along the stream 
lines of the basic flow, we can use a crude representation 
of the basic flow for the calculation of the integral. For 
this purpose the motion is assumed to be uniform in 
front of and behind the cascade, with h = | there; 
within the passage we assume that the basic motion is 
the irrotational flow in a circle about a vortex at the 
For this latter motion, 


k= 


center of the circle. 


a = 760+ const., B = 7% log (r/r), r/o 


where 7, 79, and @ are defined as in Fig. 3. 


(0/08) (log h) = 1/1 


Hence, 


and 


da = ry d0 
in the passage, and thus 
S h2(0/08) (log h) da = (r/ro)?0 


As a further approximation we shall assume that the 
width of the passage is small compared with the radius, 
so that r does not differ greatly from 7» across the 
passage. We may then put r/ro = 1 in the above ex- 


pression, which becomes 


JS h(0/0p) (log h) da = 6 (14) 


For the complete basic motion represented approxi- 
mately by uniform approach flow, curved flow in the 
passage bounded by two circular ares, and uniform flow 
downstream, we find that the value of ¢ in the outflow 
is, from Eqs. (13) and (14), 
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FIG.4. VELOCITY DISTRIBUTION, APPROACHING 
FIRST CORNER OF 4'x3' WINDTUNNEL. 


r. Ow) Ov “ OU 


—%e — (15) 
op Oz Oz 
where ¢ (Fig. 3) is the total angle of deflection of the 
stream; we have put #4 = | downstream of the bend, 
and 2, w, are the 8, s components of velocity far down- 
stream. Thus, the axial vorticity in the stream leaving 
the cascade 1s equal to the velocity gradient in the ap- 
proaching stream multiplied by twice the angle of deflec- 
tion. 

If the motion settles down to a steady secondary 
flow far downstream superposed on the general stream, 
we shall have 0u/Oa = 0 in Eq. (7), so that the equa- 
tion of continuity becomes 


(07/08) + (Ow,/d0z) = O (16) 


£1 


Introducing a stream function of the secondary motion 
defined by 7; = Oy/0z and w, = —Oy/08, Eq. (16) is 
satisfied and Eq. (15) becomes 


(O°y/0B2) + (0% /dz?) = 2€(0U/dz) (17) 


It is assumed that there is no flow across the extensions 


of the cascades denoted by 6 = +6 (Fig. 3) or at the 
end walls z = 0, z = 2/, so that the boundary condi- 
tions may be taken to be y = O for 8 = +6 and for 
$= 0.5: = 27. 


The theory given above is a development of the theory 
of von Karman and Tsien for small deflections;’ there 
are differences between the two theories with regard to 
the retention of some second-order terms. 


(2.2) Solution for a Bend of High Aspect Ratio 


If the width of the passage is small compared with 
the height, the secondary flow will be mainly parallel 
to the span of the airfoils except near the end walls. 
The solution of Eq. (17) corresponding to this type of 


flow is 
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yy = (8? — b*)e(OU/Oz) | 
4h = (B? — 5b?) €(0°?U/dz") (18 
w, = —2Be (OU /0z) 


. 


provided that 6°(0°U/0z*) is negligible in comparison 
with 0U/0z; for the velocity distribution considered 
below, this condition is satisfied except near the sta. 
tions marked A in Fig. 4. 

The expression (18) may be regarded as a particulg 
integral of Eq. (17) and is the required solution away 
from the end walls and from A (Fig. 4). Near the end 
wall z = 0, we have to add a complementary function 
to this particular integral to satisfy the boundary con. 
dition there. This complementary function is the well. 
known ‘“‘surface-wave’’ type of solution of Laplace's 
equation, and the complete integral is 


( ° b2) ol’ 
v = B- mara 


in 


2 1) rz (2 l 
z B, exp | - al : | cos | et ef) (19 
n=0 


2b 2b 
This form satisfies the condition y = O for B = +b, 
and the coefficients B, are determined by the condition 
¥Y = Oforz = 0. Putting z = 0, Y = 0 in Eq. (19), we 
obtain 


oe r ou : ; (2n + 1)rB ; 
ee & )_.+ > By cos | ; 2b — 


n=0 


Hence, 


32 (—)"b'e | 
(2m + 1)* r® \Oz/z=0 


so that the coefficients decrease rapidly with increase 
of n. 

We may thus, to the order of accuracy with which we 
are concerned, neglect all except the first term in the 
series, in which case we may replace the form (19) by 


the expression 


v = (8? — Bb’) 7 ee ee (20 
a * Os " 2b : 


The values of the two component velocities are, then, 


ae 2 ou ; 4 
ri a (B* — O*)e oe ent — + 


2b 


c », oU 12 

(Bp? — be 5s? [ — exp (- =) (21 

oU ° - 

WwW, = —2Be [ — exp (- = )| (2? 
QZ 2b! 


Eqs. (20), (21), (22) are valid for the upper half of the 


channel. 


(3) EXPERIMENTS 


An experimental check of the theory was made 
the first bend of the R.A.E. 4- by 3-ft. Wind Tunnel. 
The velocity distribution approaching the core! 
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SECONDARY FLOW IN A 


(Fig. 4) is such as to give an almost linear variation of 
velocity for the first 15 in. away from the wall. The 
center passage of the cascade, which has a span of 6 ft. 
and a width 2b = 4.8 in., was used. The experiments 
consisted of measurements with a yawmeter of the 
angles of deflection of the stream in the 6 and gz direc- 
tions. The traverses covered the upper part of the 
span of the turning vanes over which the approach ve- 
locity variation was appreciable. Because of the un- 
certainty of the main stream velocity, the angles were 
obtained by reading the pressure difference between 
the two tubes of the yawmeter for several angular 
settings and interpolating for zero difference. The 
slopes of the curves then also give the local velocity. 

To avoid interference from neighboring passages, 
the trailing edge of the outer of the two turning vanes 
enclosing the center passage was extended as shown in 
Fig. 5. The heads of ‘the yawmeter were inserted 
about !/,in. into this passage. 


(4) COMPARISON OF EXPERIMENTS WITH ‘THEORY 


A direct comparison of the measured secondary ve- 
locities with the theoretically predicted values is given 
in Figs. 6 and 7. Fig. 6 shows the velocities in the 
spanwise direction; Fig. 7 shows those normal to the 
span. In both cases the velocities are made nondi- 
mensional by dividing by the mean speed Uy across 
the tunnel section approaching the corner. 

The theoretical values of 7; and w, are given by Eqs. 
(21) and (22), using the velocity distribution given in 
Fig. 4. 

It will be seen in Fig. 6 that agreement between 
theory and experiment is good. In Fig. 7, the agree- 
ment in the nonlinear region is not so good, though one 
might say that the experimental curve is a smoothed- 
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FIG.7. SECONDARY FLOW NORMAL TO SPAN 
EXPERIMENTS. 


THEORY 4& 


out version of the theoretical curve. This really 
means that second-order effects are coming into play. 
The exponential decay law predicted near the end 
wall is followed closely. 

In Fig. 8, a complete comparison of experiment with 
theory over the whole region considered is given by 
plotting the stream lines of the lateral motion—1.e., the 
lines Y = constant. The values of y appropriate to 
each stream line have been made nondimensional by 
dividing up Ub. The theoretical stream lines are 
given by Eq. (20) using the velocity distribution of Fig. 
4. 

The experimental stream lines have been obtained 
by integration of the secondary flow velocities. Ex- 
cept in the region near the end wall, the steps in the 
process are We choose the stream line 
along the boundaries of each passage as y = 0. Then, 
the value of y at any point (8, z) is given by 


- 
—b 
where the integration is performed at z = 


order that both vanes bounding a passage shall be 
stream lines, it is necessary that 


4 
/ w, dB = 0 
—b 


constant, and this condition was satisfied by 
The maxi- 


as follows: 


w, dB 


constant. In 


for z = 
adjusting the zero of the measurements. 


MEASURED FROM TUNNEL ROOF 
ZERO AT CENTRE OF GAP & POSITIVE TOWARDS 


— COMPARISON OF 


mum adjustment required was a change of the zero of 
Wy Un of 0.04. 
Near the end wall the stream lines were obtained 


from the integral f v, dz. To join the stream lines the 
0 


values of y given by the w; measurements were plotted 
against z and the value of z taken, for which y is a maxi- 
mum. Call these values 2, Yo. At 20 since Oy /dz = 0) 
the value of v; must be adjusted to zero and y will then 


v=Wt [vo dz 


The maximum adjustment required was a change 0 


be given by 


zero of v;/ Uy, of 0.04. 





The differences between theory and experiment are 


more apparent in Fig. 8, the maximum value of —y Uv! 
being 0.255 from the measurements and 0.204 from the 
calculations. However, the general characteristics 
of the flow are the same and the experiments confirm 
that the physical basis of the theory is sound. The dis 
crepancies can be attributed to the second-order effects, 
which have been neglected in the theory. For example, 
consider an angular deflection of the flow of 25° at the 
end of a passage. Assuming this deflection to grow 
linearly round the curved part of a passage, then the 
lateral deviation of a stream line through the passagt 


is about 8 in. This deviation is sufficient to modily 


the velocity gradient appreciably, and, in fact, the 
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SECONDARY FLOW 


measurements show that the velocity gradient along the 
centerline of the passage becomes reversed for the first 
5in. away from the wall. 

(5) INDUCED DRAG 


The power loss associated with the secondary flow 


p ~2l *b 
J / U(w? + 
2 Jo —b 


therefore define an induced drag pressure 
drop coefficient 


ft f U ( ay 4 (=¥) 13 d 
tb/ Jo busy Uy Uy Pe 


C,; has been calculated both by integration of the 


isgiven by 


w,") dB dz 


We may 


Ci = 


pi 


measured secondary flow and from the theoretical dis- 


tributions. 
The measured value is thus of the order of half the 


The respective values are 0.032 and 0.019. 


profile drag pressure drop coefficient and should be 
taken into account in estimating corner loss. 

It is possible that there is some recovery of the energy 
in the secondary flow when mixing between the flows 
from adjacent passages takes place, and in the case of a 
wind tunnel there may be some recovery at subsequent 
corners. However, it is unlikely that this recovery 
would exceed 50 per cent. 


(6) CONCLUSION 


The experiments show that the theory gives a good 
approximation to the secondary flow in a cascade with 
a 90° corner and a large variation of approach ve- 
locity. 
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1945 











Free and Forced Vibrations of Rotating Blades 


ROBERT PLUNKETT* 
The Rice Institute 


SUMMARY 


The matrix equations governing transverse vibration of a rotat- 
ing twisted cantilever beam, such as an airplane propeller or tur- 
bine blade, are developed. From these equations the variation 
in natural frequency with rotational speed may be found, and 
some simplified methods for finding these curves are shown. 
The same approach is used for the case of forced vibrations of a 
blade due to oscillating loading, and the case of vibration in an 
airplane propeller due to yaw is considered as an example. It is 
shown that this approach can be expected to give results within 
about 10 per cent of experimentally measured ones if a rather 
simplified determination of the aerodynamic forces is permissible, 
as it is in normal flight. 


INTRODUCTION 


fee VARIATION OF NATURAL FREQUENCY of lateral 
vibration of a propeller or turbine blade with rotat- 
ing speed is of evident interest in design. This rela- 
tionship is both complicated and extremely tedious to 
-alculate. The methods in present use are the Ray- 
leigh or energy method" ? and the Myklestad method.’ 
The first of these suffers from being only an approxima- 
tion to the frequency and from the necessity of guessing 
The second may well be compared with 
However, insofar as is known 


a mode shape. 
the proposed approach. 
to the author, this is the first attempt to analyze an 
arbitrary blade* with due regard to rotating speed, 
twist, and blade angle at the same time, although either 
of the two methods named above could be modified to 
take these factors into account. The principal draw- 
back of the Myklestad method is that it inherently 
loses significant figures—that is, it is necessary to carry 
more than six significant figures in the untwisted case 
and more than twelve in the twisted one® to get two- 


figure accuracy in the answer. 
The Southwell approach leads to the relationship 
w? = wo” + a? 


where w is the natural frequency of vibration, wo is the 

natural frequency at standstill, and Q is the rotational 

speed in rad. per sec.; @ is sometimes taken to be a con- 

stant® 7 but should be considered as a function of Q. 

Fortunately, this turns out to be practically the same 
Received April 29, 1950. Revised and received July 28, 1950. 
* Assistant Professor of Mechanical Engineering. 


u + wd 
QO m, 0 


Mm, — My 


relationship for all similar blades, so that once the curves 
are drawn for one blade the results may be easily applied 
to all blades of that family. It should be recognized, 
however, that a, the Southwell factor, may vary within 
wide limits even under present-day operating condi- 
tions, and, as blades become more limber and are ro- 
tated at faster speeds, this will become of more impor- 


tance. 


EQUATIONS 


In a recent paper,’ the author has derived a set of 
matrix equations that closely describe the deflection of 
a rotating blade caused by a known transverse force. 
By a well-known application of D’Alembert’s principle, 
these may be used to find the equations governing the 
behavior of a rotating blade vibrating transversely. 
Since the vibration, to the first order, is sinusoidal, the 
maximum deflection, ignoring damping, is one that 
would be caused by transverse forces equal to the mass 
times the square of the vibration frequency times the 
local amplitude of vibration.*® 

In the ensuing paper, the same notation as that of the 
author’s previous paper will be used, and the equations 
will be numbered consecutively with it; all equations 
numbered less than 10 or X appear in the first paper. 
In accordance with the above, the transverse force in the 
x-direction [Eq. (4f) ] is 


F 26 = w*Uzjm; (10a) 


Where u,; is the amplitude of vibration in the x-direc- 
tion at the point j, w is the frequency of vibration in 
radians per unit time, and m, is the mass per length 6. 
The transverse force in the y-direction is, similarly, 

F,j6 = w*uyym; (10b) 
From this it is found from Eq. (4f) that the column 
vector of bending moments due to inertia forces caus- 
ing vibration is 
Mz,—- My, 0 


My = wd u (10c) 

0 MmM,—™, 

where m, and m, are defined by Eqs. (4e)._ Now from 
this and Eq. (IV), it may be seen that 


u = 6(Q?m_g + w*m,) u (10d) 


M:;— My, 
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Variation of Natural frequency with rotating speed. 








Fic. 1. 


where mg, of this paper corresponds to m of the previous 


paper. By the same process as was used in finding 


Eq. (V), 


M, = ' 6[(Q?2m_g + w*m,) AJO-T, 1M, = (X&) 


U, = (w’C + @D)M, 


4 p 


(XA) 


which defines the 27 X 2n square matrices C and D. 
In using Eq. (X) to find the natural frequencies, given 
numerical values of &, it will be found that the calcula- 
tions for those vibrations that involve displacement pri- 
marily in the limber direction are made simpler by the 
fact that the coefficients for the stiff directions are so 
small as to be negligible. Thus, for these cases, the 
effective matrices reduce to n by n in size, reducing the 
work by a factor of 4.° 

As a practical matter, only the smallest values of w? 
are of interest; to find these, it is simpler to divide Eq. 


(X) through by ”. 


(1/w,2)M yp, = [C + (22/w,2)D]M,, (XB) 


where for a given value of Q, w, is one of the natural fre- 
quencies and \/,, is the corresponding column vector of 
moments. 

Then, for any value of the ratio 2?/w’, iteration will 
lead to the largest values of 1/w? or the smallest of w?. 
Some of these values may be negative, and these are of 
no physical significance; they may be eliminated by 
the usual trick of adding a corresponding constant to 
the main diagonal terms of [C + (Q?/w?)D]. If 


Qu = hu 
(0+ gk)u = (A+ k)u 


and all negative roots may be made smaller in absolute 
value than any positive roots by a suitable choice of k. 
To find the next lower value of 1/w?, it is necessary to 
use the orthogonality condition 
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XpsM pr = M>,'I-"M yy, = 0 


which will be proved later, Eq. (XIIB).® ' 


2 


w? — Q* contours 


In certain cases it is sufficient to know the natural 
frequencies at a given rotational speed, in which case 
the above analysis is sufficient; but often it is desired 
to know the functional relationship between w’? and 
Q?.6 From physical considerations it is possible to 
show that this relationship is a series of curves or con- 
tours such as are shown in Fig. 1. It should be empha- 
sized that these curves are purely qualitative. 

The requirement for a solution of Eq. (XA) is that 
the determinant be zero or 


lw2C + 2D — g| = 0 (lla) 


where 9 is the appropriate unit matrix. This is a 
(2n)th order polynomial in w? and 2? and their cross- 
products. It is not difficult to show that Eq. (lla) or, 
alternatively, Eq. (XA) isa positive definite form, so that 
for any positive value of Q® this must have 2n positive 
roots of w’, where one or more of the roots may be multi- 
ple. The multiple roots correspond to the points where 
the curves of Fig. 1 cross; thus, there are a total of 2” 
such curves. Physical reasoning shows that J), JJ,, and 
IIT, are the curves for the first three modes of vibration 
in the limber direction and J, and JJ, are the first two 
modes in the stiff direction. 

While it would be possible to calculate the values of 
w* for a series of values of 2° and thus plot the curves, 
there are ways of making this process a little easier. A 
more accurate curve may be drawn through a given 
number of points if the slopes are also known at these 
points. This slope is Ow*?/0?; in finding it, it is neces- 
sary to deal with the row vector, as well as the column 


vector, M,—i.e., 


M or = (w,?C + 2°D)M,, (11b) 
defines the column vector and 
aor = X pr(w,?C + 2?D) (lle) 


defines the row vector corresponding to a certain nat- 
ural frequency w, for a given 2°. Differentiating Eq. 


(11b) with respect to 2°, 


OM or Ow,” 7 e ° OM », 
se (= C+ p) My + (w?C + 2D) 
(11d) 
Premultiplying by x>,, 
, r OM >, 
[xor — Xpr(w2C + 2D)] = © 
on 
or (= C+ D) M >, (lle) 
The left-hand side is zero by Eq. (llc) and 
0w,?/02? = XprDM pr XprCM », (XI) 
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The numerator and denominator of the right-hand side 


are both scalors. 
Xp may be found directly from ,. 
transposed of Eq. (11b), 


M pr’ = Mor'(w,?C + 27D)! 
But, Eqs. (X) and (XA), 
(w?C + QD) = (63/E)O(w2m, + Q?mg) Ad-'T—} 


Forming the 


(12a) 


or 
(w?C + 2?D)’ = (68/E)I—-'0-"' [(w?m,, + Q?mg) A]’0’ 
(12b) 
however, 
I- = J-! (12c) 
thus, 
[(w?m,, + Q?mg) A]! = [(w*m, + Omg) A] (12d) 


since these are symmetrical matrices; and, Eq. (3c), 


or 


(w*C + QD)’ = (63/E)I—'0 X 


[(w?m,, + 22mg) AlO-! =I-"(w?C + OD) (12e) 
Substituting from Eq. (12e) into Eq. (12a) 
M,,'I- = M,,’I-*(w*C + 27D) (12f) 
By a comparison of Eq. (10f) with Eq. (9c), 
te = M,,'T-*, M pr’ = Xprl = (XITA) 


It may be shown that [reference 10, p. 77, Eq. (11)] 


Xpell or = 0 (r # Ss) (12g) 


where x,s iS the row vector corresponding to the 
natural frequency w,” and M,, is the column vector cor- 


responding to w,”; or 


M ysl (XIIB) 


IM,, = 0 


ASYMPTOTES 


Since the equation for w? vs. 2? is a (2n)th order 
polynomial leading to 2m curves, these curves must 
be asymptotic to straight lines for large values of Q?. 
Since these asymptotes are rather easily found, they 
make the plotting of the curves simpler. It may be 
shown that the curves of the actual propeller will be 
asymptotic to lines that have the same slope as our 
finite solution but somewhat different intercepts; how- 
ever, this does not lessen their value for practical plot- 
ting purposes. It would be possible to find the slopes 
of these asymptotes by evaluating Eq. (XI) for large 
values of 2°, but there is an easier way. The asymp- 
totes have the equations 


w,? = b,Q? + k, 


From the definition of an asymptote, these equations 


(13a) 
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must be approached by the polynomial as 2? « 
Substituting this expression into Eq. (11c), 


Xpr = Xpr[22(b,C + D) + k,C] (13) 
As Q?— o, this becomes 
Xpr(b,C + D) > 0 (13¢) 
—xprb,C = xpD 
—XxXprCD— = (1/b,)x pr (Q2?—> w~) = (134) 
Letting [Eq. (IIIA) ] 
‘i= a, xX = x,0 (13e) 


Then, from Eqs. (lle), (X), and (XA), 
Xr = x,(6°/E)(w,2m,, + Q?mg) AO-T-6 (13h) 


So that Eq. (13d) becomes 


—X,M,mMg' = (1/0,)x, (XIIIA 
mM, | is rather easily found since m, is a triangular 
matrix, Eq. (4e). Since 
m,—* 0 
z 
mg”! = 
—1 
0 mM, 
then 
(§ — mym,—") 0 
m,Mg ' = 
0 (m.m,—' —8) 


m, ‘ may be found from Eq. (4e) by normal methods 


and gives 


(m,—)yz = 0 (j < 2) 
l ; ‘ 
“> = (j = 4) 
km, 
PFs (13g) 
—1 l at 
sa n + n (7 > 1) 
dDkmy, Ddokm, 
(2,7 = 0,1,2,. . ,2—1) 


and m,mg~! is found by the usual matrix multiplica- 
tion. An additional advantage of using x, instead of 
Xpr arises from the fact that, since these refer only to 
large values of 2? and the beam stiffness is inoperative, 
the only possible directions of vibration are in the plane 
of rotation or out of the plane of rotation. Thus, let- 
ting 
Xr = bXzr Xv} 


then either x,, = 0 or x,, = 0, and Eq. (XIIIA) be- 
comes 
') = (1/by) x21 


—XerlF — MyM, 


(XIIIB 
—Xyr(mzm,—! — 8) = (1 — 
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FREE AND FORCED 


It may be seen that since 


i+ g= 1— g) 


—(§ — mym,—"') —(m,m, 


then 


and 


by, = ber +1 (13h) 


x, is found from Eq. (XIIIB), x», may be found from 
Bq. (13e), and .7,, may be found from Eq. (XITA). 

The above give the values of b,, x»,, and M,, which 
are used to find the k,. The determinant of the square 
matrix in Eq. (13b) must be zero, or 


07(b,C + D) + k,C — 9) = 0 (14a) 

Expanding this in powers of Q°, 
Ho” + H, — 0?" ~ +H, "-? +...4+H=0 
(14b) 


Itis easily seen that 


H, = \m,C + D| 


But this is zero by Eq. (12c) for the selected values },. 
Thus, Eq. (14b) may be written 


IT, , + 0(2-2) = 0 


i. =%;, oF oe (l4c) 
But 


1 
(14d) 


I, ] 


Oj 


nN l n nN 1 
Saat om Gi + k, : >» CijG ji = ( 
7=0 1 j3=0 
where G;; is the element in the jth row of the 7th column 
of the adjoint of (6,C + D) as defined in reference 10, 
page21. However, 


Gi; — M ix pv) (l4e) 


(reference 10, page 61, Eq. (3.5D)]. Thus, substituting 


from Eq. (14c) into Eq. (14d), 
Re = XorM or/XorCM oy (XIV) 
Arnoldi'!? has worked out a similar method to find the 
elect of small changes in blade angle on the natural fre- 


quency. The final equation is 
Ow,?/OB = —-2xprIM o:/XorCM », (XIVA) 
where 8 is the pitch angle and 
0 J 
J= 
—g 0 


It will be noted that the denominator is already used in 
Eq. (XIV) and that the numerator involves the scalor 
product of rows and columns only. To be explicit, 


f ' j 
.Xar Xowris _ 1 Xur Xerrf 
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FORCED VIBRATIONS 


Another important problem is to find the response of 
a rotating blade to aerodynamic forces that vary sinu- 
soidally with time. If the intensity of these forces per 
unit length of blade is known, it is not difficult to calcu- 
late the variable bending moments in the x and y direc- 


tions. Letting 
M (2, t) = M,0(z) cos wt | 
(15a) 
M (2, t) = M wo(z) cos wot 9 
then, in matrix notation, as before, 
M . 
{Maoh (15b) 


_ UMS 


where J, is the bending moment due to the amplitude 
of the applied vibratory forces only. Then it may be 
seen that [Eq. (XA) ] 


M, = (wC + 2°D)M, + 0Mo (XV) 
where w is now the frequency of the applied forces. 
This may be solved by iteration—i.e., assume a value 
of \/,, substitute in the right hand side, and take this 
solution as a new assumption for \/,. Unfortunately, 
this simple process will not always converge but 
may be made to converge by an averaging technique.'* 
Eq. (XV) gives only the vibratory moments; the 
may be found by a similar 


steady-state values 


technique.° 


A practical example of the above is given by an air- 
plane whose propeller axis makes an angle of yaw with 
the direction of the wind velocity. Because of this 
angle, the lift on the propeller varies with a frequency 
equal to the frequency of rotation; the most difficult 
job is to find Mo. An extensive discussion of how to 
find the variable lift on the propeller due to the yaw 
angle above is found in reference 13, but it may be 
simply stated as finding the lift, assuming a succession 
of steady-state lift conditions. It is necessary to con- 
sider the change in lift caused by twist of the blade, 
since the aerodynamic forces act at approximately the 
quarter chord point and not at the elastic axis; this 
twist effect can change the total vibratory lift by as 
much as 25 per cent for a normal propeller in use today. 
A method of analyzing this condition is demonstrated 
in reference 13, and a numerical solution to the forced 
vibration problem is also worked out and compared with 
experimentally determined values. The difference be- 
tween the calculated value and the experimental value 
is 6 per cent. A total of nine cases calculated from 
these principles by a slightly different procedure shows 
an error of less than 20 per cent in six of them and 
errors between 20 per cent and 31 per cent in three 
cases where the aerodynamic analysis, assuming steady- 
state conditions in a free airflow, is open to considerable 


question. 
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On the Theory of Heat Transfer Through a 
Laminar Boundary Layer 


Arthur N. Tifford 

Department of Aeronautical Engineering, The Ohio State 
University, Columbus 

January 9, 1951 


INTRODUCTION 


WwW: HAVE BEEN VERY MUCH INTERESTED in M. J. Lighthill's 
recent paper,' because of its close connection with some 
of the work supported here at The Ohio State University by the 
Office of Air Research of the Air Materiel Command, U.S.A.F. 
The present note is concerned with a relatively simple alteration 
of a portion of the analysis bringing about remarkable agreement 
with more exact analyses of the rate of heat transfer in the pres- 
ence of a surface pressure gradient. 

By way of review: in the analysis referenced, Lighthill as 
sumes, in addition to the usual assumptions, a linear velocity 
profile, « = [7o(x)/u]y. After an interesting application of oper 
ational calculus, he arrives at the following: 


l/s / 
= p ) V ato(x) ° 
Qu? (1/3)! 


Qix) = 
x me —1/3 
V o7(z) dz dT(x,) (29) 
0 x 
and 
l 1 
k/2 3p\ 
Qo( x) dx = = = 
0 (1/3)! \ wp? 
1 l ; 
V oro) dz dT(x) (380) 
wt @ x 
which simplifies, for isothermal surfaces, to 
1 l 
- Pp ‘ 
= ().807 ( ) V ori(x) dx (34) 
Me 0 


The comparison of the results obtained by means of Eq. (34) in 
the case of isothermal wedges, with the values determined by de- 
tailed numerical integration of the differential equations, has been 
given in his Table I for « = 0.7. 


TABLE I. (hl 


m — 0.0904 — 0.0654 0 
Lighthill’s values 0 0.497 0.601 
“Correct’’ values 0.438 0.541 0.585 
Error, per cent — 100 —8 +3 

TABLE II. (hi/k)(v 
m - 0.0904 — 0.0654 
Values from revised analysis 0.437 0.548 
“Correct’’ values 0.438 0.541 
Error, per cent —0.2 +1.3 
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ANALYSIS 


The error percentages shown in Table I arise from two distinct 
sources: the distortions of the boundary-layer temperature field 
close to the surface caused by (1) the presence of a surface pres- 
sure gradient and (2) the neglected boundary-layer velocity pro- 
file at some distance from the surface. Both distortions may be 
expected to vary with the Prandtl Number. 

The larger discrepancies tabulated are primarily caused by the 
presence of a surface pressure gradient. We shall therefore first 
try to include its effect in the analysis. 

The velocity profile in any laminar boundary layer is given by 


Eq. (1), 


; ot » & 
ee. te ee (1) 


) 
dx Qu p dx ~ 


To accomplish our purpose, we must retain the second, as well as 
the first, term of the infinite series in the approximate representa- 
tion of the velocity profile. Unfortunately, however, the analysis 
under these circumstances is extremely difficult, except in the 
i.c., dp/dx = O—and in the case of the 
i.e., t = 0. For the present, there- 


case of the flat plate 
separating boundary layer 
fore, let us seek an effective distance, 5, which will convert the 
parabolic velocity profile into the particular linear velocity profile 
having the same rate of heat transfer. 


effective u = (y/u)[ro + (dp/dx)é| (2) 


The effective distance is found to be approximately four-thirds of 
the momentum thickness of the boundary layer [6 = (4/3)6**]. 

As matters stand at this point, the rate of heat transfer would 
be proportional to o'/* regardless of the pressure gradient. How- 
ever, as has been shown in reference 2, although a one-third 
power law is appropriate to a linear boundary-layer velocity pro- 
file, a one-fourth power law is appropriate to a y? velocity profile. 
In order to cover all cases, the product, a7, should be replaced 
throughout Lighthill’s results by its effective value given by Eq. 
(3). 


effective (oro) = actual (ot») + (4/3)6**(dp/dx)o*/' (3) 


As for the second source of error, it has been found that in both 
the case of the flat plate and the case of the separating boundary 
layer,t percentage errors of +3 per cent (o = 0.7) are caused. 
This error may be separated into two parts—namely, a +2 per 
cent error occurring when the Prandtl Number is unity and a 
slight Prandtl Number correction associated with the magnitude 
of the Prandtl Numbers of gases (their values lie between 0.5 
and 2). The error may be eliminated by multiplying the results 
obtained after use of Eq. (3) by the factor, 0.980°-°?. 

Table II shows the result of the application of the revised 
analysis to the isothermal wedges considered in Table I. 

As can be seen, almost exact agreement with the ‘ 
values is obtained. Table III demonstrates that this agreement 
extends to the full range of Prandtl Numbers associated with 
gases in the two cases in the literature* in which the Prandtl 
Number effect has been studied in detail. 


‘correct”’ 


t Analyzed by Research Associate, John Wolansky 


k)\(v/ul)' 


1/9 1/3 1 4 
0.648 0.653 0.587 0.398 
0.596 0.576 0.495 0.325 
+9 +13 +184 4221/2 
ul)’. [¢ = 0.7] 

0 1/9 1/3 1 4 
0.585 0.591 0.571 0.495 0.329 
0.585 0.596 0.576 0.495 0.325 
0 —0.8 —0.8 0 +1.3 
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TABLE III. (hl/k)(v/m/)'? 


Flat Plate Forward Stagnation Region 


Values Values 

from from 

revised Pohlhausen’'s revised Squire's 
o analysis values analysis values 
0.6 0.554 0.552 0.466 0.466 
0.7 0.585 0.585 0.495 0.495 
0.8 0.613 0.614 0.522 0.521 
0.9 0.639 0.640 0.547 0.546 
1.0 0.664 0.664 0.570 0.570 
a3 0.686 0.687 0.592 0.592 


It is recommended that the revised analysis presented here be 
put to general use. To facilitate its application, attention is 
called to the graph of Fig. 1, developed by Falkner in reference 
4. This curve enables one to readily determine the approxi- 
mate momentum thickness of the laminar boundary layer at any 


point. 
SYMBOLS 
h = average heat transfer coefficient 
k = thermal conductivity of fluid 
l = total distance along surface from forward stagnation point 
m = constant exponent (uax"') 
p = static pressure 
Qo = local rate of heat transfer per unit area at the surface 
To = local temperature difference between the wall and the local free 
stream 
u = fluid velocity along surface 
ui = local free stream fluid velocity along surface 
x = distance along surface 
AY = distance perpendicular to surface 
2Z = dummy variable 
v) = effective distance 
6** = momentum thickness of boundary layer 
u» = absolute viscosity of fluid 
v = kinematic viscosity of fluid 
p = density of fluid 
o Prandtl Number of fluid 
TO local surface shearing stress 
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Calculation of Deflections in Indeterminate 


Structures 


H. F. Michielsen 
Chief of Structural Research, Royal Aircraft Factory Fokker, 
Amsterdam, Nederland 


February 1, 1951 


pa OF THE CONSIDERATIONS from Boley and Moore 
(Readers’ Forum, JOURNAL OF THE AERONAUTICAL Scr. 
ENCES, Vol. 17, No. 8, p. 526, August, 1950) and from Cox and 
Sopwith (Readers’ Forum, JOURNAL OF THE AERONAUTICAL 
ScIENCES, Vol. 17, No. 12, p. 815, December, 1950) gives rise to 
the impression that there must be some confusion regarding th 
essential features of both procedures under consideration. It js 
hoped that the following developments may somewhat clarify the 
situation. 


The internal loads L (stresses, bending moments, etc.) acting 
on the elements of a structure, can be written 


L = Lo + FL; l 


where Ly balances the known external loads, Lr is self-equilibrat 
ing (without external loads), and F is a statically indeterminate 


factor. If dS denotes the inverse stiffness of a structural element, 
F can be determined from Castigliano’s Theorem of Least Work, 
yielding 

(0/dF) f (L2/2)dS = 0 2 
hence 


S LLe dS = f[LoLey dS + Ff Ly? dS = 0 3 
From Eq. (3) the redundant quantity F is determined. 


Now, if the component of the deflection at a given point ina 
given direction is to be determined, the dummy load P’ = 1 is 
assumed to work in that point and in that direction. The dummy 
load would yield internal loads 


L’ = Lo’ + F'Lr | 
where again 
S Lo'Lr dS + F' f Le? dS =0 5 


Making use of the principle of virtual displacements, the de- 
flection 6 is then calculated from 


6= fLL'ds ‘ 


The purpose of the note of Boley and Moore now seems to kx 
to draw attention to the fact that Eq. (6) can be written as 


6= [ LL,’ ds 


which equation obviates the computation of F’ from Eq. (5 
and the determination of the integrals rit’ dS, as far as struc- 
tural elements not stressed by Lo’ are concerned. Especially i 
the case of highly redundant and complicated structures (for 
instance, wings and fuselages), this seems to be an extremely im- 
portant advantage. The load Ly’ can be chosen in any suitable 
way (involving the least number of elements) able to equilibrate 


the dummy load. 


Eq. (7) can be derived from Eq. (6), making use of the con- 
sideration that, according to Eq. (4) and the outer members 0 
Eq. (3), 


SLL’ dS — [Ll dS = F' f LLy dS = 0 8 
In exactly the same way it can be proved that 
6 = [LoL' dS 9 


which equation obviates the computation of F from Eq. (3) an¢ 


can be important if, in a structure of known cross-sectional di- 
mensions, the deflection in a single point at several different 
loadings Ly is to be determined. 
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READERS’ 


Cox and Sopwith, on the other hand, set up the general ex- 
pression for the strain energy U, which can be written, with the 


notations used here, 
U = J Lot dS + 2P’ f Lolo’ dS + 2F [Loker dS + 
Pr? [Lo 2dS + 2P'F[ Ly'Lr dS + P°L'Le? dS (10) 


Equating 0U/OF to zero with P’ = 0 leads to 


S LoLe dS + Ff[/ Lr? dS = 0 (11) 
in accordance with Eq. (3). 
The deflection is then determined from 6 = OU/OP’, with 
Pp’ = 0, 
6 = [Lolo dS + FS Lo'Lr dS (12) 


and, because of Eq. (1), Eqs. (12) and (7) are identical. Never- 
theless, Eq. (7) seems to be preferable to Eq. (12) for practical 
computation. Its form is extremely simple and easy to keep in 
mind, and it attacks more directly the problem involved—for 
instance, if stresses are known but no complete analysis is avail- 
able. 

Unfortunately Messrs. €ox and Sopwith do not mention the 
texts in which the procedure, outlined in their note, can be 
found. Some well-known textbooks deal with deflections of 
determinate structures and analysis of redundant structures but 
without drawing special attention to deflection of redundant 
structures. 

Although the method, described by Messrs. Boley and Moore, 
has been practiced in our stress offices for many years, our ex- 
perience is that, in its simple form as expressed in Eqs. (7) and 
(9), it is not as well known as it ought to be. We therefore 
certainly agree with the purport of the first paragraph of their 


note. 


The Calculation of Deflections of 
Indeterminate Structures 


Bruno A. Boley and Roland H. Moore 
The Ohio State University, Columbus, and Goodyear Aircraft 

Corporation, Akron, Respectively 
January 30, 1951 
— OF CALCULATING DEFLECTIONS of indeterminate 

structures was recently discussed in reference 1; 
questions on it have been brovght up by H. L. Cox and D. G. 
Sopwith in reference 2. The treatment of reference 2 differs 
from that of reference 1 principally in that it makes use of Castig- 
liano’s theorem rather than of the principle of virtual displace- 
For any linearly elastic structure, the two methods are 
equivalent; the main difference is one of procedure. When 
Castigliano’s theorem is used, it is customary to include the 
dummy loads with the actual loading and the redundant quanti- 
ties at the outset in calculating the strain energy U; the desired 
deflection is then the partial derivative of U with respect to the 
appropriate load. When the principle of virtual displacements is 
employed, the redundant quantities are first evaluated in terms 
of the actual loading; then, as a separate step, the dummy 
system is applied to the structure, the internal stresses due to it 
are found, and the deflection is obtained by the equation of refer- 
ence 1. 

It was stated in reference 1 that in the calculation of the in- 
ternal stresses due to the dummy load, it is permissible to disre- 
gard the redundant members of the structure. This simplifica- 
tion is used automatically when the calculations are performed 
by Castigliano’s theorem (because of the use of the partial de- 
rivative); not so when the principle of virtual displacements is 
The latter procedure is more com- 


some 


ments. 


applied as discussed above. 
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monly used in the United States, and for this reason a special 
mention of the subject was thought desirable. It might be men- 
tioned that, if, in using the Castigliano’s theorem approach, the 
calculations were (for convenience) to be split in two parts, in the 
second part (which would consist of evaluating the total energy 
stored due to the action of the dummy loads) the redundant mem- 
bers could again be neglected. 

The authors appreciate the criticisms of Messrs. Cox and Sop- 
with. They also wish to acknowledge the comments of G. A. 
Starr, of Chance Vought Aircraft, who stated that reference 3 
contains applications (related to beam problems) of the simplifi- 
cation here considered, and the remarks of Dr. N. J. Hoff, of the 
Polytechnic Institute of Brooklyn, indicated that the 
present procedure appears in his mimeographed notes, soon to be 
published under the title of Minimal Principles of Structural 


Analysis. 


who 
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On Plastic Buckling of Plates and a Theory of 
Plastic Slip 


P. Cicala . . 
Instituto Aerotecnico and University of Cordoba, Argentina 


January 8, 1951 


N A READER’S ForuM NOTE (JOURNAL OF THE AERONAUTICAL 

Scrences, Vol. 17, No. 6, p. 378, June, 1950), the writer 
showed that the predictions of a flow theory of plasticity can 
agree with test results on buckling of a compressed sheet if cer- 
tain initial eccentricities are supposed to exist in the test speci- 
mens. In an interesting discussion (JOURNAL OF THE AERONAU- 
TICAL ScrENCES, Vol. 17, No. 11, p. 742, November, 1950), Pro- 
fessor Bijlaard argued that these eccentricities are relatively 
large and that the tested sheets probably do not show such a 
degree of irregularity. If this is the case, their low buckling loads 
give experimental evidence against the Prandtl-Reuss theory. 
As other tests are in conflict with deformation theory, it may be 
inferred that neither theory can be confidently used in the analy- 
sis of plasticity problems involving marked changes in the direc- 
On the contrary, a new ‘“‘plastic slip” 
Some con- 


tion of the stress tensor. 
theory! appears to agree with known test results. 
clusions deduced® * from this analytically rather complex theory 
are presented here. 

If a state of uniaxial stress 6, = o><a, (x, y, z orthogonal axes, 
¢» proportional limit) is reached through a continued increase of 
only this component and, subsequently, a small change 62 in the 
stress tensor is produced, the following cases are to be dis- 
tinguished: 

(a) If 6E contains a sufficiently large increment 60 of o,, the 
strain increments are related to the stress increments by equa- 
tions of the same form as for previous theories. The coefficients 
appearing in these equations depend only upon the initial stress 
ao; it was shown? that they can be deduced from the simple ten- 
sion test without requiring the calculation of the characteristic 
function F, which in this theory represents the unit plastic slip. 
For instance, denoting by 1/E, the coefficient of 5¢, in the expres- 
sion of the strain increment de,, we find 


o 
1/E, = (3/8E,) + (1/4E,) + (3 86) de/E, (1) 
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where £; and &, are the tangent and secant moduli. If only the 
components og, and go, are incremented, Eq. (1) holds when 
50/50,>0/o, if 0, has the same sign as éo and o or when — do + 
60:>(a/cp) — 1 if 50, has the opposite sign. 

(b) If 52 contains a sufficiently large decrease of o,, all plastic 
strain components disappear. 

(c) If d¢ lies between certain limits, the coefficients of the men- 
tioned equations depend not only on the initial stress o but also 
upon the ratios of stress increments; their calculation implies 
the function F, whereas the limiting values do not depend on F. 
Reversing all the signs of the stress increments which correspond 
to a limiting configuration between cases (a) and (c) leads to a 
limiting configuration between cases (b) and (c). 

In the particular case that only a tangential stress 67 in a plane 
perpendicular to z is introduced, the relation between 60, 67 and 
the corresponding shear strain 6y may be summarized in a graph 
(see Fig. 1). The horizontal lines through A and B correspond to 
conditions (a) and (b), respectively. The curve AMB was ob- 
tained for 0 = 1.25¢,, assuming for F the simplest linear expres- 
sion. 

The buckling of a strip of thickness s, width 6, length L>»d, 
with one of its long sides free and the other simply supported, 
subjected to longitudinal compression «, may be analyzed ac- 
cording to the above results. The case of vanishing initial 
eccentricities is considered. We suppose that the values of the 
incremental ratios (Fig. 1) calculated for the conditions at onset 
of buckling can be used to represent the corresponding deriva- 
tives during the whole buckling process. Suppose that the ratio 
t/y takes the same values for all the points of a transversal sec- 
tion; then we obtain? 


1 \ 9 
ao = (s*v/b*u) v au 
0 ~ 


y = (su/b)? — o(u/v)? 


whereu = \ ‘dr/do and v = V dr ‘dy. By means of the diagram 
in Fig. 1, vcan be deduced as a function of u, and, thus, from Eqs. 
(2), the relationship between the axial load and the strain y 
(which we may measure at the most stressed fiber) is obtained in 
a parametrical form. 

For a given a, Eqs. (2) yield the same value of y and, hence, of 
the unit twist, for any section. Thus, neglecting the bending 
and warping stresses leads to a triangular variation of the angle of 
twist, the maximum being attained halfway between the sup- 
ported end sections. 

With fair accuracy, for «>, we may write 


v/VG, =1 +k (u2 — uo?) "u-™ (3) 
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where uw = 1/ V mo is the value corresponding to point A, G, and 
Gm are the moduli defined in Fig. 1, ” is a positive integer that 
may be given the value = 3, and k = VW Gm 'G, — 1. 

With Eq. (3) andv = / G, for u < uo, Eqs. (2) can be readily 
evaluated. For small values of y, the solution may be expressed 
in the form 


o/os = 1 + (-yb?/s2u?) — [kn/(n + 1)] (yb?2/ks2uo2) "+ 1)/n 
+ higher powers of y 


where o, = G,s?/b? is the stress at which buckling starts, coincid- 
ing with the known critical stress deduced on the basis of defor- 
mation theories of plasticity. 

By increasing y, the stress o tends to the value G,,s?/b?, which, 
here, takes the same role as the reduced modulus stress in column 
buckling. 

The approximation may be improved by taking into account 
the variation of r/y over the section. Denoting, as before, by y 
the maximum shear strain and by defining 7 as the ratio of the 
twisting moment to bs?/3, Eqs. (2) hold unchanged, whereas a 
new relation v(u) must be found. According to calculations 
carried out for s<b and k<1, the same Eq. (3) with m = 4 may 
be used with sufficient accuracy. 

The above results evidently need experimental checks. If 
differential relations between twist, torque, and axial load similar 
to the one between y, 7, and o plotted in Fig. 1 are obtained 
experimentally, the torsional buckling in the elastoplastic range 
can be described by the above equations in cases where the effect 
of warping stresses is negligible. 

The solution presented, though roughly approximate, gives a 
qualitative outlook of the buckling process of plates, evidencing 
its analogy to the Shanley phenomenon. 
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The Effect of Prandtl Number on the 
Theoretical Shock-Wave Thickness 


Paul A. Libby 

Assistant Professor of Aeronautical Engineering, Polytechnic 
Institute of Brooklyn 

January 10, 1951 


T° REFERENCE 1, L. Meyerhoff has presented the results of an 
extensive investigation of the theoretical structure of a one- 
dimensional shock wave according to the usual hydrodynamic 
approximations. Among other interesting results he has indi- 
cated that, as the Prandtl Number (NV ,) is changed from 0.75 to 
0.657 in the equations for constant gas properties [Eq. (8) or 
(8c) of reference 1], the shock-wave thickness is considerably 
increased. Since for air, the Np, is indeed somewhat less than 
3/4, this result would indicate that the hydrodynamic equations 
might describe the structure of a shock wave in air with more ac- 
curacy than is usually attributed to them on the basis of Becker’s 
exact solution for Np, = 3/;. However, this accuracy is usually 
gaged by comparing the shock-wave thickness as defined by 
Prandtl or by Taylor and Maccoll with the mean free path of the 
gas molecules. The comparison in reference 1 between the struc- 
ture with Np, = 0.657 and Np, = */, has been made with differ- 
ent conditions upstream of the wave and, therefore, with differ- 
ent mean free paths. It is of interest, therefore, to replot Fig. 2 
of reference 1 in terms of the mean free path of the gas molecules. 
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For convenience, the mean free path of the gas molecules far 
ahead of the wave is used as the reference mean free path. For 
the initial conditions given in Fig. 2 of reference 1, one finds, for 
the Np, = 0.657 solution, in the usual notation 

po = 0.457 X 10-* Gm./cm.’ } 

To = 500°C. . (1) 

wo = 3.56 X 10 ‘Gm. /cm.-sec.*| 
where the subscript zero denotes conditions far ahead of the 
wave. The mean free path ) is defined by the equation 

A = pw/dpc (2) 

where 6 = a constant depending on the accuracy of the theoretical 
analysis (1/3 < 6 <0.499) (the value 0.35 due to Tait is used 
here) and where c = mean molecular velocity, V/ (8 a) RT. 

From the values given in Eq. (1), 


rN = 2.96 X 10-5 cm. (3) 
For the Becker solution (Np, = */;) as shown in Fig. 2 of refer- 
ence 1, one finds that 
po = 0.752. 107-3 Gm./cm.3 } 
Ty = 0°C. ‘ (4) 
Mo = 1.72 XK 107-4 Gm. /om.-sec.*) 
ind that 
dh = 1.46 X 10-5 cm. (5) 


If the two curves are plotted in terms of — = x/Xo, one obtains 
the result shown in Fig. 1 herein. It will be seen that, except for 
the region downstream of the wave where the computational 
instability discussed by Meyerhoff causes his numerical solution 
to deviate from the Rankine-Hugoniot asymptote, the computed 
points fall on a single curve. From this result (which, of course, 
could not be obtained without the solution for Np, = 0.657), it 
may be concluded: (1) that the increase in shock-wave thickness 
with Np, shown in reference 1 is really due to the change in initial 
conditions; (2) that the analyses based on Np, = 4/, are valid for 


* The values of uw have been taken from the source listed in reference 1 
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Fic. 1. Shock-wave structure for two Prandtl Numbers. 


other values of Np, of interest for air; and (3) that the doubts 
that arise concerning the accuracy of the hydrodynamic equa- 
tions for the description of the structure of a shock wave at high 
Mach Numbers and which are based on the exact analysis with 
Np, = 3/4should be considered to remain for the actual N p, values 
exhibited by air. The discussion here is based on the solution 
for constant gas properties, but it can be expected to apply for 
the case of variable gas properties, variable Vp,, as well. 
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